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implanted  with  ‘‘He  and  ^Be  have  been  studied  using  PL  and  SIMS. 

Following  60  min,  750°C  anneals,  indium  is  detected  on  the  surfaces  of  SI3N4 
encapsulant  layers  by  both  AES  and  SIMS.  SIMS  results  also  suggest  that  indium 
may  outdiffuse  through  PSG  caps  during  similar  anneals.  Indiffusion  of  Si  from 
the  cap  to  the  underlying  InP  is  found  to  be  significant  in  annealed  Si02  and  S^Np 
capped  samples.  Little  such  Si  contamination  is  observed  after  annealing  InP  with 
PSG  encapsulants.  PL  measurements  indicate  that  the  spectral  features  and  lumino¬ 
sity  levels  of  capped  and  annealed  InP  are  least  altered  by  PSG  encapsulants. 

From  these  results  it  is  concluded  that  PSG  encapsulation  best  preserves  the 
initial  characteristics  of  encapsulated  InP  during  furnace  anneals. 

PL  measurements  indicate  that  following  51  keV,  1013-1015  cm  2  implants  of 
‘‘He,  thirty  minute  anneals  at  750°C  appear  sufficient  to  restore  the  lattice  to 
approximately  pre-implant  levels  of  quality.  For  30  min  anneals  of  100  keV,  9Be^ 
implanted  InP,  maximum  Be  activation  appears  to  occur  for  all  doses  (1013-1015cm  2 ) 
at  ~  750  C,  1015  cm  2  9Be  implanted  InP  exhibits  saturation  of  Be  activation  be¬ 
tween  30-60  min  with  decreasing  Be  associated  luminscence  resulting  from  longer 
anneals.  The  Be  acceptor  ionization  energy  is  estimated  to  be  41.3  +  3  meV  with 
characteristic  emission  at  1.382  eV  (5°K) . 

Depletion  zones  of  the  compensating  impurity,  Cr  or  Fe,  are  observed  with 
SIMS  in  SI  InP  capped  with  S^Np  and  annealed  at  800°C  for  30  min.  Cr  exhibits 
larger  concentration  gradients  in  these  zones  than  does  Fe. 

Under  typical  annealing  conditions  for  InP  (T  >  700°C,  t  =*  15-30  min)  it  is 
observed  using  SIMS  that  implanted  9Be  is  a  rapid  diffusant  in  SI  InP.  High  dose 
(1015  cm  2) ,  100  keV  implants  of  9Be  into  SI  InP  result  in  structured  and  corre¬ 
lated  Be  and  Cr  or  Fe  atomic  profiles.  In  addition,  flat  tails  of  9Be  extending 
over  distances  on  the  order  of  a  micron  are  observed  in  high  dose  implanted  InP. 
Such  constant-concentration  tails  of  9Be  are  not  observed  in  comparably  processed, 
undoped  vapor  phase  epitaxial  (VPE)  InP. 

Drastic  redistribution  of  Mg  and  Fe  occurs  during  30  min  anneals  of  SI  InP 
implanted  with  1015  cm  2,  250  keV  26Mg.  As  in  the  previous  Be  studies  using  SIMS, 
constant-concentration  tails  of  Mg  extending  over  distances  of  1  pm  or  more  are 
observed  in  the  bulk  following  annealing  and  the  redistribution  of  Fe  is  again 
correlated  to  the  diffusion  profile  of  the  implanted  acceptor. 

With  SIMS,  little  redistribution  of  28Si  occurs  for  30  min,  750°C  anneals 
of  5  x  1014  cm  2,  240  keV  28Si  implants  of  SI  InP.  This  behavior  of  the  implante 
donor  Si  differs  significantly  from  the  rapid  redistribution  behavior  of  the 
implanted  Be  and  Mg  acceptors  studied  here. 
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Successful  ion  implantation  of  InP  will  depend  on  detailed  Infor¬ 
mation  regarding  the  protective  qualities  of  prospective  encapsulants  and 
also  the  thermal  annealing  parameters  (i.e.,  time  and  temperature)  required 
to  reorder  the  implant  damaged  lattice  and  to  activate  the  implanted 

impurity.  Using  Auger  electron  spectroscopy  (AES) ,  secondary  ion  mass 

spectrometry  (SIMS),  and  low  temperature  photoluminescence  (PL)  the  InP 

encapsulation  properties  of  CVD  Si02»  CVD  PSG  (phosphosilicate  glass),  and 

rf  plasma-deposited  Si^  have  been  examined.  In  addition,  the  annealing/ 

implantation  characteristics  of  semi-insulating  (SI)  InP  implanted  with 
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He  and  Be  have  been  studied  using  PL  and  SIMS. 

Following  60  min,  750°C  anneals,  indium  is  detected  on  the  surfaces 
°f  Si3N4  encaP8ulant  layers  by  both  AES  and  SIMS.  SIMS  results  also  suggest 
that  Indium  may  outdiffuse  through  PSG  caps  during  similar  anneals.  Indif¬ 
fusion  of  Si  from  the  cap  to  the  underlying  InP  is  found  to  be  significant 
in  annealed  Si02  and  Si^  capped  samples.  Little  such  Si  contamination  is 
observed  after  annealing  InP  with  PSG  encapsulants.  PL  measurements  indi¬ 
cate  that  the  spectral  features  and  luminosity  levels  of  capped  and  annealed 
InP  are  least  altered  by  PSG  encapsulants.  From  these  results  it  is  con¬ 
cluded  that  PSG  encapsulation  best  preserves  the  initial  characteristics 
of  encapsulated  InP  during  furnace  anneals. 
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PL  measurements  Indicate  that  following  51  keV,  10  -10  cm 

implants  of  ^He,  thirty  minute  anneals  at  750°C  appear  sufficient  to  restore 

the  lattice  to  approximately  pre-implant  levels  of  quality.  For  30  min 
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anneals  of  100  keV,  Be  implanted  InP,  maximum  Be  activation  appears  to 
occur  for  all  doses  (10^-10^  cm  at  “  750®C.  At  750°C,  lO^cm  ^  ^Be 
implanted  InP  exhibits  -  saturation  of  Be  activation  between  30-60  min  with 
decreasing  Be  associated  lumlnscence  resulting  from  longer  anneals.  The  Be 
acceptor  ionization  energy  is  estimated  to  be  41.3  ±  3  meV  with  character¬ 
istic  emission  at  1.382  eV  (5°K). 

Depletion  zones  of  the  compensating  impurity*  Cr  or  Fe,  are 
observed  with  SIMS  in  SI  InP  capped  with  and  annealed  at  800°C  for 

30  min.  Cr  exhibits  larger  concentration  gradients  in  these  zones  than 
does  Fe. 

Under  typical  annealing  conditions  for  InP  (T  >  700°C,  t  -  15-30 
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min)  it  is  observed  using  SIMS  that  implanted  Be  is  a  rapid  diffusant 

in  SI  InP.  High  dose  (10^  cm  ^) ,  100  keV  implants  of  ^Be  into  SI  InP 

result  in  structured  and  correlated  Be  and  Cr  or  Fe  atomic  profiles.  In 
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addition,  flat  tails  of  Be  extending  over  distances  on  the  order  of  a 

micron  are  observed  in  high  dose  implanted  InP.  Such  constant-concentration 
9 

tails  of  Be  are  not  observed  in  comparably  processed,  undoped  vapor  phase 
epitaxial  (VPE)  InP. 

Drastic  redistribution  of  Mg  and  Fe  occurs  during  30  min  anneals 
15  —9  26 

of  SI  InP  implanted  with  10  cm  ,  250  keV  Mg.  As  in  the  previous  Be 
studies  using  SIMS,  constant-concentration  tails  of  Mg  extending  over 
distances  of  a  1  pm  or  more  are  observed  in  the  bulk  following  annealing 
and  the  redistribution  of  Fe  is  again  correlated  to  the  diffusion  profile 


of  the  Implanted  acceptor. 
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With  SIMS,  little  redistribution  of  Si  occurs  for  30  min,  750°C 
14  —2  28 

anneals  of  5  x  10  cm  ,  240  keV  Si  implants  of  SI  InP.  This  behavior 
of  the  implanted  donor  Si  differs  significantly  from  the  rapid  redistribu¬ 
tion  behavior  of  the  implanted  Be  and  Mg  acceptors  studied  here. 
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1.  INTRODUCTION 

In  the  last  several  years  the  semiconductor  community  has  witnessed 
the  emergence  of  InP  (and  its  related  alloys)  into  the  position  of  the  third 
most  significant  semiconductor  after  Si  and  GaAs.  This  has  occured  because 
of  potential  device  performance  benefits  offered  by  InP  over  either  Si  or 
GaAs  in  several  areas. 

One  area  in  which  InP  devices  are  expected  to  have  major  impact 
and  also  find  their  earliest  application  is  in  the  field  of  microwave  and 
millimeter  wave  solid  state  components  [1-4].  Much  of  the  motivation  of  the 
current  research  in  InP  originated  with  a  theoretical  proposal  in  1970  by 
Hilsum  and  Rees  [5]  of  a  "three-level"  electron  transfer  effect  in  InP  lead¬ 
ing  to  superior  microwave  characteristics  in  InP  Gunn  diodes  in  comparison 
to  similar  GaAs  devices.  Subsequent  experimental  results  refuted  the  exis¬ 
tence  of  a  "three-level"  effect  in  InP  [6,7];  however,  the  predictions  by 
Hilsum  and  Rees  were  an  i'-portant  catalyst  in  initiating  serious  research 
efforts  examining  the  device  potentials  of  InP.  As  a  result  of  such  efforts 
it  is  now  known  that  because  of  several  transport  and  material  properties 
(i.e.,  energy  relaxation  time,  peak-to-valley  ratio,  electron  diffusion  co¬ 
efficient  at  fields  above  threshold,  and  thermal  conductivity)  that  InP 
microwave  and  millimeter  wave  devices  should  operate  at  higher  frequencies, 
with  higher  efficiency,  and  with  lower  noise  than  comparable  GaAs  devices 
[1,4]. 

Another  area  in  which  InP  is  expected  to  make  significant  contribu¬ 
tions  is  in  the  field  of  electro-optics  [8].  Large  efforts  are  underway  at 
several  laboratories  to  develop  sources,  detectors,  and  even  integrated 
optoelectronic  circuits  for  the  1.1  -  1.6  ym  wavelength  region  where  state-of- 
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the-art  optical  fibers  exhibit  their  lowest  attenuation  and  dispersion 
[9,10].  In  contrast  to  Si  or  GaAs/AlGaAs  based  electro-optical  devices, 
InP/lnGaAsP  based  devices  easily  span  the  wavelength  region  from  0.9  -  1.7 
ym.  Besides  microwave  and  electro-optical  applications  there  is  also 
speculation  that  InP  may  become  important  in  the  fabrication  of  high  speed 
Integrated  logic  circuits  [11,12].  InP  is  an  attractive  material  for  such 
circuits  because  the  interface  state  density  in  InP  metal-insulator-semi¬ 
conductor  (MIS)  devices  is  intermediate  between  that  of  Si  and  GaAs  [13-15]. 
This  property  means  that  while  retaining  the  high  speed  characteristics  of 
InP  or  GaAs,  InP  integrated  logic  circuits  could  employ  normally-off  de¬ 
vices,  in  marked  contrast  to  the  present  limitations  of  metal-gate  active 
devices  in  GaAs  IC  technology. 

Despite  such  research  interest  in  the  technological  potentials  of 
InP,  the  material  and  devices  currently  fabricated  from  it  are  in  a  primi¬ 
tive  stage  of  development  even  in  comparison  to  the  still-evolving  GaAs 
technology.  The  main  factors  inhibiting  larger  research  efforts  into  InP 
devices  are  the  availability  and  quality  of  single  crystal  InP  and  the  lack 
of  a  processing  technology  appropriate  to  the  requirements  of  InP  device 
fabrication.  The  availability  of  commercial,  single  crystal  InP  Is  presently 
confined  to  a  small  number  of  vendors  with  the  majority  of  high  quality, 
high  purity  InP  still  grown  by  government  or  corporate  research  labs  for 
"in-house"  use.  Since  the  early  reports  of  liquid  encapsulated  Czochralski 
(LEC)  growth  of  InP  by  Mullin  et_  al_.  [16,17]  major  improvements  in  the 
quality  of  bulk,  single  crystal  InP  have  occured.  In  comparison  to  ordinary, 
single  crystal  Si  even  the  best  currently  available  InP  is  poor  in  terms  of 
purity  and  crystalline  defects.  Further  refinement  in  InP  growth  techniques 


and  a  lowering  of  Che  material  costs  await  the  economic  impetus  provided 
by  a  wider  acceptance  of  InP  by  the  device  community.  The  present  inertia 
in  developing  InP  devices  is  primarily  due  to  the  fact  that  only  a  limited 
body  of  work  exists  concerning  the  processing  characteristics  of  this  semi¬ 
conductor. 

Studies  presented  in  this  dissertation  address  issues  associated 
with  one  specific  area  of  InP  material  processing,  that  of  ion  implantation. 
It  is  envisioned  that  ion  implantation,  due  to  its  control  and  reproduci¬ 
bility,  will  be  used  in  the  manufacture  of  InP  devices  such  as  it  is  now 
routinely  employed  in  Si  or  GaAs  fabrication.  Compared  to  other  aspects  of 
the  newly  emerging  InP  technology,  research  concerning  InP  ion  implantation 
is  a  recent  development.  The  first  general  surveys  of  ion  implantation  in 
InP  were  published  by  Davies  jlJ..  [18]  and  Donnelley  et  al .  [19]  in 
1977.  In  the  years  since  these  first  reports  a  large  number  of  questions 
have  arisen.  Specifically,  a  few  such  questions  are:  Which  encapsulant 
best  protects  and  preserves  the  characteristics  of  encapsulated  InP  during 
furnace  anneals?  What  are  the  annealing  parameters  (i.e. ,  time  and  tempera¬ 
ture)  most  effective  in  reordering  the  implant-damaged  lattice  and  in 
activating  implanted  Impurities?  What  is  the  nature  of  atomic  redistribu¬ 
tion,  if  any,  exhibited  by  implanted  or  other  intentionally  present  impuri¬ 
ties  during  furnace  anneals?  Study  of  these  questions  forms  the  core  of 
this  thesis.  As  might  be  expected  in  experimental  research,  during  the  span 
of  this  work  additional  issues  and  questions  have  unexpectedly  been  encoun¬ 
tered  and  are  also  described  here.  The  results  presented  in  this  work  should 
contribute  to  the  current  understanding  of  ion  implantation  in  InP,  and 
thereby  encourage  a  broader  utilization  of  this  material. 
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2.  IMPLANTATION  AND  ANNEALING  PROCEDURES 

2.1  Material  Characteristics 

The  bulk  InP  (100)  substrates  used  in  this  work  were  grown  by  the 
liquid-encapsulation  Czochralski  (LEC)  method  at  the  Naval  Research  Labora¬ 
tory  (NRL)  [20],  Sn-doped  (100)  InP  (NRL  #2-13-H)  was  used  during  Auger 
electron  spectroscopy  (AES)  and  secondary  ion  mass  spectrometry  (SIMS) 
studies  of  the  annealing  encapsulants  to  minimize  charging  effects  during 

analysis  of  these  dielectrics.  The  Sn  concentration  in  this  n-type  material 
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was  determined  by  NRL  to  be  approximately  2  x  10  cm  .  The  photolumines¬ 
cence  measurements  were  made  on  semi-insulating  (SI) ,  Fe-doped  (100)  InP 
(NRL  //2-10-H) .  This  material  had  a  resistivity  of  ~  8  x  10"^  cm  ^  as  de¬ 
termined  by  SIMS  [21].  The  SIMS  studies  of  atomic  redistribution  during 
anneals  were  performed  with  both  Fe-  and  Cr-doped  InP  (100) .  The  Fe-doped 
material  for  these  redistribution  studies  was  the  same  as  that  employed  in 

the  photoluminescence  measurements.  The  SI,  Cr-doped  (100)  InP  (NRL  #l-48-H) 

3  52 

had  a  resistivity  of  ~  5  x  10  ficm  and  a  Cr  atomic  concentration  of 
"  6.5  x  10"^  cm  ■*,  as  determined  with  SIMS  [21]. 

A  detailed  discussion  of  the  LEC  growth  procedure,  together  with 
an  examination  of  boron  contamination  in  these  crystals,  is  presented  in 
Chapter  10. 

2.2  Implantation 

Ion  Implantation  is  now  a  standard  technique  for  doping  Si  and 
GaAs,  and  the  general  properties  of  the  method  have  been  extensively  re¬ 
viewed  in  the  literature  [22-27].  Therefore,  the  remarks  in  this  section 
will  be  confined  to  experimental  details  most  relevant  to  this  work. 
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The  Implantation  process  originates  with  a  plasma  containing  a 
number  of  ionized  impurities.  These  ions  are  then  accelerated  to  high  energy 
and  the  desired  impurity  is  selected  from  the  others  using  mass  separation 
(post-acceleration  analysis) .  This  desired  ion  beam  is  then  raster  scanned 
across  the  target.  The  total  impurity  dose  is  obtained  by  monitoring  the 
current  flowing  from  the  target  or  surrounding  Faraday  cups  to  ground.  Upon 
penetration  of  the  target  the  ions  lose  their  energy  via  three  processes: 

(i)  collisions  between  the  ions  and  the  screened  nuclei  of  the  target  atoms 
(nuclear);  (ii)  interactions  with  the  bound  electrons  of  the  target  (elec¬ 
tronic);  (ill)  charge  exchange  between  ions  and  host  atoms  (exchange)  [28, 

29].  It  is  the  nuclear  collision  process  which,  by  displacing  the  target 
atoms  from  their  lattice  sites,  produces  various  forms  of  radiation  defects 
that  must  be  removed  by  annealing.  Following  implantation,  the  "as-implanted" 
impurity  distribution  is  a  function  of  the  ion's  mass,  the  implantation 
energy,  target  orientation,  and  the  material  characteristics  of  the  target. 

The  basic  features  of  the  implanted  profile  such  as  mean  range  and  standard 
deviation  can  be  fairly  accurately  predicted  from  the  theory  of  atomic  stop¬ 
ping  developed  by  Lindhard,  Scharff,  and  Schiott  (hereafter  referred  to  as 
LSS) [30-33].  Usually  the  implant  profile  can  be  adequately  described  by  a 
gaussian  distribution  using  the  mean  range  and  standard  deviation  obtained 
from  LSS  [22,24].  It  should  be  noted,  however,  that  certain  combinations  of 
projectile,  target,  and  implant  energy  can  result  in  skewed  implant  profiles 
which  can  be  accurately  represented  only  by  using  higher  moments  of  the 
distribution.  Techniques  based  on  LSS  theory  for  calculating  these  higher 
moments  are  described  in  the  literature  [34-37].  The  use  of  such  higher 
moments  as  it  pertains  to  this  work  will  be  discussed  in  Chapter  8. 


Figure  2.1:  Schematic  diagram  of  the  300  keV,  A. I.  300-MP  ion  implanter. 
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Implantations  for  this  work  were  performed  with  a  300  keV  Accel¬ 
erators  Inc.,  300-MP  Ion  implanter  shown  schematically  In  Figure  2.1.  All 
Implantations  were  done  at  room  temperature  with  the  sample  normal  Inclined 
7°  from  the  incident  beam  to  minimize  channelling  effects.  The  schedules 
for  these  implants  are  given  in  Table  2.1  together  with  relevant  range  and 
damage  statistics.  Values  for  the  theoretical  atomic  profiles  were  obtained 
from  a  computer  program  utilizing  the  Gibbons,  Johnson,  and  Mylroie  formula¬ 
tion  [37]  of  LSS  theory  (full  program  output  for  these  implants  is  given  in 
Appendix  1) .  The  statistics  for  the  damage  profiles  have  been  obtained  from 
Winterbon's  tables  of  data  [35].  It  is  assumed  that  the  profile  of  energy 
deposited  into  atomic  collisions  is  directly  proportional  to  the  amount  of 
damage  (i.e. ,  Frenkel  defects,  vacancies,  and  divacancies)  produced.  The 
values  for  the  energy  deposition  density  listed  in  Table  2.1  are  approxima¬ 
tions  based  on  the  product  of  the  dose  and  the  nuclear  energy  loss  per  unit 
length  [37].  Shown  in  Figures  2.2  and  2.3  are  the  SIMS  profiles  of  the 
unannealed  implants  together  with  the  LSS  Gaussian  and  the  Gaussian  damage 
profiles  constructed  from  the  statistical  values  of  Table  2.1.  Gaussian 
approximations  are  used  for  the  damage  profiles  because  the  detailed  dis¬ 
tributions  are  not  known. 

2. 3  Encapsulation  Procedures 

To  protect  the  surface  of  InP  from  thermal  degradation  during  the 
necessary  high  temperature  anneal  following  implantation,  various  thin  films 
were  studied  as  encapsulants.  A  detailed  study  of  the  protective  qualities 
of  these  caps  is  given  in  Chapter  U  ;  only  the  experimental  details  of  the 
film  depositions  are  discussed  in  this  section. 


11 


Prior  to  thin  film  deposition,  the  samples  were  rinsed  in  ultra¬ 
sonic  baths  of  trichloroethylene,  acetone,  and  isopropyl  alcohol.  The 
unimplanted  samples  used  in  the  encapsulant  studies  of  Chapter  4  were 
additionally  etched  in  0.1%  bromine-methanol  and  rinsed  prior  to  encapsula¬ 
tion. 

2.3.1  Silicon  Dioxide  (SiC>2)  end  Phospho silicate  Glass  (PSG)  Films 

The  chemical  vapor  deposited  (CVD)  Si02  layers  (more  correctly, 

SiO  )  were  deposited  by  oxidation  of  silane  (SiH, ) .  In  this  system  the 
gases  were  pre-mixed  for  10  min  in  a  nitrogen-flushed  cold  wall  reactor. 

The  InP  was  then  heated  rapidly  (<  15  sec)  by  a  graphite  heater  strip  to  a 
temperature  of  400° C  to  initiate  deposition.  The  deposition  rate  for  the 
S102  films  used  in  this  work  was  ~  1500  £  min  \  Both  sides  of  these 
samples  were  encapsulated  with  ~  3500  &  of  S102. 

The  CVD  phosphosilicate  glass  (PSG)  layers  were  deposited  in  the 
same  reactor  as  the  Si02»  by  the  addition  of  phosphine  (PH^)  to  the  reactant 
gases.  The  deposition  temperature  for  the  PSG  layers  was  300°C  and  the 
deposition  rate  was  ~  2500  £  min-^.  Estimates  based  on  AES  data  indicate 

that  these  PSG  layers  contained  4-6  at.  %  P  (6.2  -  9.3  wt.  %  P) .  Both 
sides  of  these  samples  were  encapsulated  with  ~  4000  £  of  PSG. 

The  average  pinhole  density  in  the  S102  films  was  determined  by 

-2 

Cu-decoration  [38]  to  be  130  cm  .  Although  corresponding  measurements  of  the 
PSG  film  pinhole  density  were  not  made,  it  is  expected  that  this  value 
should  be  comparable  to  that  of  the  S102  films  as  the  same  deposition 
system  was  used  for  both  films. 


Figure  2.4:  Schematic  diagram  of  rf  plasma  silicon  nitride  deposition 
system  [39]. 
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2.3.2  Silicon  Nitride  (SiJij)  Films 

The  rf  plasma-deposited  Si^N^  layers  used  in  this  work  were  de¬ 
posited  in  the  system  shown  schematically  in  Figure  2.4.  The  details  of 
this  system  have  already  been  extensively  discussed  in  other  publications 
[39,40].  Measurements  based  on  Rutherford  backseat tering  and  AES  have 
shown  that  Si^N^  films  deposited  in  this  system  contain  oxygen  levels  on 
the  order  of  1  at.  Z  or  less,  which  has  been  shown  to  be  Important  in  GaAs 
encapsulation  [39,41].  These  Si^N^  layers  were  deposited  by  heating  the 
samples  quickly  (<  20  sec)  to  a  temperature  of  350°C  using  a  low  mass  carbon 
strip  and  simultaneously  initiating  the  rf  discharge.  The  deposition  rate 
was  ~  110  £  min  \  All  Si^N^  samples  were  encapsulated  with  ~  2500  X  of 
Sl^N^  on  both  sides. 

Using  Cu-decoration  techniques,  the  average  pinhole  density  in 

_2 

these  Si^N^  films  was  found  to  be  15  cm  . 

2.4  Annealing  Procedures 

Thermal  anneals  were  performed  in  a  Trans-Temp  silica  furnace 
-3  -3  -1 

tube  with  320  cm  volume  and  a  flow  of  200  cm  min  of  desiccated  forming 
gas  (4%  H,  in  N^) .  The  annealing  temperature  was  monitored  by  a  chromel- 
alumel  thermocouple  directly  beneath  the  sample  holder.  To  minimize  thermal 
shock  effects,  samples  were  allowed  to  slowly  (~  2  min)  reach  250“ C  in  a 
cool  furnace  zone  before  insertion  into  the  equilibrated  temperature  zone. 
Timing  of  the  anneals  was  started  when  the  sample  temperature  was  within 
5%  of  the  stated  anneal  temperature.  The  anneals  were  terminated  by  rapidly 
removing  samples  to  the  cold  end  of  the  furnace  and  quenching  them  with  a 
high  velocity  flow  of  forming  gas. 

The  special  procedures  used  for  annealing  samples  in  evacuated, 
silica  ampoules  will  be  discussed  in  Section  4.4.1. 


Auger  Electron  Spectroscopy  (AES)  was  used  in  this  work  to 


analyze  the  encapsulants  discussed  in  Chapter  2  for  traces  of  indium  or 

phosphorus  which  might  have  diffused  through  the  cap  layer  during  the 

annealing  process.  During  AES  a  sample  is  held  in  high  vacuum  and  is 

irradiated  by  a  1-10  keV  electron  beam  of  a  few  yA.  This  excitation 

produces  a  number  of  secondary  processes  of  which  the  most  relevant  is 

the  ionization  of  atoms  within  the  irradiated  sample  volume.  In  such 

atoms,  the  primary  beam  causes  the  excitation  of  inner  shell  (K,L,  and  M) 

-15 

electrons.  Within  a  short  time  (~10  s)  an  outer  shell  electron  decays 
to  fill  this  inner  shell  vacancy.  The  recombination  energy  in  this  decay 
is  either  carried  away  be  a  photon  (X-ray  emission)  or  by  another  outer 
shell  electron  which  is  ejected  from  the  atom  (Auger  electon) .  The  energy 
of  the  Auger  electron  is  determined  by  the  atomic  energy  levels  involved  in 
the  process  and  is  therefore  characteristic  of  the  emitting  atom.  For 
atoms  within  3-20  X  of  the  surface  there  is  a  high  probability  that  the 
Auger  electron  will  not  suffer  an  inelastic  collision,  thereby  maintaining 
its  characteristic  energy.  Since  the  Auger  peaks  are  small  compared 
to  the  broad  background  in  the  electron  energy  spectrum  (arising  from 
inelastically  scattered  electrons  from  deep  atoms)  the  energy  density 
curve  N(E)  is  differentiated,  dN(E)/dE,  to  improve  detection.  Most 
elements  have  several  Auger  transitions  over  a  wide  energy  range,  which 
allows  one  to  identify  the  element  after  consulting  a  table  of  characteristic 
Auger  electron  energies.  Depth  profiling  is  accomplished  by  combining 
Auger  electron  spectroscopy  with  ion  beam  (Ar  ,Xe  )  sputtering  of  the  surface 
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A  schematic  diagram  of  a  typical  AES  system  is  shown  in  Figure 
3.1.  A  cylindrical  mirror  analyzer  is  the  most  commonly  used  dispersing 
element.  Detection  and  amplification  are  provided  by  an  electron  multi¬ 
plier.  A  lock-in  amplifier  or  a  pulse  counter  is  used  to  discriminate 
the  low-level  Auger  peaks  from  the  background  continuum.  To  monitor  the 
intensity  of  several  Auger  transitions  during  depth  profiling  a  multi¬ 
plexer  is  employed.  More  detailed  reviews  of  the  Auger  process,  AES 
elemental  sensitivity,  and  other  AES  system  capabilities  may  be  found  in 
the  literature . [28,  43-45]. 

In  the  course  of  the  encapsulant  studies  Auger  electron  spectra 
were  acquired  with  either  a  Model  545  or  a  Model  595  Physical  Electronics 
Industries,  Inc.  (PHI)  scanning  Auger  microprobe  (SAM).  With  the  PHI  545, 
derivative  mode  spectra  and  conventional  multiplexer  depth  profiles  were 
obtained  using  primary  electron  beam  energies  of  5  keV,  primary  currents 
of  0.4-1. 0  yA,  and  typical  beam  diameters  of  5-10  yra.  The  same  informa¬ 
tion  was  procured  with  the  PHI  595  using  primary  electron  beam  energies 
of  10  keV,  primary  currents  of  0.09  yA,  and  typical  beam  diameters  of 
1  ym.  A  1  keV  Xe+  sputtering  beam  was  used  in  the  PHI  545  and  a  1.3  keV 
Ar+  beam  in  the  PHI  595.  Derivative  mode  spectra  were  taken  from  each 
surface  and  at  several  points  in  the  film  bulk  by  interrupting  the  ion 
sputtering  beam.  To  prevent  the  loss  of  information  from  the  surfaces, 
no  "ion  dusting"  of  the  surface  was  performed.  Wherever  possible  surface 
spectra  were  recorded  from  areas  which  had  been  minimally  exposed  to  the 
primary  electron  beam  so  as  to  avoid  desorption  effects.  Depth  profiles 
were  obtained  by  monitoring  the  Si  KLL,  N  KLL,  0  KLL,  P  LW,  and  In  MNN 
Auger  transitions.  By  using  the  standard  imaging  capabilities  of  the 
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systems,  every  effort  was  made  to  obtain  spectra  and  depth  profiles  from 
areas  free  of  pinholes.  Images  of  an  element’s  distribution  across  a 
sample  surface  were  obtained  by  scanning  the  primary  beam  across  the  surface 
and  intensity-modulating  the  z-axis  of  a  cathode  ray  tube  display  with  the 
signal  from  the  element's  Auger  transition. 

3.2  Secondary  Ion  Mass  Spectrometry  (SIMS) 

A  natural  distinction  can  be  made  between  depth  profiling 
techniques,  based  upon  the  function  of  the  sputtering  process.  Some  techniques 
(such  as  AES)  use  ion  beam  sputtering  to  expose  underlying  material  to 
a  second  probing  beam  which  generates  the  analytical  signal.  Other  depth 
profiling  techniques  use  ion  beam  sputtering  to  generate  sputter-ejected 
material  to  analyzed.  Secondary  ion  mass  spectrometry  (SIMS)  is  an  example 
of  the  latter  technique.  SIMS  was  used  in  this  work  predominantly  to 
examine  the  atomic  redistribution  of  intentionally  present  impurities  such 
as  Cr  or  Fe  and  of  implanted  impurities  in  InP  during  thermal  anneals. 

SIMS  was  also  used  as  a  complementary  method  to  AES  for  characterizing 
the  thin  film  encapsulants. 

In  SIMS  a  sample  is  bombarded  by  a  1-30  keV  ion  beam  (Ar+  ^O^, 

^0  ,  Cs+)  at  low  (nA-yA)  current.  Within  a  depth  of  5-50  X  this  bombard¬ 
ment  produces  sputtered  material  in  a  number  of  forms:  neutral  atoms; 
molecules;  negative  and  positive  secondary  ions;  multiply  charged  secondary 
ions;  and  dimers.  In  standard  depth  profiling  one  is  most  often  concerned 
with  the  first  three  forms.  The  yield  of  secondary  ions  can  vary  from 
0.1-10%  of  the  total  sputtered  material  [45].  Ion  bombardment  of  a  solid 
produces  secondary  ions  by  kinetic  [46]  and  chemical  [47,48]  processes. 


i 


L8 


The  kinetic  process  results  from  the  kinetic  energy  transfer  from  the 
bombarding  ion  to  the  sample  atom.  Sample  atoms  near  the  surface  are  then 
either  ejected  from  the  surface  or  are  excited  to  ionized  and  metastable 
states.  Unbound  electrons  near  the  excitation  volume  neutralize  most 
atoms  before  they  can  escape  the  surface.  Most  secondary  ions  used  for 
mass  analysis  result  from  an  uncharged,  metastable  particle  escaping  from 
the  surface  and  emitting  an  electron  in  vacuum.  The  chemical  process 
relies  on  the  reduction  of  the  unbound  electrons  which  neutralize  the 
ionized  atoms  in  the  bulk  by  the  introduction  or  presence  of  a  chemically 
reactive  species  in  the  analyzed  material.  It  is  assumed  that  both  the 
kinetic  and  chemical  processes  occur  concurrently.  Published  work[47-49] 
has  established  that  secondary  ion  yields  are  enhanced  by  the  presence 
of  an  oxide  layer  on  the  analyzed  surface.  In  order  to  utilize  this 
chemical  effect  oxygen  is  either  introduced  into  the  sample  chamber  or, 
more  commonly,  oxygen  ion  bombardment  is  used  which  results  in  the 
formation  of  a  surface  oxide  by  the  shallow  implantation  of  oxygen  atoms 
into  the  sample.  Enhanced  negative  secondary  ion  yields  have  been  observed 
with  the  use  of  cesium  ion  bombardment  and  are  believed  to  result  from 
a  reduction  in  the  surface  work  function  by  the  chemical  effect  [50,51]. 

More  detailed  reviews  of  the  physical  processes  associated  with  SIMS 
analysis  may  be  found  in  the  literature  [44,45,50,52]. 

A  schematic  diagram  of  the  CAMECA  IMS-3f  used  during  the  SIMS 
studies  is  shown  in  Figure  3.2.  This  model  is  a  direct  imaging  (ion 
microscope)  mass  analyzer  as  originally  proposed  by  Castaing  and  Slodzian[54] . 


By  feeding  the  appropriate  gas  into  a  high  intensity  duoplasmatron  ion 
source,  primary  ions  such  as  Ar+  or  ot  can  be  produced.  Cesium  (Cs"*") 
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Figure  3.2:  Schematic  representation  of  a  CAMECA  IMS-3f  direct  imaging 
secondary  ion  mass  spectrometer (after  CAMECA  [53]). 


ions  are  produced  in  a  separate  ion  source  which  vaporizes  solid  cesium 
and  ionizes  the  gaseous  atoms.  Such  primary  ions  are  then  extracted  from 
the  ion  source  and  accelerated  into  the  keV  energy  range.  Scanning  plates 
then  raster  the  primary  ion  beam  across  the  sample  surface.  Usually  the 
primary  ion  beam  is  scanned  over  a  much  larger  area  than  that  used  for 
analysis,  with  typical  crater  dimensions  varying  from  50-500  ym  on  a  side. 
By  proper  mechanical  or  electronic  aperturing  only  those  ions  originating 
from  the  flat-bottomed  center  of  the  sputtered  crater  give  rise  to  the 
detected  signal.  Ions  from  the  sloping  crater  walls  are  thereby  rejected 
as  these  would  adversely  degrade  depth  resolution [ 50 ] .  The  sample  itself 
is  maintained  at  a  repulsive  polarity  to  the  desired  secondary  ions  for 
analysis.  These  secondary  ions  are  extracted  by  an  electrostatic  im¬ 
mersion  lens  which  maintains  a  point-to-point  location  image  of  the  origin 
of  each  secondary  ion.  A  special  set  of  lenses  referred  to  as  transfer 
optics  transport  the  extracted  secondary  ions  to  a  double-focusing  mass 
spectrometer.  Sputtered  secondary  ions  have  an  energy  distribution  which 
peaks  in  the  range  of  ~  5-10  eV,  but  can  have  low-level  tails  extending 
up  to  several  hundred  eV[52].  The  double- focusing  mass  spectrometer  is 
advantageous  since  it  can  accommodate  both  the  directional  spread  and 
energy  distribution  of  ions,  resulting  in  a  high  transmission  efficiency. 

It  accomplishes  this  by  combining  ion  optical  elements  which  focus  the 
divergent  beams  in  such  a  way  that  the  velocity  dispersions  cancel[55,56] . 
Ion  detection  can  be  recorded  photographically  or  observed  on  a  flourescent 
screen.  Quantization  of  the  ion  signal  is  accomplished  by  either  a 
Faraday  cup  or  electron  multiplier.  The  entire  apparatus  is  interfaced 
to  a  small  desk-top  microcomputer. 
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In  the  studies  discussed  here,  the  0^  primary  beam  (5. 0-7. 7 
keV  Impact  energy)  was  rastered  over  areas  of  ~  250  ym  and  standard  pre¬ 
cautions  were  taken  to  reject  spurious  signals  arising  from  the  crater 
walls.  A  Sloan  Dektak  stylus  instrument  was  used  to  measure  the  resultant 
crater  depths.  Our  experience  and  that  of  others [57,58]  has  shown  the 
uncertainty  in  such  measurements  to  be  +  5-6%.  These  depth  measurements 
together  with  the  assumption  of  a  constant  sputtering  rate  were  used  to 
convert  the  sputtering  time  scale  to  a  depth  scale.  The  SIMS  ion  count 
data  were  calibrated  in  terms  of  concentration  by  integrating  the  area 

under  the  unannealed  profile  (corrected  for  background  effects)  and  set- 
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ting  it  equal  to  the  implanted  dose.  The  P  matrix  ion  signal  was  used 

to  calibrate  the  concentrations  in  the  annealed  sample  to  the  concentration 

scale  of  the  unannealed  implant.  The  concentration  scales  for  ^6Fe  and 
52 

Cr  were  obtained  from  a  standard  prepared  by  implanting  Fe  and  Cr  into 
semi-insulating  InP.  The  conversion  factors  relating  ion  intensity  to 
concentration  were  corrected  for  background  effects. 

3.3  Photoluminescence  (PL) 

Low  temperature  photoluminescence  (PL)  is  a  versatile,  non¬ 
destructive  characterization  technique.  Because  radiative  recombination 
processes  depend  critically  on  the  degree  of  lattice  perfection,  PL  is  an 
extremely  sensitive  technique  for  studying  either  native  or  processing- 
related  defects  at  levels  which  other  characterization  methods  would  fail 
to  detect.  Furthermore,  the  typical  binding  energies  of  shallow  acceptors 
and  the  chemical  shifts  between  them  are  large  enough  in  GaAs  and  InP  that 
PL  can  be  readily  applied  to  the  characterization  of  shallow  acceptors  in 


these  materials  (the  characterization  of  shallow  donors  requires  more 
exacting  techniques  using  far  infrared  spectroscopy).  As  examples  of  such 
work,  PL  has  been  used  extensively  by  members  of  this  laboratory  to 
characterize  implant  damage,  annealing  defects,  and  implanted  impurities 
in  Si[59] ,  ZnSe[60],  GaAsP  [61,62],  GaP[62],  AlGaAs  [62],  and  GaAs[63]. 

In  this  work  PL  has  been  applied  to  various  InP  studies.  It  has  been  used 
here  to  characterize  the  effectiveness  of  the  thin  film  encapsulants.  In 
addition,  it  has  been  used  to  determine  the  annealing  parameters  of  temper¬ 
ature  and  time  required  for  InP  lattice  recovery  from  implant  damage  and 
for  optical  activation  of  Implanted  impurities. 

In  general,  photoluminescence  is  the  recombinative  luminescence 
emitted  from  a  crystal  excited  to  a  nonequilibrium  state  by  irradiation 
from  light  with  an  energy  greater  than  the  crystal's  band  gap.  Because 
the  absorption  coefficient  for  such  light  is  large,  the  excitation  is 
confined  to  a  region  very  near  the  irradiated  surface.  A  large  non¬ 
equilibrium  population  of  electrons  and  holes  is  created  within  this  excited 
region  by  the  absorption  of  energy.  Since  this  population  of  excess 
carriers  is  nonuniformly  distributed  near  the  surface,  the  photo-generated 
carriers  diffuse  away  from  the  excitation  region.  While  diffusing  the 
excess  carriers  are  depleted  by  radiative  and  nonradiative  recombination 
processes  such  that  the  sample  volume  probed  is  confined  to  roughly  a 
diffusion  length  from  the  irradiated  surface.  By  appropriate  choice  of 
excitation  wavelength,  the  absorption  profile  and  generated  carrier 
population  profile  can  be  chosen  to  yield  recombination  luminescence  from 
near-surface  regions  of  particular  interest.  As  the  excess  photo-generated 
electrons  and  holes  diffuse  through  the  sample,  a  variety  of  radiative 
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decay  processes  can  occur,  including  free  exciton,  exciton  bound  to  neutral 
donor  or  acceptor,  donor  to  free  hole,  acceptor  to  free  electron,  donor 
to  acceptor,  or  impurity  to  defect  transitions.  These  transitions  are 
useful  in  implantation  studies  because  by  carefully  monitoring  them  it 
is  possible  to  characterize  the  optical  activation  of  implanted  dopants 
and  also  to  determine  the  extent  of  defects  created  by  post-implantation 
processing.  In  addition,  nonradiative  processes  indirectly  have  a  role  in 
in  implantation/photoluminescence  studies  by  decreasing  the  total  integrated 
sample  luminescence  and  thereby  indicating  the  degree  of  residual  lattice 
damage  following  implantation  and  annealing.  Due  to  self-absorption  effects 
it  is  usually  most  useful  to  collect  the  recombination  luminescence  emitted 
from  the  irradiated,  front  surface  of  the  sample. 

Figure  3.3  is  a  schematic  illustration  of  the  photoluminescence 

system  employed  in  these  studies.  PL  spectra  were  obtained  from  samples 

mounted  in  a  strain-free  manner  in  a  gas-exchange  liquid  helium  cryostat 

(Janis  "Super  Varitemp") .  The  temperature  of  the  sample  mount  during  photo- 

o 

excitation  was  maintained  at  5  K,  as  monitored  by  a  Ge  resistance  thermometer 
embedded  in  the  copper  sample  holder.  Excitation  was  provided  by  the 
5145  X  line  of  a  Coherent  Radiation  Ar-ion  laser  with  an  incident  power 
density  of  ~  10  W/cm  and  a  focused  beam  diameter  of  ~  40  urn.  At  this 
wavelength  the  1/e  point  for  the  penetration  depth  of  the  excitation  beam 
is  939  %  [64].  The  front  surface  luminescence  signal  was  collected  and 
focused  onto  the  entrance  slit  of  a  0.5mSpex  Model  1302  double-grating 
spectrometer.  The  optical  output  was  detected  by  a  cooled  (77°K)  photo¬ 
multiplier  with  S-l  response  using  conventional  lock-in  amplifier 


Figure  3.3:  Schematic  diagram  of  the  photoluminescence  system. 


techniques.  The  spectral  resolution  used  in  this  work  was  6.5  %  (or 
1.1  meV) .  The  (near-bandedge)  spectral  range  examined  was  from  1.44  to 


1.30  eV  (or  8608  to  9535  X).  The  emission  spectra  presented  here  are 
uncorrected  for  system  response. 

Throughout  the  entire  series  of  spectral  examinations  a  single 
control  sample  was  maintained  to  calibrate  for  possible  differences  in 
system  response  between  runs.  All  samples  were  indivudally  examined 
prior  to  processing  to  determine  their  optical  characteristics.  Using 
this  data,  the  integrated  intensity  values  of  the  processed  samples  were 
normalized  to  minimize  differences  in  the  spectral  characteristics  that 
would  result  from  sample  variations.  Each  data  point  in  the  semi-log 
plots  represents  the  average  value  from  at  least  two  samples  that  were 
processed  under  identical  conditions. 

3.4  Electron  Channelling  Patterns  (ECP) 

In  scanning  electron  microscopy,  patterns  very  similar  to  the 
standard  Kikuchi  electron  diffraction  patterns  seen  in  transmission  electron 
microscopy  can  be  obtained  by  rocking  a  collimated  electron  beam  about 
a  point  on  a  crystal  and  monitoring  the  backscattered  electron  intensity. 
Such  patterns  were  first  observed  by  D.G.  Coates [65]  and  are  usually 
referred  to  as  Coates-Kikuchi  lines  or  electron  channelling  patterns(ECP) . 
The  ECP  are  a  function  of  the  electron  beam  orientation  with  respect  to 
the  crystal  lattice  symmetry  planes.  Due  to  this  dependence  on  lattice 
planes,  the  ECP  are  useful  for  a  determination  of  crystalline  perfection 
because  any  change  or  loss  in  lattice  periodicity  will  result  in  an 
altered  or  degraded  channelling  pattern.  The  technique  is  particularly 
useful  for  qualitative  assessments  of  lattice  perfection  in  implantation 


work,  where  the  process  greatly  disorders  and  sometimes  even  amorphizes 
the  implanted  crystals. 

Briefly,  the  physics  of  electron  channelling  patterns  may  be 
explained  as  follows.  When  an  electron  beam  is  incident  on  a  solid  a 
variety  of  energy  loss  events  occur,  some  of  which  result  in  the  back- 
scattering  of  the  incident  electrons.  As  the  angle  of  incidence  from 
sample  normal  increases,  the  fraction  of  backscattered  electrons  also  in¬ 
creases.  Besides  this  usual  backscattering  phenomena  there  occurs  in 
crystalline  solids  an  interaction  between  the  incident  electrons  and  the 
periodic  atomic  lattice  planes.  In  the  diffraction  theory  of  low  energy 
electrons  (below  100  keV)  the  SchrSdinger  equation  for  the  electron 
probability  distribution  in  a  periodic  lattice  yields  a  solution  which  is 
a  superposition  of  two  types  of  Bloch  waves[66,67] .  One  set,  say  type  1 
Bloch  waves,  have  maxima  at  atomic  sites,  implying  strong  atomic  interactions 
with  resultant  high  backscattering  probabilities.  The  other  set,  type  2 
Bloch  waves,  have  maxima  between  atomic  planes,  connoting  weak  atomic  inter¬ 
actions  with  high  transmission,  low  backscattering  probabilities.  The 
changes  in  backscattered  electron  intensity  observed  when  the  incident  beam 
is  rocked  through  various  angles  from  normal  occur  because,  except  at  certain 
angles,  either  type  1  Bloch  waves  (high  backscattering)  or  type  2  Bloch 
waves  (low  backscattering)  are  predominantly  excited.  The  particular  angles 
at  which  type  1  and  type  2  Bloch  waves  are  equally  excited  are  the  Bragg 
angles  9_  satisfying  the  classical  Bragg  diffraction  relation 

D 

nA  *  2d  sin  6g 

For  angles  less  than  9  type  1  Bloch  waves  are  preferentially  stimulated, 


27 


resulting  in  high  intensity  because  of  the  large  backscattering.  For 
angles  greater  than  6  type  2  Bloch  waves  predominate,  resulting  in  lower 

O 

backscattering  intensity.  Therefore,  as  the  incident  electron  beam 
sweeps  through  the  Bragg  angle  for  each  order  of  crystalline  lattice  planes 
there  is  a  discernible  change  in  backscattered  electron  intensity.  It 
is  this  periodic  variation  in  backscattered  intensity  which  defines  the 
electron  channelling  pattern. 

In  the  work  discussed  here,  electron  channeling  patterns  (ECP) 
were  obtained  using  a  JEOL  model  JSM-35  scanning  electron  microscope 
equipped  with  components  (JEOL  35-SDC)  which  permit  selected  area 
channeling  patterns  to  be  observed.  Beam  voltages  of  35  kV  were  used 
to  analyze  areas  with  ~  20  urn  diameters.  The  backscattered  electron 
signals  were  usually  differentiated  to  enhance  the  sharpness  of  the 
normal  channeling  pattern.  To  examine  the  electron  channeling  patterns 
at  various  depths,  some  samples  were  chemically  etched  in  bromine-methanol. 
The  depths  of  these  etched  surfaces  were  measured  with  a  Sloan  Dektak 
stylus. 


4.  ENCAPSULANT  STUDIES 


4.1  Introduction 

If  implantation  is  to  become  a  standard  doping  technique  in  InP, 
it  will  be  necessary  to  develop  reliable  annealing  methods.  Annealing  is 
required  to  activate  the  implanted  impurity  and  also  to  reorder  the 
damaged  crystalline  lattice.  Either  beamed-energy  (laser  or  electron 
beam)  anneals  or  thermal  (furnace)  anneals  can  be  employed,  but  furnace 
annealing  is  the  most  prevalent  method.  Thermal  anneals  of  III-V  compounds 
such  as  GaAs  and  InP  require  care  to  prevent  surface  deterioration  at  the 
common  annealing  temperatures. 

The  thermal  degradation  of  unprotected  InP  surfaces  has  been 
described  in  articles  by  Bayliss  and  Kirk  [68,69]  and  by  Lum  and  Clawson[70]. 
This  surface  damage  results  from  the  preferential  evaporation  of  the  column 
V  component  at  elevated  temperatures.  At  800°C  the  P  vapor  pressure  over 
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InP  is  ~  10  atm.  and  is  generally  orders  of  magnitude  greater  than  that 
of  As  over  GaAs  for  the  temperature  range  considered  here  [71].  Therefore, 

InP  surface  degradation  during  anneals  can  be  much  more  severe  than  in  GaAs. 
Mass  spectroscopic  studies  [72]  indicate  that  the  (001)  InP  surface  evapo¬ 
rates  congruently  up  to  362  +  5°C,  with  higher  temperatures  resulting  in  a 
disproportionate  loss  of  P.  In  addition,  electrical  measurements  of  un¬ 
protected,  n-type  InP  annealed  for  15  min  at  450°C  show  that  a  2  urn 
surface  depletion  zone  is  formed,  presumably  by  P  vacancies  [73].  These 
temperatures  are  comparatively  low  to  the  thermal  anneal  temperatures  of 
700-800°C  which  photoluminescence  [74,75]  and  electrical  measurements [19 , 76, 77] 
indicate  are  needed  to  successfully  reorder  the  damaged  lattice  and  to 
activate  the  implanted  impurity. 


To  prevent  thermal  damage  at  elevated  temperatures  it  is  neces¬ 
sary  to  anneal  InP  in  a  carefully  controlled  environment  [70,78-81]  or  to 
encapsulate  the  surface  with  a  thin  dielectric  layer  having  suitable 
mechanical  and  chemical  properties  to  suppress  evaporation  and  yet  not 
interact  significantly  with  the  material  [74-77,  82-84],  Several  recent 
reviews  of  ion  implantation  III-V  compounds  have  discussed  the  major  role 
encapsulants  can  have  in  determining  the  electrical  characteristics  of 
these  materials  [25-27],  In  this  work  Auger  electron  spectroscopy  (AES), 
secondary  ion  mass  spectrometry  (SIMS)  ,  and  low  temperature  photoluminescence  (PL) 
were  used  to  characterize  the  annealing  properties  of  the  most  commonly 
employed  InP  encapsulants  -  chemical  vapor  deposited  (CVD)  SiC^,  rf  plasma- 
deposited  Si^N^,  and  CVD  phosphosilicate  glass  (PSG) . 

Before  discussing  the  results  of  these  studies,  it  is  necessary 
to  point  out  that  the  anneal  temperature  limits  for  the  SiC^,  PSG,  and 
Si^N^  layers  described  here  were  approximately  650°C,  800°C,  and  840°C 
respectively.  Detailed  analysis  of  the  high  temperature  failure  of  these 
caps  is  presented  in  Chapter  5. 

4. 2  AES  Results 

AES  examination  of  as-deposited  and  annealed  (600°C,  60  min)  Si02 
caps  did  not  indicate  the  presence  of  P  o  In  either  on  the  surface  or  in 
the  film  bulk.  It  has  been  reported  [85]  that  15  min  anneals  of  CVD  SiO^ 
caps  in  InP  similarly  showed  no  AES  P  or  In  within  the  films  for  temperatures 
in  the  range  from  500  -  750°C.  Wilmsen  et_  al_.  [86]  have  reported,  however, 
that  based  on  ESCA  (electron  spectroscopy  for  chemical  analysis)  measure¬ 
ments  In  was  detected  throughout  as-deposited  CVD  Si02  layers  on  InP  when 
the  deposition  temperature  is  340°C  or  greater.  ESCA  analysis  by  Wilmsen' s 


group  of  one  of  our  annealed  (600°C,  30  min)  CVD  Si02  layers  suggested 

that  small  amounts  of  In  and  P  were  distributed  throughout  the  film. 

Because  of  the  large  area  analyzed  by  ESCA  and  because  of  the  pinhole 

density  in  our  CVD  Si02  films,  it  is  unclear  if  these  In  and  P  signals 

represent  actual  outdif fusion  or  a  loss  through  pinholes.  AES  examination 

of  as-deposited  and  annealed  (750°C,  60  min)  PSG  caps  did  not  reveal  the 

presence  of  In  or  P  (above  the  initially  deposited  level)  either  on  the 

film  surface  or  in  the  bulk.  Auger  spectra  from  as-deposited  Si^N^  caps 

also  did  not  show  any  levels  of  P  or  In,  nor  was  P  detected  anywhere  in 

annealed  (750°C,  60  and  30  min)  Si^N^  caps  of  InP.  However,  In  was  found 

on  the  surfaces  of  annealed  (750°C,  60  and  30  min)  Si^N^  caps.  This  In 

was  confined  exclusively  to  the  near  surface  region  and  quickly  dropped 

to  background  levels  within  the  film  bulk.  These  results  agree  with 

AES  measurements  by  Kasahara  jet  _al.  [82]  who  found  In  outdif fusion  and 

no  indication  of  P  outdiffusion  from  Si_N  caps  of  InP  annealed  for  20 

3  + 

min  at  900°C.  It  appears  on  the  basis  of  AES  measurements,  therefore,  that 
In  loss  from  InP  surfaces  does  occur  with  Si^N^  encapsulation. 

Figure  4.1  shows  the  Auger  derivative  mode  spectra  of  the  In  MNN 
transitions  from  several  Si^N^  caps  annealed  for  30  min  at  750°C.  Comparison 
of  Figures  4.1(a)  and  4.1(b)  indicates  that  "thick"  (~  2000  %)  Si^N^  caps 
are  superior  to  thinner  layers  in  their  ability  to  suppress  In  loss  from 
InP  surfaces.  The  captions  "exposed"  and  "unexposed"  in  Figure  4.1  refer 
to  whether  the  Si^N^  cap  was  in  contact  with  the  forming  gas  during  the  anneal 
("exposed")  or  whether,  in  the  same  forming  gas-flushed  furnace,  the  cap 
was  covered  with  a  thin,  demetallized  silica  slide  ("unexposed").  Suppressed 
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Figure  4.1:  Indium  MNN  peaks  in  various  Si^N^  caps  annealed  at  750°C  for 
30  min.  "Exposed"  means  the  cap  was  in  contact  with  forming  gas  during  the 
anneal.  "Unexposed"  means  the  cap  was  covered  by  a  demetallized  silica  slide 
during  the  anneal  [87]. 
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In  loss  has  consitently  been  observed  in  such  covered  Si^N^  films  annealed 
at  temperatures  above  700°C,  as  Figure  4.1(c)  indicates.  Others  have 
reported  using  this  covering  technique  [74],  and  it  has  even  been  used  to 
successfully  anneal  bare  InP  surfaces  up  to  760°C  [81]. 

AES  examination  of  the  encapsulants  on  InP  has  several  limitations. 
These  include:  (1)  overall  atomic  sensitivity  of  only  ~  0.5  at.  %; 

(2)  phosphorus  desorption  by  both  electron  and  ion  beam  irradiation  [88] j 
and  (3)  indium  desorption  stimulated  by  electron  beam  irradiation  [89].  As 
a  result  of  these  sensitivity  limitations,  In  and  P  must  be  present  in 
relatively  large  concentrations  in  the  film  to  be  detected  by  AES,  allowing 
only  gross  characterization  of  the  encapsulants.  For  these  reasons,  SIMS 
was  used  in  addition  to  AES  to  characterize  the  encapsulants. 

4.3  SIMS  Results 

To  supplement  the  information  obtained  with  AES,  the  annealing 
characteristics  of  the  InP  caps  were  examined  using  secondary  ion  mass 
spectrometry  (SIMS) .  The  primary  advantage  of  SIMS  over  AES  is  its  much 
higher  elemental  detection  sensitivity.  Because  of  this,  over  the  last 
several  years  there  has  been  a  growing  acceptance  and  use  of  SIMS  in  the 
analysis  of  dielectric  thin  films  [90-92],  The  technique  is  not  without 
its  problems,  however,  and  before  discussing  the  results  it  would  be 
worthwhile  to  point  out  some  difficulties  in  the  use  of  SIMS  for  dielectric 
thin  film  analysis. 

The  major  problem  which  occurs  during  SIMS  depth  profiling  of  di¬ 
electric  films  is  surface  charging.  This  can  promote  field  enhanced  migra¬ 
tion  of  the  profiled  ions  [93],  More  importantly,  by  affecting  the  energy 
distribution  of  secondary  ion  currents,  such  charging  effectively  degrades  the 
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the  technique's  ultimate  sensitivity  194-96].  To  minimize  this  sensitivity 
degradation,  several  methods  of  reducing  charging  effects  have  been  used 
with  varying  degrees  of  success.  These  methods  include  deposition  of  a 
conducting  film  onto  the  surface  [93,97];  use  of  conducting  grids  in 
contact  with  the  surface  [91,92,94];  and  flooding  of  the  surface  with  an 
electron  beam  [95,96].  None  of  these  methods  has  been  generally  adopted 
as  a  standard  procedure.  Also,  since  surface  charging  is  seldom  totally 
absent  from  SIMS  depth  profiles  of  dielectric  films,  another  problem  lies 
in  finding  a  valid  and  convenient  method  of  reprocessing  data  to  correct 
for  the  sputter  yield  variations  caused  by  charging.  The  prevalent 
technique  is  to  normalize  the  profiles  of  interest  to  the  profile  of  a 
host  element  or  another  element  within  the  film  of  known  constant  concen¬ 
tration  [91,92,96].  To  our  knowledge,  however,  it  has  not  been 
demonstrated  that  the  secondary  ion  yields  of  all  impurities  in  a  given 
matrix  will  respond  uniformly  to  charging. 

For  the  profiles  presented  here  no  special  procedures  were 
performed  to  reduce  surface  charging.  Several  areas  were  analyzed  from 
each  film  and  the  results  were  averaged  to  construct  representative 

impurity  profiles.  After  constructing  such  profiles  the  data  was  re- 
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processed  by  normalizing  all  signals  levels  to  the  Si  profile. 

Examination  of  the  annealed  caps  revealed  that  none  of  these  films 
31 

contained  P  at  levels  above  those  observed  in  as-deposited  films  (where 
31 

the  P  level  may  be  real  or  an  instrumental  background  level) .  It  seems 

30  + 

likely  that  molecular  interference  from  silicon  hydrides  such  as  (  SiH) 

31 

may  account  for  the  difficulties  in  detecting  low  P  levels.  The  hydrogen 
for  these  hydrides  may  originate  within  the  film  itself  or  from  residual 


water  vapor  contamination  In  the  analysis  chamber.  The  In  profiles 
from  as-deposited  and  annealed  (600°C,  60  min)  Si02  caps  also  did  not  differ 
significantly.  The  ^^In  "real  time"  (unnormalized)  profiles  of  PSG  and 
Si^N^  caps  before  and  after  anneals  (750°C,  60  min)  are  shown  in  Figures 
4.2  and  4.3,  along  with  the  normalized  profiles.  It  is  important  to  note 
the  charging  effects  observed  during  these  sputter  profiles  before  con¬ 
sidering  the  ^^In  profiles.  Based  on  AES  measurements  described  in 

section  4.3  it  is  known  that  the  Si  levels  in  these  layers  are  fairly 

28 

uniform.  Therefore,  the  large  variations  in  the  Si  signals  evident  in 

the  unnormalized  profiles  of  Figures  4.2(a)  and  4.3(a)  are  presumably 

the  result  of  charging  effects  during  sputtering  (comparable  charging 
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effects  have  been  reported  for  Si  in  SIMS  profiles  of  SiC^  films)  [91,92]. 

Of  significance  here  are  the  observed  differences  between  unannealed  and 

28 

annealed  films  in  their  response  to  charging.  The  Si  levels  before 
and  after  annealing  in  Figures  4.2(a)  and  4.3(a)  suggest  that  our  as- 
deposited  and  annealed  Si^N^  caps  charge  similarly,  whereas  our  PSG  caps 
showed  considerable  variation  in  their  charging  behavior  depending  on  whether 
or  not  they  were  annealed.  Such  variation  in  charging  behavior  was  also 
observed  between  unannealed  and  annealed  SiC^  films.  These  variations  in 
charging  did  not  affect  our  interpretation  of  Si02  results,  as  identical 
conclusions  were  drawn  from  both  real  time  and  normalized  profiles;  but 
they  do  affect  our  interpretation  of  the  ^^In  profiles  from  PSG  caps. 

The  unnormalized  ^^In  profiles  in  Figure  4.2(a)  suggest  that  there  is 
a  large  accumulation  of  ^^In  near  the  surface  of  the  annealed  PSG  cap. 

The  normalized  ^^In  profiles  of  Figure  4.2(b)  show  a  similar  surface 
accumulation  of  In  but  contain  a  confusing  depressed  ^^In  level  in  the 
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Figure  4.2;  SIMS  profiles  of  ^^In  in  annealed  (750°C,  60  min)  PSG  films  on 
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InP:  (a)  unnormalized  ("real-time")  data  and  (b)  normalized  to  Si  [87]. 
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Figure  4.3:  SIMS  profiles  of  In  in  annealed  (750°C,  60  min)  Si^N^  films 
on  InP:  (a)  unnormalized  ("real-time”)  data  and  (b)  normalized  to  28Si  [87] 


bulk  compared  to  the  unannealed  PSG  film.  It  appears  that  at  levels 
detectable  with  SIMS,  In  outdif fusion  through  PSG  caps  of  InP  is  pos¬ 
sible;  however,  more  work  will  be  needed  to  clearly  resolve  this  issue. 

In  comparison  to  the  PSG  ^^In  profiles,  the  ^^In  profiles  from 
annealed  Si^N^  caps  shown  in  Figure  4.3  are  not  ambiguous.  Both  un¬ 
normalized  and  normalized  profiles  show  an  outdiffusion  of  In  in  annealed 
Si^N^  caps  of  InP,  in  agreement  with  our  Auger  results. 

In  a  more  conventional  application,  SIMS  was  used  to  analyze  the 
InP  surfaces  beneath  the  caps  of  annealed  samples.  This  was  done  to 
determine  the  extent  of  Si  contamination,  if  any,  originating  from  the  cap. 

Si  indiffusion  into  InP  from  the  encapsulants  is  of  significance  since 
Si  is  a  donor  in  InP  [98]  and  will  therefore  influence  the  near  surface 

electrical  characteristics  of  such  contaminated  InP.  To  avoid  sputter 
28 

knock-in  of  Si  from  the  caps  into  the  analyzed  region,  the  caps  were 

removed  from  the  sample  surfaces  prior  to  SIMS  analysis  with  buffered 

28 

hydroflouric  acid.  The  resulting  SIMS  Si  profiles  are  shown  in  Figure  4.4. 

Despite  the  low  annealing  temperature.  Si  indiffusion  into  InP  from  our 

SiC^  cap  is  the  most  pervasive  of  the  three  caps  examined,  with  levels 

above  background  extending  to  ~  1.0  ym  for  60  min  anneals  at  600°C. 

o  28 

Si^N^  caps  annealed  for  60  min  at  750  C  show  higher  than  normal  Si 

levels  extending  into  the  bulk  a  distance  of  ~  0.7  ym.  Of  the  three 

28 

films  analyzed  here,  only  our  PSG  caps  indicated  no  significant  Si 
indiffusion  into  the  underlying  InP  during  annealing.  Besides  ordinary 
diffusion,  it  is  possible  that  the  driving  mechanism  for  the  Si 
contamination  arises  from  strain  at  the  interface  between  the  cap  and 
the  InP  substrate.  This  would  be  consistent  with  the  observation  of  less 
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Si  contamination  from  the  PSG  cap  because  a  PSG  cap  with  5  at.  %  P  is 
better  matched  to  InP  in  terms  of  thermal  expansion  coefficient  than  the 
other  two  caps  [19]. 

4.4  Photoluminescence  Results 

From  the  AES  and  SIMS  results  already  discussed,  no  information 
was  obtained  concerning  possible  P  outdif fusion  through  the  caps.  To 
examine  this  effect,  low  temperature  photoluminescence  measurements  were 
made  of  InP  samples  annealed  in  controlled  atmospheres  and  InP  samples 
annealed  with  the  three  caps.  A  previous  study  of  GaAs  encapsulants  by 
our  group  has  demonstrated  the  usefulness  of  PL  in  characterizing  the 
effectiveness  of  encapsulants  [41]. 

Figure  4.5  shows  the  PL  spectra  from  virgin  samples  of  the  semi- 
insulating,  Fe-doped  InP  (100)  used  in  these  studies.  This  material  exhibits 
three  luminescence  bands  at  1.414,  1.384,  and  1.341  eV.  These  are 
attribured  respectively  to  bandedge  luminescence  (BE),  band  to  acceptor  and 
donor  to  acceptor  (BA/ DA)  luminescence,  and  luminescence  from  the  1-LO 
phonon  replica  of  the  BA/DA  band.  A  fuller  discussion  of  these  features 
and  their  origins  appears  in  Chapter  6. 

4.4.1  Controlled  Atmosphere  Anneals 

In  order  to  interpret  the  PL  spectra  of  InP  annealed  with  our 
caps,  it  was  first  necessary  to  examine  and  characterize  the  spectra  from 
bare  InP  annealed  in  controlled  environments  (vacuum.  In-source,  and  P- 
source).  Such  heat  treatments  should  favor  the  creation  of  annealing 
defects  which  might  subsequently  lead  to  defect  related  luminescence  peaks 
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Figure  4.5:  Photo luminescence  spectra  of  virgin  Fe-doped,  semi-insulating 
InP  (100) .  Numerical  factors  are  amplification  relative  to  that  used  for 


Annealing  ampoules  were  fabricated  from  silica  tubing.  A  small 
dimple  was  formed  in  each  ampoule  to  prevent  contact  between  the  samples 
and  source  materials  during  the  anneals.  The  silica  tubes  were  cleaned 
in  standard  organic  solvents  and  placed  in  aqua  regia  for  24  hrs.  The 
ampoules  were  then  rinsed  in  deionized  water  and  electronic  grade  isopropyl 
alcohol.  Immediately  prior  to  insertion  of  the  sample  and  source  materials, 
the  tubes  were  outgassed  by  heating  them  with  an  oxyhydrogen  flame  while 
being  evacuated  with  a  diffusion  pump.  The  source  materials  were  five-nines 
pure  In  and  six-nines  pure,  red  P  purchased  from  Alfa  Ventron.  After  in¬ 
sertion  of  the  samples  and  source  material  the  ampoules  were  evacuated  to 
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10  Torr  and  sealed  off  to  a  volume  of  about  5  cm  .  The  amount  of  red 

phosphorus  used  in  these  anneals  was  calculated  to  create  a  P  vapor 

pressure  of  0.5  atm.  at  750°C,  sufficient,  according  to  published  vapor 

pressure  curves  of  P  over  InP  [77,99],  to  insure  an  overpressure  of  P. 

For  In-source  anneals,  excess  amounts  of  In  were  inserted  in  the  ampoules 

to  yield  saturation  overpressure  during  the  anneals.  After  annealing,  the 

ampoules  were  quenched  by  inserting  the  tube  end  furthest  from  the  sample 

in  cold  water.  Two-point  probe  electrical  tests  of  sample  surfaces  following 

anneals  at  750°C  in  any  of  the  three  environments  revealed  that  in  all 

cases  the  originally  semi-insulating  InP  had  formed  a  conductive  layer. 

Figure  4.6  exhibits  the  BA/DA  bands  from  samples  annealed  in  the 
three  environments  for  the  temperature  range  from  450-750°C.  Within  the 
limits  of  resolution  used  in  this  work,  no  shifts  were  observed  in  the  peak 
positions  of  the  BE  bands.  Changes  were  noted  in  the  relative  intensities 
of  the  BE  bands  but  these  changes  generally  varied  in  accordance  with  the 
decreases  or  increases  shown  by  the  BA/DA  bands.  Because  of  this  and  for 
simplicity,  only  emission  in  the  region  of  the  BA/DA  band  is  shown  in 
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Figure  4.6:  Photoluminescence  spectra  of  the  BA/DA  band  after  annealing  SI 
InP  at  the  indicated  temperatures  and  times  in  the  following  sealed  ampoule 
environments:  (a)  to  (d) ,  vacuum;  (e)  to  (h) ,  In  vapor;  and  (i)  to  (1),  P 
vapor.  Numerical  factors  are  amplification  relative  to  the  bandedge 
luminescence  of  Figure  4.5  [101]. 


TA 


5 


Figure  4.6.  Three  significant  features  of  these  spectra  are:  (1)  varia¬ 
tions  in  the  BA/DA  band  intensities  with  annealing  temperature;  (2)  a 
new  emission  peak  at  1.378  eV  which  becomes  the  dominant  emission  of  the 
BA/DA  band  for  T  >  650°C;  and  (3)  the  appearance  of  a  new,  small  emission 
peak  at  1.393  eV  in  450°C  and  550°C  vacuum  and  In-source  anneals. 

Comparing  the  amplification  factors  of  the  various  BA/DA  bands 
from  samples  annealed  between  450°C-650°C  to  that  of  the  BA/DA  band  from 
virgin  material,  it  is  clear  that  the  annealed  samples  suffer  a  loss  in 
luminosity.  Presumably  these  intensity  losses  result  from  competition 
between  the  usual  radiative  processes  operative  in  virgin  material  and  new, 
non-radiative  processes  involving  defects  created  by  the  anneals.  This 
luminosity  loss  continues  and  worsens  for  samples  annealed  in  In-source 
ampoules  at  least  up  to  anneal  temperatures  of  750°C  (Figure  4.6(h)).  The 
general  trend  is  reversed,  however,  in  750  C  vacuum  and  P-source  anneals 
(Figures  4.6(d)  and  4.6(£)).  The  BA/DA  band  intensity  in  these  samples  is 
now  greater  than  that  observed  in  the  virgin  spectra  of  Figure  4.5.  The 
increase  of  BA/DA  band  integrated  intensities  above  virgin  levels  is 
correlated  to  the  increasing  dominance  of  the  1.378  eV  transition.  This 
is  certainly  the  case  when  one  compares  the  featureless  BA/DA  bands  of 
Figures  4.6(d)  and  4.6 (A)  to  the  structured  BA/DA  emission  band  of  Figure 
4.6(h).  We  will  discuss  possible  causes  for  these  increases  in  intensities 
when  the  spectra  from  capped  samples  are  shown  in  the  following  section. 

It  is  perhaps  surprising  that  the  vacuum  and  P-source  anneals 


behave  similarly  at  the  high  anneal  temperature  of  750  C.  This  result  is 


understandable,  however,  in  view  of  the  physical  surface  deterioration 
of  the  samples.  As  expected,  the  samples  annealed  in  the  P-source  ampoules 
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showed  no  visible  deterioration  and  retained  their  mirror-like  finish. 

Samples  annealed  in  the  In-source  ampoules  had  the  most  thermally  damaged 
surfaces  of  the  three  sets.  This  may  be  attributable  to  the  In-source 
acting  as  a  large  sink  for  the  phosphorus  which  evaporates  from  the  sample 
surfaces.  Surprisingly,  the  samples  annealed  in  vacuum  exhibited  much 
better  surfaces  than  did  the  corresponding  In-source  annealed  samples.  With¬ 
out  excess  In  present  in  the  ampoule  to  absorb  the  evaporated  phosphorus, 
it  appears  that  the  sample  surface  is  eventually  capable  of  achieving  P 
vapor  phase  equilibrium  within  the  sealed  environment,  preventing  further 
surface  deterioration. 

In  the  650°C  anneals  of  Figure  4.6,  a  weak  shoulder  appears  on 
the  BA/DA  band  at  1.378  eV  which  at  750°C  becomes  the  new,  dominant  emission, 
peak  of  the  band.  This  new  spectral  peak  at  1.378  eV  may  originate  from 
several  sources,  including  Si  contamination  from  the  silica  ampoule  walls 
(Si  is  a  known  contaminant  of  liquid  phase  epitaxial  InP  grown  between 
550°C-715°C)  [100],  annealing-induced  radiative  defects,  or  near-surface 
gettering  of  impurities  from  the  bulk.  At  present  the  specific  origin  of 
this  spectral  feature  at  1.378  eV  cannot  be  identified. 

The  new  feature  of  particular  interest  here  is  the  small  emission 
peak  at  1.393  eV,  which  is  most  pronounced  in  the  450°C  and  550°C  vacuum 
and  In-source  anneals.  In  comparable  P-source  anneals,  this  1.393  eV 
feature  appears  only  as  a  weak  shoulder  on  the  BA/DA  band.  Only  a  few  of 
the  many  published  studies  of  InP  photoluminescence  have  mentioned  emission 
peaks  in  this  energy  region  [102-104].  To  our  knowledge  Williams 
ej:  al.  [102]  were  the  first  to  report  a  feature  at  1.39  eV.  Although  they 
did  not  speculate  on  the  origin  of  this  emission,  they  noted  that  their 
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1.39  eV  line  Increased  In  Intensity  when  the  samples  were  annealed  for 
5  min  in  Hj  at  500°C.  They  also  found  that  the  intensity  of  the  1.39  eV 
line  increased  as  the  epitaxial  layer  from  vapor  phase  epitaxial  InP  was 
etched  off  and  became  strongest  near  the  layer-substrate  interface  (where 
thermal  damage  to  the  substrate  surface  occurs  prior  to  growth  of  the 
epitaxial  layer).  More  recently,  Yamazoe  et  al.  [103],  have  reported 
observing  a  small  emission  line  at  1.39  eV  in  1  hr  anneals  of  InP.  On 
the  basis  of  these  reports  and  from  the  near  suppression  of  the  1.393  eV 
line  in  450°C  and  550°C  P-source  anneals,  we  believe  the  1.393  eV  emission 
line  results  from  P-related  defect  luminescence.  Some  loss  of  P  can  be 
exptected  in  the  P-source  anneals  as  the  atmosphere  is  equilibrated.  To 
minimize  such  transients  of  the  phosphorus  vapor  pressure  in  the  ampoule, 
others  have  reported  converting  red  phosphorus  to  white  phosphorus  by 
heating  the  ampoule  prior  to  annealing  [105]. 

The  annealing  behavior  of  the  1.393  eV  line  is  of  interest,  since 
for  T  >  550°C  the  intensity  of  this  emission  line  decreases  dramatically. 
This  reduction  could  result  from  radiative  recombination  competition  with 
the  newly  emerging  line  at  1.378  eV.  It  is  also  possible,  however,  that 
the  decreased  intensity  may  result  from  a  reduction  in  the  concentration 
of  the  P-related  defects  responsible  for  the  1.393  eV  emission  at  anneal 
temperatures  greater  than  550°C.  This  would  imply  that  the  associated 
defect  is  a  vacancy  complex  rather  than  isolated  P  vacancies,  since  the 
latter's  concentration  would  be  expected  to  increase  with  temperature 
instead  of  decreasing.  Yamazoe  et  al.  [103]  have  reported  similar  de¬ 
creases  in  the  PL  intensities  of  other  InP  annealing  related  defects  for 
T  >  550°C.  In  addition,  their  deep  level  transient  spectroscopy  (DLTS) 
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measurements  of  four  annealing  related  deep  level  traps  indicate  that 
all  four  traps  achieve  their  maximum  concentrations  at  500°C  anneal 
temperatures. 

4.4.2  Encapsulant  Anneals 

The  BA/OA  bands  from  samples  encapsulated  and  annealed  with  our 
three  encapsulants  are  shown  in  Figure  4.7.  All  encapsulant  layers  were 
removed  from  the  annealed  samples  with  hydroflouric  acid  prior  to  PL 
measurements.  Two-point  probe  electrical  tests  of  sample  surfaces  following 
750°C  anneals  with  PSG  or  Si^N^  caps  indicated  that  conductive  surface 
layers  had  formed  during  the  anneal.  Similar  effects  following  anneals 
with  PSG  and  Si^N^  caps  have  also  been  reported  by  others  [19,82]. 

On  the  basis  of  the  appearance  of  the  1.393  eV  line  in  Figures 
4.7(a)  to  4.7(d)  it  is  clear  that  P-vacancy  related  defects  are  created 
in  all  our  SiC^  capped  samples.  It  may  also  be  observed  that  at  temperatures 
between  600-650°C,  the  1.378  eV  emission  is  of  comparable  intensity  to  the 
1.384  line  present  in  virgin  material.  (The  spectra  of  Figure  4.7(d)  was 
obtained  by  analyzing  regions  away  from  the  cracks  which  normally  develop 
in  our  Sit^  caps  at  this  temperature.)  Although  not  evident  in  Figure  4.7, 
for  samples  annealed  at  T  <  650°C,  we  often  observe  the  BA/DA  band  peak 
at  1.381  eV  instead  of  at  the  usual  location  of  1.384  eV.  Local  variations 
in  the  single  crystal  or  in  the  cap  may  account  for  the  inconsistent 
position  of  the  BA/DA  band  peak.  We  believe  that  this  1.381  eV  peak  is 
related  to  the  1.384  eV  line  and  that  the  shift  in  energy  results  from 
the  overlap  of  the  1.384  eV  peak  with  the  emerging  1.378  eV  peak. 
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Figure  4.7:  Photoluminescence  spectra  of  the  BA/DA  band  after  annealing  SI 
InP  at  the  indicated  temperatures  and  times  with  the  following  encapsulants : 
(a)  to  (d),  Si02;  (e)  to  (h) ,  Si^N^;  and  (i)  to  (1),  PSG.  Numerical  factors 
are  amplification  relative  to  the  bandedge  luminescene-  of  Figure  4.5  [101]. 
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By  examing  the  spectra  of  our  PSG  capped  and  annealed  samples 
shown  in  Figures  4.7(i)  to  4.7(a)  it  appears,  because  of  the  presence  of 
shoulders  at  1.393  eV,  that  P-vacancy  defects  are  also  created  to  some 
extent  in  these  samples.  This  could  be  due  to  pinholes  or  to  insufficient 
phosphorus  content  in  our  PSG  to  totally  suppress  phosphorus  outdif fusion. 
Further  work  with  various  types  of  PSG  films  will  be  required  to  clarify 
this  issue.  Comparing  Figure  4.7(0  to  Figures  4.7(d)  and  4.7(h)  it  can 
be  concluded  that  these  PSG  caps  are  much  more  effective  than  SiC^  or 
Si^N^  in  preventing  the  formation  of  the  1.378  eV  emission  feature. 

Judging  by  the  intensity  of  the  1.393  eV  emission,  our  Si^N^  caps 
seem  comparable  to  the  PSG  caps  in  their  ability  to  prevent  substantial 
phosphorus  outdiffusion  (Figures  4.7(e)  to  4.7(h)).  Note  that  for  anneal 
temperatures  T  >  650°C  the  1.378  eV  emission  comes  to  dominate  the  BA/DA 
band  of  our  Si^N^  capped  and  annealed  samples.  Furthermore,  the  luminosity 
of  the  BA/DA  band  in  the  750°C  anneal  (Figure  4.7(h))  is  now  significantly 
greater  than  that  found  in  virgin  samples.  Also,  as  previously  discussed 
with  the  SiC^  spectra,  we  sometimes  find  the  BA/DA  band  peak  in  Si^N^ 
capped  samples  annealed  at  T  <  650°C  at  1.381  eV  rather  than  at  1.384  eV. 

It  was  mentioned  in  the  previous  section  that  the  emission  at 
1.378  eV  may  originate  with  Si  indiffusion,  annealing-induced  radiative 
defects,  or  near  surface  gettering  of  impurities  from  the  bulk.  A 
comparable  downward  energy  shift  of  6  meV  by  the  BA/DA  band  peak  has  also 
been  recently  reported  by  Shanabrook  et  aJL.  [106]  for  Si^N^  capped  and 
annealed  SI  InP(Fe) .  From  temperature  dependent  PL  measurements  they 
attribute  the  new  emission  to  donor-acceptor  pair  recombination,  with 
the  implication  that  the  donor  concentration  has  increased.  They  conclude 


that  such  donors  could  be  provided  by  Si  indiffusion  from  the  Si^N^  cap 

or  from  accumulation  and  redistribution  of  donors  from  the  bulk  to  the 

surface.  From  the  appearance  of  the  spectra  in  Figures  4.7(h)  and  4.7(£) 

28 

and  from  the  SIMS  Si  profiles  shown  in  Figures  4.7(b)  and  4.7(c)  of  the 

preceding  article  it  certainly  appears  likely  that  Si  could  play  a  role 

in  the  1.378  eV  peak.  It  is  not  clear  at  this  point,  however,  why  the 

1.378  eV  peak  does  not  similarly  dominate  the  BA/DA  band  of  Sit^  capped 

28 

and  annealed  samples  when  the  SIMS  Si  profiles  from  these  samples 
showed  the  deepest  indiffusion. 

In  our  high  temperature  ampoule  anneals  (Figures  4.6(d)  and 
4.6(£))  and  again  in  our  750°C,  Si^N^  anneal  (Figure  4.7(h))  we  have 
observed  BA/DA  band  and  BE  band  intensities  greater  than  those  found  in  virgin 
material.  Such  increases  in  bandedge  luminosity  have  been  reported  in  the 
past  for  annealed,  semi-insulating  GaAs(Cr)  [107],  Recent  reports  have 
also  described  increases  in  the  bandedge  emission  of  InP  samples  capped 
with  Si^N^  and  annealed  [74,106].  It  has  been  noted  here  and  elsewhere  [101] 
that  these  increases  in  luminosity  are  correlated  to  the  emergence  of  the 
1.378  eV  feature.  Such  increases  in  luminosity  may  result  from  the  an¬ 
nealing  out  of  native  defects  in  bulk  InP,  surface  accumulation  of  im¬ 
purities  with  greater  radiative  efficiency  than  those  normally  present, 
or  from  a  reduction  of  the  surface  space  charge  layer  due  to  the  presence 
of  the  impurity  and/or  defect  involved  in  the  1.378  eV  emission. 

Figures  4.8(a)  and  4.8(b)  plot  the  integrated  intensities  of  the 
BE  and  BA/DA  bands  for  Si^N^  and  PSG  capped  and  annealed  samples.  Before 
pointing  out  the  significance  of  these  curves  we  must  comment  on  a  feature 
not  discussed  in  this  article  which  affects  our  Si0N.  data.  We  have 
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Figure  A. 8:  Integrated  intensities  of  (a)  the  BE  band  and  (b)  the  BA/DA  band 
from  SI  InP  capped  with  either  Si^N^  or  PSG  and  annealed  for  60  min  at  the 
indicated  temperatures  [101]. 


observed  that  the  spectra  from  samples  annealed  with  Si^N^  at  T  >  650  C 
(and  also  from  some  high  temperature  P-source  anneal  spectra)  contain  a 
discernible  tail  of  luminescence  on  the  low  energy  side  of  the  BA/DA  band 
(Figures  4.6(Jl)»  4.7(g),  and  4.7(h)).  This  broad  tail  extends  from  ap¬ 
proximately  the  low  energy  side  of  the  BA/DA  band  to  the  position  of 
the  1-LO  phonon  band,  with  weak  structures  sometimes  appearing  at 
1.363  eV  and  1.355  eV.  In  this  spectral  range,  Heim  [108]  has  reported 
a  band  centered  at  1.358  eV  which  is  associated  with  long  radiative  decay 
times,  and  he  attributes  it  to  a  DA  transition  involving  a  slightly 
deeper  acceptor  than  the  one  responsible  for  the  peak  of  the  BA/DA  band. 
White  et^  al.  [109]  have  observed  a  transition  in  this  region  at  1.367  eV, 
which  they  identify  as  a  DA  transition  to  Cd  acceptors.  It  is  feasible  that 
the  luminescence  we  observe  in  this  region  originates  with  an  acceptor 
gettered  to  the  surface  from  the  bulk,  but  at  present  our  data  does  not 
support  specific  indentification.  Whatever  the  origin,  the  effect  of 
this  luminescence  tail  has  not  been  included  in  the  measurements  of  the 
Si^N^  BA/DA  integrated  band  intensities  plotted  in  Figure  4.8.  Considering 
the  Si^N^  spectra  of  Figure  4.7  together  with  the  results  of  Figure  4.8 
it  becomes  clear  that  the  BE  and  BA/DA  intensities  in  the  Si^N^  capped  and 
annealed  samples  exceed  virgin  levels  when  the  1.378  eV  emission  dominates 
the  BA/DA  band.  More  importantly,  it  can  be  concluded  from  Figure  4.8  that 
our  PSG  caps  maintain  the  PL  properties  of  the  underlying  InP  samples 
with  much  higher  integrity  than  do  our  Si^N^  caps. 
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4.5  Summary 

The  results  of  these  studies  of  InP  annealing  encapsulants  can 

be  summarized  as  follows: 

(1)  Both  AES  and  SIMS  indicate  that  In  diffuses  through  our  Si^N^  caps 
when  annealed  at  750°C  for  60  min.  Few  conclusions  could  be  drawn 
concerning  outdiffusion  in  either  SiO^  or  PSG  caps.  Some  SIMS  results 
suggest  that  In  outdiffuses  through  PSG  caps  for  60  min,  750°C 
anneals,  but  more  work  will  be  required  to  substantiate  this. 

(2)  It  is  found  from  AES  that  thick  (~  2000  X)  Si^N^  caps  are  superior  to 
thinner  layers  in  their  ability  to  suppress  In  loss  from  InP  surfaces. 
Additionally,  AES  measurements  consistently  indicate  that  In  loss  from 
Si^N^  encapsulated  samples  can  be  reduced  by  covering  the  Si^N^ 
layer  with  a  thin  silica  slide  during  the  anneal. 

(3)  On  the  basis  of  SIMS  results.  Si  indiffusion  is  found  to  be  significant 
in  our  SiO^  and  Si^N^  capped  and  annealed  samples.  The  PSG  caps 
examined  here  showed  little  Si  contamination  of  the  underlying  InP 
after  annealing. 

(4)  PL  spectra  from  samples  annealed  in  controlled  environments  (vacuum. 
In-source,  and  P-source)  indicate  that  a  spectral  line  at  1.393  eV 
results  from  P  vacancy-related  defect  luminescence. 

(5)  The  presence  of  the  1.393  eV  line  suggest  that  at  least  up  to  anneal 
temperatures  of  650°C,  our  Si02  caps  allow  P  loss  from  the  sample 
surface.  A  weak  feature  at  1.393  eV  is  also  evident  in  the  spectra 
of  PSG  capped  and  annealed  samples,  suggesting  some  P  loss.  This  may 
occur  due  to  pinholes  or  to  inadequate  phosphorus  in  these  PSG  films. 
Our  Si^N^  caps  seem  comparable  to  the  PSG  layers  in  their  ability  to 
prevent  substantial  P  outdiffusion. 
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(6)  High  temperature  anneals  (T  ~  750°C)  with  Si^N^  result  in  the  emergence 
of  a  new  peak  at  1.378  eV  which  dominates  the  BA/DA  emission.  This 
emission  at  1.378  eV  may  originate  with  Si  indiffusion,  annealing- 
induced  radiative  defects,  or  near  surface  getter ing  of  impurities 
from  the  bulk.  Associated  with  this  are  BA/DA  band  and  BE  band  inte¬ 
grated  intensities  greater  than  those  found  in  virgin  material. 

(7)  In  comparison  to  SiC^  or  Si^N^,  PL  measurements  indicate  that  PSG 
caps  best  preserve  the  spectral  features  and  the  luminosity  levels 
of  InP  after  annealing. 
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5.  HIGH  TEMPERATURE  LIMITS  ON  ENCAPSULATION 

5.1  Introduction 

Perhaps  the  most  prevalent  problem  with  thin  film  encapsulants 

is  their  tendency  to  fail  at  high  anneal  temperatures  by  cracking,  peeling, 

or  blistering  [19,110].  Such  film  failure  is  usually  attributable  to  thermal 

expansion  mismatch  with  the  underlying  material,  but  other  factors  such 

as  chemical  interactions  can  play  a  role  [110].  In  addition,  there  is 

considerable  variation  in  the  properties  of  thin  films  deposited  by  the 

same  method  at  various  laboratories.  The  anneal  temperature  limits  for  the 

o 

SiO^,  PSG,  and  Si^N^  layers  discussed  here  were  approximately  650  C, 

800°C,  and  840°C  respectively.  Our  CVD  SiO^  films  fail  by  cracking  at 
temperatures  in  excess  of  650°C  whereas  CVD  SiO^  films  on  InP  prepared  at 
other  labs  have  been  reported  to  survive  temperatures  up  to  700°C  before 
cracking.  Our  CVD  PSG  films  also  fail  by  cracking  in  a  fashion  similar 
to  our  SiO^  layers.  Contrary  to  other  reports  [84]  cracking  of  Si^N^ 
caps  on  InP  is  not  a  universal  phenomena  for  T  >  750°C.  We  routinely 
observe  Si^N^  layers  surviving  60  min  anneals  at  800°C  without  cracking. 

The  causes  of  the  high  temperature  limits  on  the  Si02  and  Si^N^  caps  are 
discussed  in  this  chapter. 

5.2  Si02  Cracking 


5.2.1  Experimental  Observations 

The  major  difference  observed  in  this  work  between  Si02  encap¬ 
sulated  InP  annealed  below  650°C  and  that  annealed  above  650°C  is  the 
cracking  phenomena  shown  in  Figure  5.1.  This  cracking  behavior  was  seen 
in  all  our  samples  annealed  above  650°C.  Other  results  indicate  that  the 
temperature  for  cracking  of  Si02  films  is  not  unique  but  can  vary  in  the 
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range  from  650  to  750°C  [19,111].  To  determine  the  role  of  film  thickness, 
samples  were  prepared  with  Si02  layers  of  ~  1000  X,  1500  X,  2300  X,  and 
3000  X  on  both  sides  and  annealed  at  750°C  for  30  minutes.  Each  of  these 
samples  showed  cracking  of  the  Si02  encapsulant  similar  to  that  shown  in 
Figure  5.1. 

It  was  possible  to  observe  the  evolution  of  these  cracks  through 
the  semi-transparent  liner  of  the  furnace.  Following  insertion  into  the 
equilibrated  temperature  zone,  small  cracks  appear  at  the  edges  of  the 
samples  within  5  minutes.  Generally,  samples  are  cleaved  into  a  rectangu¬ 
lar  shape.  The  cracking  begins  at  the  long  edges  of  such  samples  and 
propagates  essentially  parallel  to  the  short  edges.  Within  this  time 
period  no  cracks  originate  in  the  central  portions  of  the  samples  unless 
a  gross  pinhole  develops  there.  As  time  progresses  these  cracks  increase 
in  number  and  propagate  from  the  edges  toward  the  center.  Cracks  running 
parallel  to  the  long  axes  of  the  samples  generally  develop  after  the  major 
cracks  meet  near  the  center  of  the  sample. 

Figure  5.1(b)  shows  three  distinct  regions  in  the  neighborhood  of 
a  crack  network.  The  region  which  comprises  the  majority  of  this  photograph 
is  the  uncracked  SiC^  encapsulant.  We  shall  label  the  dark  areas  within 
the  cracks  region  I  and  the  light,  metallic-appearing  areas  in  these  cracks 
region  II.  In  addition,  close  examination  of  the  uncracked  Sit^  in  Figure 
5.1(b)  reveals  the  presence  of  small,  hairline  features  between  the  major 
cracks.  It  is  assumed  that  with  further  heat  treatment  these  would 


develop  into  major  cracks  also. 
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A  JEOL  scanning  electron  microscope  (SEM)  was  used  to  examine 
the  cracks  in  greater  detail.  Figure  5.2  is  a  SEM  micrograph  showing 
the  three  regions  mentioned  above.  The  central  portion  of  this  figure 
corresponds  to  region  II  and  the  darkened  edges  to  region  I. 

Two  other  features  can  also  be  seen  with  improved  clarity  in 
Figure  5.2.  One  of  these  features,  seen  as  a  line  in  Figure  5.1(b),  is 
a  clearly  visible  track  running  through  the  large  crack  structure.  The 
width  of  this  small  track  is  ~  1.5  pm,  while  that  of  the  larger  crack  is 
~  25  pm.  The  path  of  this  1.5  pm  feature  does  not  appear  to  be  associated 
with  any  major  crystalline  axes,  in  contrast  to  the  25  pm  feature  which 
is  oriented  piecewise  in  (100)  directions  as  seen  more  clearly  in  Figures 
5.1(b)  and  5.3(a).  It  is  also  clear  from  Figures  5.1(b)  and  5.2  that  the 
25  pm  feature  is  not  symmetrically  centered  about  the  1.5  pm  track. 

The  other  feature  revealed  in  Figure  5.2  is  the  apparent  shearing 
of  the  SiO^  layer.  What  appears  to  be  a  tearing  of  the  Si0„  film  can  be 
seen  in  the  vicinity  of  the  1.5  pm  track.  Examination  of  the  1.5  pm 
feature  in  areas  corresponding  to  region  I  show  this  same  shearing  of  the 
Si02  layer. 

A  series  of  photographs  were  taken  after  the  Si02  layer  was 
removed  with  buffered  HF.  Two  photographs  from  this  series  are  shown  in 
Figure  5.3.  Figure  5.3(a)  is  a  photograph  taken  with  the  microscope  focused 
on  the  surface  plane  of  the  sample.  The  features  visible  in  this  photo¬ 
graph  correspond  to  region  II  in  Figure  5.1(b).  The  speckled  surface  of 
these  regions  closely  resembles  the  surface  of  bare  InP  after  annealing 
for  30  minutes  at  650°C  in  an  In  overpressure.  This  suggests  that  these 
regions  are  P-depleted  InP  islands.  The  dark  regions  in  Figure  5.3(a) 


Figure  5.2:  SEM  micrograph  of  a  cracked  region  in  Si02  encapsulated  InP 
annealed  at  750  C  for  30  min.  The  central  region  is  a  P-depleted,  In-rich 
material  buildup  in  a  thermal  etch  trough  [110]. 
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Figure  5.5:  Auger  spectra  from  the  sample  analyzed  in  Figure  5.4  but  with 
the  SiC^  layer  removed.  [110]. 
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correspond  to  the  areas  labeled  region  I  in  Figure  5.1(b).  With  the  Sit^ 
layer  now  removed,  we  can  observe  that  these  areas  are  actually  troughs 
in  the  sample  surface.  Focusing  the  microscope  into  these  troughs  yields 
the  image  of  Figure  5.3(b).  It  can  be  seen  that  the  bottoms  of  these 
troughs  (corresponding  to  region  I  )  are  populated  by  In  droplets.  Such 
droplets  arise  when  phosphorus  evaporates  from  the  InP  substrate.  In  this 
situation  In  is  left  on  the  surface  and  this  In  melt  continues  to  etch  the 
InP  substrate  in  order  to  maintain  equilibrium  saturation  with  phosphorus 
at  the  substrate  temperature  [70,99].  The  rectangular  shapes  of  these 
droplets  is  characteristic  of  thermal  etching  of  (100)  InP. 

Figure  5.4  exhibits  Auger  spectra  taken  from  a  sample  annealed 
for  30  minutes  at  750°C.  The  ~  2300  X  SiC^  layer  was  left  on  this  sample 
to  permit  analysis  of  the  three  previously  mentioned  areas  using  a 
stationary  electron  beam.  Spectra  taken  from  the  uncracked  regions  of 
the  SiO^  revealed  the  presence  of  In  peaks  at  298,  344,  404,  and  410  eV. 

We  also  observed  P  peaks  at  120  and  1859  eV.  These  peaks  do  not  appear 
in  the  spectra  of  either  as-deposited  Si02  or  Sit^  annealed  at  600°C  for 
30  minutes,  unless  the  analyzed  region  includes  pinholes  [87].  It  is 
unclear  whether  these  signals  originate  from  In  and  P  which  has  diffused 
through  the  film,  or  from  material  deposited  on  the  Si02  surface  from 
nearby  cracks.  Examination  of  the  Auger  spectra  from  region  II  (mottled, 
metallic  colored  areas  in  figure  5.3(a))  reveals  several  differences 
from  the  spectra  of  the  uncracked  SiC^*  The  Si  and  0  signals  from  region 
II  are  the  weakest  obtained  from  the  three  regions.  Very  little  Si  is 
actually  detected  within  this  region,  as  can  be  seen  in  Figure  5.4.  The 
weakened  0  signal  may  arise  from  an  oxide  on  the  In-like  surface  of  region 


500  1000  1500  2000 

Electron  Energy  (eV) 


Figure  5.3:  Auger  spectra  from  InP  annealed  at  750°C  for  30  min  with  ~  2500 
X  of  CVD  Si02«  This  series  was  taken  with  the  cracked  Si02  still  on  the  InP 
surface  [110] . 
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Figure  5.4:  InP  annealed  at  750°C  for  30  min  with  a  ~  2500  X  Si02  layer.  I 
layer  was  then  removed  in  buffered  HF.  (a)  Microscope  focused  on  the  sample 
surface  plane.  (b)  Microscope  focused  into  the  bottom  of  the  thermal  etch 
troughs.  [HOT 
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II.  Additionally,  the  ratio  of  the  P  to  In  signals  is  much  smaller  than 
that  found  in  either  the  spectra  of  the  uncracked  Si02  of  Figure  5.4  or 
the  spectra  of  the  bare  InP  surface  of  Figure  5.5.  Auger  spectra  from 
region  I  (the  region  populated  by  In  droplets  seen  in  Figure  5.3(b)) 
show  In  and  P  signals  essentially  comparable  to  region  II  spectra.  Region 
I  can  be  clearly  distinguished  from  region  II,  however,  on  the  basis  of 
the  Si  and  0  Auger  signals.  Whereas  region  II  exhibited  no  Auger  Si 
signal,  region  I  has  Si  signals  at  76  -  92  eV  and  at  1620  eV  comparable 
in  intensity  to  the  Si  signals  from  the  uncracked  Si02.  Although  the  Si 
signals  of  region  I  and  the  Si02  layer  are  comparable,  the  Auger  0  signal 
from  region  I  is  consistently  weaker  than  that  obtained  from  uncracked 
Si02. 

Following  this  analysis  we  removed  the  Si02  layer  from  the 
cracked  sample  using  buffered  HF.  Auger  spectra  were  taken  again  in  the 
three  regions  and  are  shown  in  Figure  5.5.  Essentially  no  Si  is  now 
detected  in  regions  I  and  II,  indicating  that  these  regions  were  covered 
by  a  thin  SiO^  layer  in  the  previous  Auger  analysis.  This  is  consistent 
with  the  appearance  of  these  regions  in  Figure  5.1(b),  5.2,  and  5.3.  Also, 
the  differences  in  spectra  between  the  uncracked  InP  surface  and  region  II 
indicate  that  the  latter  region  is  a  P-depleted/In-rich  region. 

5.2.2  Film  Stress 

Cracks  will  occur  in  deposited  thin  films  from  stress  which 
exceeds  the  film  strength  or  adhesive  interface  strength.  Such  stress  has 
two  components:  (1)  an  intrinsic  stress,  which  is  independent  of  the 
substrate  but  dependent  on  system  parameters  such  as  deposition  technique 
and  deposition  rate;  (2)  thermal  stress  which  arises  from  the  mismatch 
in  thermal  expansion  coefficient  between  the  film  and  substrate  at 


TABLE  5.1 

PROPERTIES  OF  CVD  Si02  FILMS  AND 
EXPANSION  COEFFICIENTS  OF  InP  AND  GaAs 


(Young's  Mo dulus)/(l-Pois son's  Ratio) 

Si02  Thermal  Expansion  Coefficient 

GaAs  Thermal  Expansion  Coefficient 

InP  Thermal  Expansion  Coefficient 

Calculated  GaAs  Thermal  Stress 
(Td  -  400°C,  Ta  -  750°C) 

Calculated  InP  Thermal  Stress 
(Td  -  400° C,  Ta  -  750°C) 


E/(l-v) 

12  -2 
1x10  dyne  cm 

ctf(Si02) 

5. 5xl0-7  °C'1 

ag(GaAs) 

6.6xl0_6  °C-1 

os(InP) 

5.0xl0-6  °C_1 

„  Q 

|ot(GaAs) | 

2 . 1x10  dyne  cm 

9 

|crt(InP)| 

1.6x10  dyne  cm 
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temperatures  differing  from  the  film  deposition  temperature.  The  total 
stress  in  the  film  is  given  by: 


where  a  is  the  intrinsic  stress,  E  is  the  Young's  modulus  of  the  film,  v 
is  the  film's  Poisson  ratio,  and  ag  are  the  thermal  expansion  coeffi¬ 
cients  of  the  film  and  substrate  respectively  [112].  As  is  commonly  done, 
we  will  use  constant  values  for  the  a's  in  our  calculations.  Values  are 
reported  in  the  literature  for  E  and  v,  although  the  ratio  E/(l-v)  is 
more  commonly  extracted  from  data.  For  purposes  of  comparison  we  have 
calculated  the  thermal  stress  in  SiO^  encapsulated  InP  and  GaAs  annealed 
at  750°C  using  average  values  from  the  literature  as  shown  in  Table  5.1  [112-120]. 
For  films  deposited  at  400°  C  and  annealed  at  750°C  the  thermal  stress  in 
the  Si02/GaAs  system  is  ~  30%  greater  than  that  in  the  Si02/lnP  system. 

To  check  the  quality  of  our  CVD  Si02  we  encapsulated  Cr-doped 
GaAs  along  with  Fe-doped  InP  and  annealed  the  samples  at  750°C  for  30  min. 
Semi-insulating  bulk  crystals  were  chosen  to  minimize  differences  in 
crystal  quality.  The  back  surfaces  of  these  samples  were  lapped  to  yield 
a  uniform  thickness  in  all  samples,  and  the  samples  were  cleaved  to 
identical  dimensions  to  minimize  differences  in  stress  that  might  result 
from  sample  size.  Depositions  were  performed  simultaneously  on  the  GaAs 
and  InP  to  insure  identical  intrinsic  stress.  In  addition  to  the 
deposition  system  described  in  Section  2.3.1,  a  second  CVD  system  with  a 
deposition  temperature  of  500  C  and  a  deposition  rate  of  ~  200  X/min  was 
employed.  Following  anneals  at  750°C  for  30  min,  the  Si02  layers  on  InP 


Figure  5.6:  InP  (a)  and  GaAs  (b)  samples  encapsulated  with  SK^  and  annealed 
at  750  C  for  30  min.  These  samples  were  encapsulated  at  500°C  with  a  Si02 
deposition  rate  of  ~  200  X/min.  An  identical  set  of  samples  encapsulated  in 
another  CVD  system  at  400°C  and  a  Si02  deposition  rate  of  ~  1500  X/min 
yielded  similar  results  [110]. 


were  cracked  in  both  sets  while  the  Si02  layers  on  GaAs  were  intact  in  both 
groups.  Figure  5.6  compares  for  one  such  set  of  samples,  the  results  of  an¬ 
nealing  at  750  C  for  30  min.  These  results  together  with  our  calculation 
of  stress  in  these  two  systems  indicate  that  the  cracking  observed  in  InP 
must  be  initiated  by  a  mechanism  other  than  thermal  stress. 

5.2.3.  Discussion 

We  believe  that  the  cracking  phenomenon  observed  in  the  annealing 
of  encapsulated  InP  is  the  result  of  an  interaction  between  mechanical 
stress,  chemical  reactions,  and  possibly  defects  in  the  film  itself  (e.g., 
pinholes).  In  this  section  we  will  discuss  the  data  in  these  terms. 

The  time  evolution  of  the  cracking  as  described  in  Section  5.2.1 
suggests  that  there  is  a  stress  mechanism  at  work.  This  stress  is  presum¬ 
ably  responsible  for  the  1.5  ym  feature  and  the  associated  shearing  of  the 
Si02  seen  in  Figure  5.2.  The  fact  that  the  1.5  ym  feature  does  not  follow 
major  crystalline  directions,  and  the  smooth  but  random  transitions  it 
makes  across  the  surface  indicate  it  is  governed  by  forces  other  than 
those  responsible  for  the  25  ym  features.  We  investigated  whether  the 
1.5  ym  fracture  could  be  caused  by  residual  mechanical  damage  of  the  sur¬ 
face  from  polishing;  however,  we  found  that  samples  with  newly  exposed 
surfaces  obtained  with  a  0.1%  bromine-methanol  etchant  exhibited  the  same 
Si02  cracking  as  the  unetched  samples.  It  is  also  unlikely  that  the  film 
shearing  results  from  a  P  pressure  buildup  at  the  interface.  This  would 
more  likely  lead  to  film  lift-off  or  peeling  by  adversely  affecting  film 
adhesion.  Experience  with  Si^N^  encapsulated  InP  has  shown  that  when 
phosphorus  evaporates  from  an  opening  in  the  film  it  does  not  cause  film 
cracking,  but  rather  a  rectangular  etch  pit  which  grows  with  time.  There¬ 
fore,  we  conclude  that  the  1.5  ym  feature  represents  the  shearing  of  the 
Si0o  film  as  a  result  of  stress. 
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According  to  our  stress  calculations,  thermally  induced  stress 
alone  cannot  explain  the  cracking  observed.  It  appears  most  likely  that  the 
cracking  is  initiated  by  a  local  reaction  of  the  SK^  film  at  the  InP 
surface.  Recently,  Wager  and  Wilmsen  have  cataloged  possible  reactions 
between  SK^  and  InP  [121].  They  conclude  that  within  the  temperature  range 
25°  -  350°C  no  reaction  is  expected  between  SiC^  and  InP  or  its  oxides. 
However,  reactions  of  Si  and  SiO  with  the  oxides  of  In  and  P  are  possible, 
as  is  a  weak  reaction  of  Si  with  InP.  The  lack  of  adequate  entropy  data 
prevented  them  from  taking  into  account  the  temperature  dependent  entropy 
terra  in  their  calculations  of  the  Gibbs  free  energy  for  the  reactions. 

This  term  might  be  expected  to  become  important  at  the  temperatures 
considered  here,  possibly  making  a  SiC^-InP  reaction  feasible.  Also  the 
stoichiometry  of  our  layers  may  not  strictly  be  SiC^,  resulting  in  reactions 
between  Si  or  SiO  with  the  native  oxide  on  the  surface  [122,123].  Recent 
profiles  taken  with  ESCA  indicate  that  there  is  a  thick  I^O^/^O^  layer 
at  the  interface  of  samples  with  CVD-S102  deposited  at  450°C  [86,124].  It  is 
presently  uncertain  if  there  also  exists  an  interfacial  layer  of  phospho- 
silicate  glass  [124].  The  presence  of  such  interface  layers  would  invalidate 
our  stress  calculations,  since  a  multi-layer  system  would  have  to  be 
considered.  Thermal  expansion  mismatch  in  such  a  multi-layered  system 
could  give  rise  to  the  stress  needed  to  initiate  the  cracking.  Alterna¬ 
tively,  a  reaction  of  Si  or  0  with  the  substrate  could  lead  to  a  mechani¬ 
cally  weakened  film,  making  cracking  possible.  Our  Auger  data  indicates 
the  possible  depletion  of  0  in  the  region  I  areas.  Figure  5.7  also  illus- 
strates  this  possibility.  An  oxygen-depleted  fault  to  the  left  of  the 
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major  crack  in  Figure  5.7(b)  is  not  visible  in  the  absorbed  current  image  of 
5.7(a),  and  is  presumably  a  crack  in  its  formative  stage.  The  question  con¬ 
cerning  possible  Si  reactions  is  presently  unresolved. 

Film  defects  such  as  pinholes  may  also  contribute  to  the  cracking 
phenomena  observed,  although  they  do  not  seem  responsible  for  the  initiation. 
It  would  be  difficult  to  explain  the  time  evolution  behavior  in  terms  of  ran¬ 
domly  distributed  pinholes.  It  is  observed,  however,  that  the  distribution 
of  cracks  across  the  surface  is  nonuniform.  This  may  be  correlated  to  the 
nonuniform  distribution  of  pinholes  we  have  observed  with  Cu-decoration  of 
these  films. 

5.3  High  Temperature  Anneals  with  SijN^  Caps 

As  mentioned  earlier,  our  Si^  caps  routinely  survive  60  min 
annea] j  ah  800°C  without  cracking;  but  another  problem  limits  us  to  anneal 
temperatures  T  <  840° C.  This  temperature  limit  is  the  result  of  the  volatile 
thermal  etching  characteristics  of  InP  and  may  therefore  similarly  affect 
other  caps  annealed  at  high  temperatures.  Figure  5.8(a)  is  a  schematic 
representation  of  a  Si^N^  encapsulated  InP  sample  prior  to  annealing.  Micro¬ 
scopic  examination  of  the  cleaved  edges  of  InP  samples  reveal  that  these  edges 
are  not  smooth  but  consist  of  families  of  parallel  faces  along  the  cleavage 
plane.  Depending  upon  the  deposition  technique,  the  layer  thickness  on 
these  edges  is  generally  less  than  that  on  the  top  or  bottom  surfaces  by 
several  hundred  angstroms.  Also,  because  of  the  rough  edges,  there  are 
areas  which  have  received  only  a  very  thin  layer  with  poor  adhesion  or  no 
layer  at  all.  Such  areas  have  a  major  effect  on  the  encapsulated  material 
during  -.nneals  above  800° C.  At  these  elevated  temperatures  phosphorus  vapor 


■■  Si3N4 


LP-2097 

Figure  5.8:  (a)  Represuntation  of  the  cleaved  edges  of  a  Si^N^  encapsulated 

InP  sample  prior  to  annealing,  (b)  Cut-away,  side  view  representation  of  the 
same  sample  following  a  high  temperature  anneal  (T  >  800°C) .  Very  rapid 
thermal  etching  due  to  preferential  P  vapor  loss  from  an  uncapped  edge  area 
has  resulted  in  undercutting  of  an  otherwise  intact  surface  layer  of  Si_N,  [87] 
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escapes  rapidly  through  thinnly  capped  or  uncapped  areas.  This  preferential 
P  evaporation  initiates  the  formation  and  agglomeration  of  In  droplets,  lead¬ 
ing  to  thermal  etching  [79-81]...  This -etching  process  of  the  In  droplets  trying 
to  maintain  equilibrium  saturation  with  phosphorus  by  reacting  with  the 
underlying  InP  continues  for  the  duration  of  the  anneal.  Of  interest  hers 
is  the  fact  that  once  this  thermal  etching  process  is  initiated  at  the  edges, 
it  will  continue  to  proceed  regardless  of  the  quality  or  integrity  of  the 
surface  encapsulating  layer.  For  anneals  performed  above  800° C,  Intact  sur¬ 
face  Si„N,  films  have  often  been  observed  which  are  undercut  by  thermal 
3  4 

etching.  This  situation  is  shown  schematically  in  Figure  5.8(b).  The  InP 
surface  not  affected  by  this  thermal  etching  appears  to  be  of  comparable 
quality  to  surfaces  annealed  at  lower  temperatures.  Therefore,  the  annealing 
temperature  limit  for  our  Si^N^  caps,  or  any  other  cap  capable  of  surviving 
these  elevated  temperatures,  depends  largely  on  the  sample  size  and  on  the 
amount  of  material  loss  around  the  edges. 

In  an  attempt  to  prevent  edge  losses  at  high  temperatures,  we  have 
examined  double  layered  caps  involving  SiC^/Si^N^  and  PSG/Si^N^.  Because  of 
the  simplicity  of  the  pyrolytic  deposition  method,  these  oxides  proved  to 
be  a  convenient  means  of  reducing  edge  loss  by  simply  increasing  the  total 
cap  thickness.  Our  experience  has  been  that  the  most  inert  encapsulant  at 
high  temperatures  (T  >  800°C)  resulting  in  the  best  InP  surfaces  is  a  double 
layered  cap  of  thick  PSG  over  our  standard  Si^N^  encapsulant.  Similar 
advantages  of  multilayer  encapsulation  have  been  discussed  by  others  for 
GaAs  [125]  and  InP  [82,126]. 
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6.  PHOTOLUMINES CNE CE  STUDIES  OF  IMPLANTED  InP 


6. 1  Introduction 

Successful  implantation  processing  of  InP  will  require  detailed 
information  regarding  the  thermal  annealing  parameters  of  temperature  and 
time  necessary  to  reorder  the  crystalline  lattice  damaged  by  the  atomic 
collisions  associated  with  implantation  and  to  activate  implanted  dopants. 

In  this  chapter,  photoluminescence  has  been  used  to  examine  the  optical 
spectra  of  semi-insulating  InP  processed  under  various  conditions  so  as 
to  characterize:  (1)  the  extent  of  lattice  reordering  which  occurs  during 
ishochronal  and  isothermal  anneals  of  InP  damaged  by  an  electrically 

4 

inactive  impurity  (  He);  and,  (2)  the  degree  of  lattice  reordering  and 

optical  activation  of  the  implanted  impurity  which  occurs  during  isochronal 
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and  isothermal  anneals  of  InP  implanted  with  an  acceptor  impurity  (  Be) . 

6.2  PL  Spectra  of  Virgin  SI  InP 

Before  discussing  the  spectra  from  implanted  and  annealed  InP  it 
would  be  useful  at  this  point  to  briefly  summarize  the  spectral  character¬ 
istics  of  the  virgin,  semi-insulating  InP  (Fe)  used  in  these  studies. 
Approximately  65  samples  from  two  neighboring  slices  of  the  same  single 
crystal  boule  were  characterized  prior  to  processing,  and  Figure  6.1  is  typical 
of  the  observed  spectra.  The  strongest  peak  in  the  spectra  occurs  at 
1.414  eV;  this  is  the  near  bandedge  luminescence  which  contains  a  large 
number  of  unresolved  emission  lines.  The  constituent  lines  of  this 
emission  peak  have  been  attributed  to  free  excitons,  excitons  bound  to 
neutral  and  ionized  donors,  and  to  excitons  bound  to  neutral  acceptors  [127-129] 
This  emission  at  1.414  eV  has  been  observed  in  the  spectra  of  InP  prepared 
under  almost  all  growth  techniques  [127-133]  and  is  most  often  attributed  to 
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Figure  6.1:  Photoluminescence  spectra  of  vLrgin  Fe-doped,  semi-insulating  (100) 
InP.  Numerical  factors  are  amplification  r  dative  to  that  used  for  the  bandedge 
luminescence  [751. 


excitons  bound  to  a  neutral  acceptor  [109].  The  width  (FWHM)  of  this 
1.414  eV  band  is  7.3  meV.  For  this  band  we  find  a  standard  deviation  of 
+  32%  in  the  values  of  the  peak  intensity  and  of  the  integrated  band 
intensity  for  the  65  samples  studied.  We  shall  refer  to  this  band  in  the 
PL  spectra  of  bulk,  SI  InP  as  the  bandedge  (BE)  peak. 

The  second  prominent  feature  in  Figure 6.1  is  located  at  1.384  eV. 
Singlet  or  doublet  emission  lines  have  been  noted  in  this  energy  range 
(1.380-1.384)  from  InP  grown  by  various  methods.  These  lines  have  been 
attributed  either  to  a  two-hole  transition  involving  an  exciton  bound  to 
the  neutral  acceptor  responsible  for  the  1.414  eV  transition  [127]  or  to 
band  to  acceptor  (BA)  and  donor-acceptor  (DA)  transitions  [108,128,133]. 
Recent  articles  identify  the  acceptor  involved  in  the  1.384  eV  peak  as 
carbon  [130,133,134].  The  width  of  this  band  in  our  samples  is  14.9  meV. 

We  find  for  this  band,  as  with  the  BE  band,  a  standard  deviation  of  +  33% 
in  the  values  of  the  peak  intensity  and  of  the  band  integrated  intensity. 

We  shall  refer  to  this  second  band  feature  in  the  PL  spectra  as  the  BA/DA 
band. 

The  third  band  in  Figure  6.1  is  located  at  1.341  eV.  This  feature 
is  separated  by  ~  43  meV  from  the  1.384  eV  band  and  therefore  is  attributed 
to  the  longitudinal  optical  phonon  replica  (1-LO)  of  the  1.384  eV  band  [131] 

Pomrenke  et  al.  [74]  have  noted  that  as  a  result  of  native 
defects  and/or  residual  impurities  the  resolution  of  peaks  in  LEC  InP(Fe) 
is  limited  to  those  separated  by  more  than  6.5  meV.  On  the  basis  of  the 
classical  Rayleigh  criteria  and  the  width  of  the  BE  band  observed  here,  we 
must  agree  with  this  observation.  It  should  be  noted,  however,  that 
shoulders  or  weak  peaks  superimposed  on  these  bands  which  were  separated  by 
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only  3-4  meV  were  often  observed.  Due  to  the  resolution  limit  and  small 
variations  in  peak  wavelengths,  the  identification  of  peak  energies  is 
confined  to  three  significant  figures. 

6.3  PL  Spectra  of  ^He  Implanted  SI  InP 

4 

In  this  section  the  lattice  reordering  in  51  keV,  He  implanted 

SI  InP  as  a  function  of  implant  f luence  and  anneal  temperature  and  time  is 

examined.  Helium  was  chosen  because  it  should  remain  electrically  and 

chemically  inactive  in  the  crystal,  thereby  affecting  the  lattice  mainly 

through  the  implantation  damage  it  creates.  For  the  energy  and  fluences 

examined  here  the  implant  energy  deposition  density  ranges  from~4.6  x  10^ 
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-  4.6  x  10  keV  cm  .  Atomic  profiles  (SIMS)  of  the  damage  related  re¬ 
distribution  of  ^Fe  in  these  samples  are  presented  in  Chapter  7.  Photo- 
luminescnece  is  useful  in  characterizing  lattice  reordering,  since  in 
incompletely  annealed  samples  nonradiative  channels  compete  with  radiative 
channels  for  electron-hole  recombination  thereby  decreasing  total  sample 
luminosity,  and  also  because  the  bandedge  luminosity  associated  with 
excitons  is  extremely  sensitive  to  lattice  perfection. 

Representative  spectra  are  shown  in  Figure  6.2  for  samples  implanted 
15  -2  4 

with  51  keV,  10  cm  He  and  annealed  for  30  min  at  the  indicated  tempera¬ 
tures.  Several  spectral  features  are  related  to  annealing  with  a  Si^N^ 
cap,  as  discussed  in  the  Section  4.4.2.  In  Figure  6.2(b)  the  twin  peaks 
at  1.381  eV  and  1.378  eV  which  occurred  in  annealed  unimplanted  samples  are 
again  observed.  Previously  this  pair  was  seen  only  in  samples  annealed 

below  650°C,  and  the  1.378  eV  peak  dominated  the  BA/DA  band  for  samples 
o  A 

annealed  above  650  C.  In  this  He  implanted  series  the  existence  of  both 
peaks  at  slightly  higher  anneal  temperatures  probably  results  from  the 
fact  that  a  shorter  anneal  time  (30  min)  was  used.  As  expected  from  the 
Si^N^  capping  results,  at  800°C  the  1.378  eV  line  dominates  the  emission  of 


Energy  (eV) 

L28  1.32  1.36  1.40  1.44 


9800  9600  9400  9200  9000  8800  8600 

Wavelength  (A)  LP-zoro  ^  ^ 

Figure  6.2:  Photo luminescence  spectra  of  SI  InP  implanted  with  10  cm  , 

51  keV  Sle  and  annealed  for  30  min  at  (a)  800°C,  (b)  700°C,  and  (c)  550°C. 
Numerical  factors  are  amplification  relative  to  bandedge  luminescence  of 
Figure  6.1  [75]. 
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the  BA/DA  band  in  Figure  6.2(a).  Although  not  shown  here,  the  behavior  of 

13  -2  14  —2  4 

the  BA/DA  band  in  the  10  cm  and  10  cm  He  implanted  and  annealed 
samples  is  similar,  i.e.,  for  T  <_  700°C  the  BA/ DA  band  contains  two  peaks 
of  comparable  intensity  at  1.381  eV  and  1.378  eV  and  for  T  >  700°C  the 
band  emission  is  dominated  by  the  1.378  eV  line.  In  Figure  6.2(a)  and 
6.2(b)  we  can  see  the  increased  luminosity  between  the  BA/DA  and  1-LO  bands 
which  was  already  noted  in  Section  4.4.2.  At  550°C  [Figure  6.2(c)]  a 
broad  feature  at  1.363  eV  appears  to  be  part  of  this  increased  luminescence 
between  the  BA/ DA  and  1-LO  bands.  Some  of  our  results  suggest  this 
feature  may  be  related  to  implant  dose.  From  Figure  6.2  the  recovery  of 
the  lattice  from  implantation  damage  can  be  seen  in  the  return  of  the  BE 
band  width  to  preimplant  values  and  also  in  the  decreasing  amplification 
required  to  observe  signal  emission  as  the  anneal  temperature  increases. 

To  monitor  this  crystalline  reordering  as  a  function  of  implant 
dose  and  anneal  temperature,  we  have  plotted  the  integrated  intensity  of 
the  BE  band  for  a  series  of  isochronal  anneals  in  Figure  6.3.  Annealing 
stages  have  been  eported  to  exist  in  implanted  InP  at  ~  200°C  and 
~  400°C  on  the  basis  of  Hall  measurements  [18].  Using  photoluminescence, 
we  have  not  been  able  to  study  samples  annealed  much  below  500°C  because 
of  low  luminosities  in  the  heavily  damaged  material.  It  is  generally 
recognized,  however,  that  in  device  fabrication,  anneals  in  recess  of  500°C 
are  required  to  activate  implanted  impurities.  Therefore,  the  measurements 
reported  here  span  the  anneal  temperature  range  of  primary  interest.  In 
agreement  with  reports  by  Davies,  jet  al.  [18,76],  we  note  rapid  recovery 
in  BE  band  emission  between  550°C  and  650°C  (Figure  6.3(a)).  At  800°C, 
all  doses  shown  BE  integrated  intensitites  above  unimplanted  levels.  As 
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Figure  6.3:  Integrated  BE  band  intensities  for  a  series  of  isochronal 
(30  min)  anneals  of  51  keV,  ^He  implanted  SI  InP  [75]. 


Relative  Integrated  Intensity 


Figure  6.4:  Integrated  BA/DA  band  intensities  for  a  series  of  isochronal 
(30  min)  anneals  of  51  keV,  ^He  Implanted  SI  InP  [75]. 
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discussed  in  Section  4.4.2,  the  increase  in  the  BE  band  integrated  intensity 
does  not  occur  until  the  1.378  eV  line  dominates  the  emission  of  the  BA/DA 
band.  An  accurate  assessment  of  the  annealing  process  is  hindered  by  the 
general  rise  in  the  intensities  above  preimplant  levels.  With  this 
qualification  in  mind,  our  data  indicate  continuing  lattice  recovery  from 
moderate  implant  damage  to  the  highest  temperature  of  ~  750°C. 

For  completeness,  we  exhibit  in  Figure  6.4  the  isochronal  annealing 
characteristics  of  the  BA/DA  band  in  these  ^He  implanted  samples.  The 
effect  of  the  tail  of  luminescence  on  the  Ioj  energy  side  of  the  BA/DA  band 
has  not  been  included  in  the  integration  of  this  band's  intensity.  As 
was  the  case  with  the  Si^N^  capping  series  of  Section  4.4.2,  the  increase 
of  the  BA/DA  band  integrated  intensity  above  virgin  levels  occur  at  a 
lower  temperature  than  does  the  comparable  increase  in  BE  band  intensity. 

At  levels  above  those  of  unimplanted  samples,  the  BA/DA  band  is  dominated  by 
emission  at  1.378  eV.  Because  of  the  similarities  of  Figures  6.3  and  6.4 
to  Figure  4.8  and  since  He  is  not  expected  to  be  optically  active  we  may 
conclude  that  the  high  temperature  anneals  (T  >  700°C)  of  this  series  are 

influenced  by  factors  related  to  annealing  with  a  Si^N^  cap. 

o  15  -2 

The  isothermal  (750  C)  annealing  characteristics  of  10  cm  , 

4 

He  implanted  SI  InP  are  shown  in  Figure  6.5  as  a  function  of  anneal  time. 

It  is  clear  that  at  this  temperature  the  integrated  intensity  of  the  BA/DA 
band  is  above  virgin  levels  for  all  t  >_  15  min.  In  the  15  and  30  min 
anneals  the  band  contains  1.381  eV  and  1.378  eV  lines  of  comparable 
intensities  while  the  1.378  eV  line  dominates  in  the  60  min  and  120  min 
anneals.  This  observation  supports  our  previous  assumption  that  not  only 
anneal  temperature  but  also  anneal  time  influences  the  point  at  which  the 
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Figure  6.5:  Integrated  intensities  of  (a)  the  BE  band  and  (b)  the  BA/DA 
band  for  a  series  of  isothermal  (750°C)  anneals  of  51  keV,  10^  cm  ^  ^He 
implanted  SI  InP  [75]. 
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1.378  eV  line  dominates  the  BA/DA  emission  band.  Comparing  Figure  6.5(a) 

to  6.5(b)  it  becomes  evident  that  the  BE  band  intensity  rises  above  virgin 

levels  at  later  anneal  times  than  does  the  BA/DA  band.  This  result  parallels 

the  relationship  of  these  two  bands  in  regard  to  the  temperature  at  which 

their  integrated  intensities  exceed  unimplanted  levels.  Figures  6.5(a) 

and  6.5(b)  indicate  no  saturation  of  the  integrated  intensities  with  time; 

instead,  the  intensities  of  both  band  increase  exponentially  with  anneal 

time.  This  result  together  with  those  of  Section  4.4.2  suggest  the  diffusion 

of  an  optically  active  impurity  into  the  analyzed  region  as  a  result  of 

annealing  with  a  Si^N^  cap. 

9 

6.4  PL  Spectra  of  Be  Implanted  SI  InP 

In  this  section,  lattice  reordering  and  impurity  activation  in 

9 

SI  InP  implanted  with  100  keV  Be  is  examined.  Several  potential  InP 

device  applications  will  require  a  controllable  method  of  acceptor  doping. 

9 

Be  is  an  attractive  implanted  acceptor  to  investigate  because  its  low 
mass  results  in  less  implant  damage  than  with  heavier  ions.  With  the 
energy  and  doses  examined  here  the  implant  energy  deposition  density  ranges 

iq  20  _  ■» 

from  ~  2  x  10  -  2  x  10  keV  cm  .  The  atomic  (SIMS)  profiles  of  the  Be 

and  Fe  redistribution  in  these  samples  are  presented  in  Chapter  7. 

15  -2  13  -2  9 

Representative  spectra  from  10  cm  and  10  cm  Be  implanted 

samples  are  shown  in  Figure  6.6.  Two  new  features  appear  which  have  not 

4 

been  seen  in  previous  Si^N^  capped  or  He  implanted  spectra.  Comparing 
Figure  6.1  to  Figure  6.6(a)  it  is  evident  that  the  peak  of  the  BE  band 


shifts  to  lower  energies.  This  BE  peak  shift  to  lower  energies  increases 

with  increasing  implant  dose.  After  an  800°C  anneal  the  BE  peak  position  for 

13  -2  14 

the  fluences  examined  here  is:  1.413  eV  for  10  cm  ,  1.410  eV  for  10 
_2  15  -2 

cm  ,  and  1.408  eV  for  10  cm  .  Since  previously  examined  spectra  showed 


Relative  Intensity  (Arbitrary  Units) 
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no  detectable  shift  In  the  BE  peak  position,  this  behavior  cannot  be 
attributed  to  capping  or  damage  effects.  We  believe  that  this  dose 
dependent  peak  shift  in  the  BE  band  most  probably  results  from  inter-site 
interactions  (Stark  interaction  effects  or  impurity  banding  effects). 

Such  effects  would  be  consistent  with  the  dose  and  anneal  temperature 
dependence  of  the  shift  since  they  are  related  to  the  concentration  of 
activated  Be  acceptors.  Considering  the  location  of  the  new  BE  peaks  in 
the  band  gap  of  InP  it  would  be  likely  that  the  peak  transitions  result 
from  exciton  recombination  at  Be  acceptors. 

9 

The  other  new  feature  in  the  spectra  of  these  Be  implanted  samples 

is  the  peak  position  of  the  BA/DA  band.  Within  our  resolution  limits  the 

position  of  this  peak  is  1.382  eV  for  all  doses  examined  here.  We  have 

15  -2  9 

surveyed  a  limited  number  of  10  cm  Be  implants  into  undoped  VPE  InP 
in  which  the  peak  resolution  is  much  better  than  in  the  SI  InP.  These  VPE 
samples  are  characterized  by  a  peak  at  1.383  eV  which  is  not  present  in 
unimplanted  material.  Additionally,  we  have  examined  the  spectra  at  various 
depths  extending  from  the  surface  to  0.75  yra  in  samples  implanted  with  10^ 
cm  ^  ^Be  and  annealed  at  750°C  for  30  min.  The  resulting  optical  profile 
of  the  BA/DA  band  emission  at  1.382  eV  shows  some  correlation  to  the  atomic 

9 

Be  profiles  obtained  with  SIMS  [138].  Also,  from  Figure  6.6  and  other 

spectra  it  is  clear  that  the  peak  intensity  of  the  1.382  eV  emission  increases 

13  15  -2 

with  implant  dose  in  the  range  from  10  -10  cm  .  Therefore,  we 

9 

believe  that  the  1.382  peak  in  Be  implanted  SI  InP  is  due  to  optically 
active  Be  acceptors.  An  estimate  of  the  Be  acceptor  ionization  energy  may 
be  obtained  if  we  ignore  possible  inter-site  effects  and  assume  the  1.382 
eV  emission  results  from  a  band  to  acceptor  transition.  Using  a  band  gap 
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energy  ■  1.4233  eV  (2°K)  [139]  we  estimate  the  Be  acceptor  ionization 
energy  to  be  41.3  +  3  meV.  This  estimate  is  within  the  range  of  reported 
binding  energies  for  other  acceptors  in  InP  [74,109,130,134]. 

Other  features  seen  in  Figure  6.6  were  discussed  above.  At 
550°C  the  10^  cm  ^  ^Be  sample  has  a  BA/ DA  band  which  contains  two  peaks, 
one  of  which  (1.378  eV)  has  been  shown  to  be  associated  with  Si^N^  capping 
and  annealing.  The  10^  cm  ^  ®Be  sample  annealed  at  550°C  does  not  show 
this  1.378  eV  peak,  presumably  because  of  the  competition  from  the  emerging 
1.382  eV  peak  associated  with  Be  acceptors.  From  Figure  6.6(c)  we  can  also 
see  increased  luminescence  between  the  BA/DA  band  and  the  1-LO  phonon  band 
as  in  previous  spectra. 

In  Figure  6.7  we  have  plotted  the  integrated  intensity  of  the  BE 

band  as  a  function  of  implant  dose  and  anneal  temperature.  Note  that  the 

10^  cm  ^  series  recovers  more  rapidly  in  the  range  from  550°C  -  650°C  than 
14  -2  15  -2 

do  either  10  cm  or  10  cm  dose  samples.  This  variation  in  annealing 

13  -2 

behavior  between  the  10  cm  dose  and  the  higher  doses  seems  to  be  much 

larger  than  might  be  expected  from  the  results  of  Figure  6.3.  There  the 

plots  showed  nearly  parallel  curves  that  were  predictably  dependent  on 

fluence.  Our  estimates  of  the  energy  deposition  density  indicate  an 

13  -2  14  -2  9 

overlap  of  this  quantity  for  the  10  cm  and  10  cm  implants  of  Be 

with  the  10^  cm  ^  and  10^  cm  ^  implants  of  Sie  (2  x  10^  -  4.6  x  10^  keV 

-3  9 

cm  ).  Although  the  differences  in  annealing  behavior  for  the  Be  doses 

considered  here  are  not  presently  understood,  it  seems  that  they  do  not 

result  from  different  levels  of  implant  damage  density. 
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As  might  be  expected  from  our  previous  isochronal  anneal  plots, 
for  T  _>  750°C  we  find  that  the  BE  band  emission  is  greater  than  in  unimplanted 
samples.  Because  of  the  strong  luminescence  associated  with  Be  acceptors 
at  1.382  eV,  the  1.378  eV  line  does  not  dominate  the  BA/DA  band  in  these 
samples  as  was  the  case  in  previous  samples  where  the  BE  band  intensity  was 
above  virgin  levels.  At  low  annealing  temperatures,  however,  (before  sub¬ 
stantial  concentrations  of  Be  are  activated)  the  1.378  eV  line  can  be 
observed  (Figure  6.6(c)).  As  stated  earlier,  we  believe  that  the  emergence 
of  the  emission  at  1.378  eV  is  related  to  the  observation  of  integrated 
band  intensitites  greater  than  those  in  unprocessed  material.  In  recent 
years  articles  have  appeared  which  report  unusual  surface  properties  for  InP 
with  high  luminosities  associated  with  minimal  surface  band  bending  [140*141]. 
It  is  possible  that  the  increased  levels  of  BE  band  luminescence  we  con¬ 
sistently  observe  at  high  anneal  temperatures  result  from  a  reduction  in  the 
surface  space  charge  layer  of  SI  InP  due  to  the  accumulation  of  the  impurity/ 
defect  involved  in  the  1.378  eV  line.  In  our  optical  depth  profiling  of  the 
Be  acceptor  emission  we  have  observed  that  the  BE  band  emission  is  above 
virgin  levels  at  all  depths  examined  from  the  surface  to  7500  X;  therefore, 
the  phenomena  is  not  restricted  to  the  near  surface.  Based  on  preliminary 

SIMS  results  this  depth  dependence  may  be  related  to  the  Increased  levels 
28 

of  Si  we  have  found  in  Si,jN^  capped  and  annealed  samples.  Further  work 
is  being  done  to  determine  if  such  a  relationship  does  exist.  A  complicating 
factor  in  the  space  charge  layer  explanation  for  increased  luminosity  levels 
is  the  finding  that  the  compensating  deep  acceptor  56Fe  getters  to  the  surface 
and  into  the  implanted  region  above  background  levels  [138].  Therefore,  it 
is  possible  that  the  BE  enhancement  results  from  the  annealing  of  native 
defects  in  bulk  InP  or  from  a  gettered  impurity  with  greater  radiative  ef- 


ficiency  than  that  of  the  other  impurities  in  virgin  material.  Obviously 
more  research  will  be  required  to  ascertain  the  origin  of  the  increased 
levels  of  luminosity  in  high  temperature  annealed  SI  InP. 

Figure  6.8  presents  the  isochronal  annealing  behavior  of  the 

BA/DA  band  in  these  Be  implanted  samples.  Again,  the  increased  luminescence 

on  the  low  energy  side  of  the  BA/DA  band  has  not  been  included  in  the 

integration  of  the  band's  intensity  .  For  purposes  of  comparison  to  Figure 

6.4  we  have  included  the  unimplanted  reference  level.  Because  the  BA/DA 

band  is  now  dominated  by  emission  associated  with  Be  acceptors,  the  curves 

of  Figure  6.8  differ  significantly  from  those  of  Figure  6.4  which  we  have 

proposed  to  be  influenced  by  near  surface  gettering  of  an  impurity  during 

annealing.  For  low  and  moderate  doses  it  appears  from  Figure  6.8  that 

saturation  of  the  Be  activation  occurs  in  the  range  700  _<  T  <_  750°C.  For 
15  -2 

10  cm  doses,  activation  of  Be  appears  to  increase  smoothly  with  temper¬ 
ature  up  to  750°C.  Between  750°C  and  800°C  all  three  fluences  exhibit  a 
small  decrease  in  Be  related  emission,  which  may  result  from  Be  migration 
away  from  the  near  surface  region. 

15  -2 

The  isothermal  (750°C)  annealing  characteristics  of  10  cm  , 

q 

Be  implanted  SI  InP  are  shown  in  Figure  6.9.  Unlike  Figure  6.5(b)  in 
which  the  BA/DA  band  integrated  intensity  increased  exponentially  with 
anneal  time,  it  is  clear  from  Figure  6.9(b)  that  the  BA/DA  band,  now 
dominated  by  Be  emission,  shows  saturation  in  the  time  range  from  30-60  min. 
For  anneal  times  t  >  60  min  a  decrease  in  Be  related  emission  is  observed, 
which  would  be  consistent  with  the  assumption  of  Be  migration  away  from  the 


near  surface  region. 
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Figure  6.8:  Integrated  BA/DA  band  intensities  for  a  series  of  isochronal 

g 

(30  min)  anneals  of  100  keV,  Be  implanted  SI  InP  [75]. 
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Figure  6.9:  Integrated  intensities  of  (a)  the  BE  band  and  (b)  the  BA/DA 
band  for  a  series  of  isothermal  (750°C)  anneals  of  100  keV,  10^  cm  ^  ^B< 
implanted  SI  InP  [75]. 


6.5  Summary 

In  this  study  of  PL  spectra  from  1.44  -  1.30  eV  we  have 

observed  the  following: 

(1)  The  virgin,  (100)  Fe-doped  InP  examined  here  is  characterized  by 
three  luminescence  bands  at  1.414,  1.384,  and  1.341  eV.  These  are 
attributed  respectively  to  bandedge  luminescence  (BE) ,  band  to 
acceptor  and  donor  to  acceptor  (BA/DA)  luminescence,  and  luminescence 
from  the  1-LO  phonon  replica  of  the  BA/DA  band.  The  widths  (FWHM)  of 
the  BE  and  BA/DA  bands  are  7.3  and  14.9  meV  respectively.  For  the 

~  65  samples  taken  from  two  neighboring  slices  of  the  same  boule, 
variations  of  ~  +  30%  were  noted  in  the  peak  intensity  and  integrated 
intensity  of  the  BE  and  BA/DA  bands. 

(2)  An  increase  in  luminosity  between  the  BA/DA  and  1-L0  band  has  been  ob¬ 
served  after  annealing  in  both  unimplanted  and  implanted  samples.  Weak 
structures  at  ~  1.36  eV  are  often  superimposed  on  this  broad  tail  of 
luminescence.  The  origin  of  this  luminescence  is  not  presently  under¬ 
stood,  but  possibilities  were  discussed  previously  in  Section  4.4.2. 

(3)  High  temperature  (T  >  700°C)  anneals  of  unimplanted  SI  InP  or  material 

implanted  with  inactive  impurities  result  in  shifting  of  the  BA/DA  peak 

from  1.384  eV  to  1.378  eV.  Preliminary  analysis  using  SIMS  indicates 

o  28 

that  in  the  Si^N^  capped  and  750  C  annealed  samples  the  Si  concentra¬ 
tion  level  is  above  virgin  levels  to  a  depth  of  ~  7000  X.  Besides  Si 
in-diffusion  from  the  cap,  annealing-induced  radiative  defects  or 
gettering  of  impurities  from  he  bulk  juld  cause  the  1.378  eV  peak. 

(4)  In  all  samples  examined,  high  temperature  (T  >  750°C)  anneals  result  in 


integrated  BE  band  intensities  greater  than  those  in  unimplanted  samples 
This  phenomenon  has  been  correlated  to  the  presence  of  the  1.378  eV 


peak.  The  luminosity  increase  may  result  from  a  reduction  in  the  sur¬ 
face  space  charge  layer  of  SI  InP  due  to  the  accumulation  of  the  active 
impurity  involved  in  the  1.378  eV  transition.  It  is  also  possible  the 
8G  band  enhancement  results  from  the  annealing  of  native  defects  in 
bulk  InP  or  from  a  gettered  impurity  with  greater  radiative  efficiency 
than  that  of  the  other  impurities  in  virgin  material. 

(5)  Isothermal  (750°C)  anneals  of  51  keV,  10*^  cm  2  Sie  implanted  SI  InP 

exhibit  exponential  increases  above  virgin  levels  with  anneal  time 

for  both  the  BE  and  BA/DA  band  integrated  intensities.  No  saturation 

effects  are  observed  for  t  _<  120  min.  We  believe  this  supports  the 

view  of  an  active  impurity  diffusing  into  the  analysis  region.  Fol- 
1.3  15  2  ^ 

lowing  51  keV,  10  -  10  cm  implants  of  He,  thirty-minute  anneals 

at  750°C  appear  sufficient  to  restore  the  lattice  to  approximately 

preimplant  levels  of  quality. 

q 

(6)  100  keV  implants  of  Be  result  in  dose- dependent  shifting  of  the  BE 
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band  peak  from  1.414  eV  in  unimplanted  samples  to  1.408  eV  in  10  cm 
implanted  InP.  This  shifting  is  probably  due  to  inter-site  interactions. 

Q 

In  the  Be  implanted  samples  the  BA/DA  band  is  dominated  by  a  new 
emission  peak  at  1.382  eV  which  we  show  to  be  associated  with  Be.  We 
estimate  the  Be  acceptor  ionization  energy  to  be  41.3  +  3  meV.  Maximum 
Be  activation  appears  to  occur  for  all  doses  at  ~  750°C  for  30  min 
anneals.  At  750°C,  10^  cm  2  ^Be  implanted  SI  InP  exhibits  saturation 
of  Be  activation  between  30-60  min  with  decreasing  Be  associated  emission 
at  1.382  eV  resulting  from  longer  anneals. 


7.  SIMS  STUDIES  OF  ATOMIC  REDISTRIBUTION  IN  ANNEALED  InP 


In  this  chapter  the  redistribution  of  intentionally  present, 
compensating  dopants  (Fe  or  Cr)  and  of  implanted  impurities  (Be,  Mg>  and  Si) 
in  annealed,  semi-insulating  InP  has  been  studied  using  SIMS. 

7.1  Fe  and  Cr  Redistribution  in  SI  InP 

Semi- insula ting  InP  is  extremely  important  in  the  fabrication  of 
devices  because  it  can  serve  as  a  substrate  for  the  subsequent  growth  of  an 
epitaxial  device  layer,  and  it  can  provide  device  isolation  when  implanta¬ 
tion  directly  into  the  material  is  employed.  Whereas  presently  almost  all 
available  semi-insulating  GaAs  is  chromium-doped,  semi-insulating  InP  can 
be  obtained  by  doping  either  with  chromium  or  iron.  Cr  has  a  much  lower 
solid  solubility  in  InP  than  in  GaAs,  which  necessitates  the  use  of  high 
purity  starting  material  to  achieve  high  resistivity  [16].  Since  such  purity 
in  InP  is  still  not  routinely  obtained,  Fe,  which  has  a  higher  solid  solu¬ 
bility  in  InP  than  does  Cr,  is  being  used  as  an  alternative  in  the  growth  of 
semi- insulating  InP  [142,143].  Recently,  a  number  of  studies  have  shown  the 
redistribution  of  Cr  in  annealed  GaAs  [144-146]  and  in  implanted  GaAs  [58,147]. 
The  redistribution  of  the  bulk  compensating  dopants,  Fe  or  Cr,  in  annealed 
and  implanted  semi-insulating  InP  is  described  in  the  following  section. 

7.1.1  Redistribution  Resulting  from  Annealing 

Figures  7.1  and  7.2  exhibit  the  profiles  of  Fe  and  Cr  resulting 
from  annealing  the  substrates  with  a  Si^N^  cap.  With  the  exception  of  a 
large  Fe  accumulation  within  the  first  1200  X  of  the  surface,  the  Fe  profile 
in  the  550°C,  30  min  anneal  Is  otherwise  indistinguishable  from  the  unannealed 
profile.  This  is  not  true  of  the  profile  resulting  after  a  800°C,  30  min 
anneal.  Again,  an  Fe  accumulation  region  with  a  maximum  concentration  of 


Relative  Intensity 


LO  2.0 

Depth  (/xm) 

'  LP-1693 


Figure  7.1:  SIMS  profiles  of  Fe  in  a  semi-insulating  InP  substrate  following 
30  min.  anneals  at  550°C  and  800°C  with  a  Si-N  encapsulant  [21]. 
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~  4  x  10  cm  extends  from  the  surface  to  a  depth  of  ~  1000  X.  However, 

following  this  zone  is  an  Fe  depletion  region.  This  broad,  flat-bottomed 

region  originates  within  1000  X  of  the  surface  and  returns  to  the  normal 

background  level  at  a  depth  of  ~  1.2  ym.  The  minimum  Fe  concentration 

16  —3 

within  this  region  is  ~  1.4x10  cm 

The  redistribution  of  Cr  resulting  from  annealing  differs  signifi¬ 
cantly  from  that  of  Fe.  As  with  Fe  there  is  no  evidence  of  major  redistri¬ 
bution  resulting  from  a  500°C,  30  min  anneal.  Qualitatively,  the  Cr  profile 
after  a  800°C,  30  min  anneal  is  similar  to  that  of  Fe,  in  that  a  surface 
accumulation  region  and  a  depletion  zone  appear.  The  surface  accumulation 
zone  ends  at  a  depth  of  ~  1000  X,  while  the  Cr  depletion  region  extends  to 

a  depth  of  ~  8000  X.  The  minimum  Cr  concentration  within  this  region  is 
14  -3 

~  9.2x10  cm  .  Although  redistribution  of  Cr  occurs  over  a  region  of 

~  8000  X  versus  "  12,000  X  for  Fe,  Cr  redistribution  is  more  severe  in  the 
sense  that  the  concentration  of  the  compensating  dopant  is  an  order  of 
magnitude  lower  in  the  depletion  zone  than  in  Fe-doped  substrates. 

The  depletion  regions  observed  in  the  Fe  and  Cr  doped  substrates 
are  in  keeping  with  reports  of  conductive  surface  layers  following  thermal 
annealing  of  InP  substrates  [19,76].  Additionally,  the  redistribution  is 
very  similar  to  recent  reports  of  Cr  migration  in  annealed  GaAs  substrates 
[58,144-146].  Originally,  the  accumulation  and  subsequent  depletion  regions  of 
Cr  in  GaAs  were  speculated  to  be  due  to  a  strain-induced  effect  caused  by  the 
thermal  mismatch  between  the  encapsulant  and  the  GaAs  [147].  However,  recent 

work  by  Vasudev  et  al. [58] ,  indicates  that  Cr  redistribution  in  GaAs  is  not  sig 
nificantly  influenced  by  the  encapsulant  and  is  in  fact  larger  under  capless 
anneal  conditions.  Moreover,  Kasahara  and  Watanabe  [146]  reported  that  Cr 
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redistribution  in  GaAs  is  controlled  by  the  partial  pressure  of  arsenic 
during  the  anneal,  which  confirms  an  early  observation  by  Tuck  et  al.  [144], 
concerning  the  role  of  As  overpressure  with  relation  to  Cr  redistribution. 

At  present,  therefore,  the  best  explanation  of  Cr  redistribution  in  GaAs 
appears  to  be  that  it  results  from  the  gettering  of  Cr  to  surface  regions 
which  have  high  concentrations  of  defects  introduced  either  by  annealing  or 
by  mechanical  damage  [148].  There  is  reason  to  believe  that  the  same  relation¬ 
ship  between  defects  and  compensating  dopant  also  explains  our  results  in  InP. 
It  is  unclear  at  present  what  importance  should  be  attached  to  the  quantitative 
differences  observed  between  Fe  and  Cr  redistribution  profiles  in  InP. 

Possibly  the  differences  can  be  explained  in  terms  of  the  surface  solubility 
and  original  background  density  of  the  compensating  dopant,  as  suggested  by 
Vasudev  et  al.  [58].  Furthermore,  differences  between  Cr  and  Fe  outdif fusion 
into  the  SijN^  cap  must  be  considered.  Obviously  further  work  is  needed  to 
determine  if  Important  differences  do  Indeed  exist  between  Fe  and  Cr  redistri¬ 
bution  in  InP  as  our  results  suggest. 

7.1.2  Redistribution  Resulting  from  the  Implantation  of  Neutral  Dopants 
Impurity  gettering  to  damaged  (i.e.,  vacancy-rich)  regions  produced 
by  ion  implantation  is  often  observed  in  Si  [149,150]  and  in  GaAs  [58,147,151]. 
In  this  section  the  gettering  of  Fe  and  Cr  to  damaged  regions  produced  by 
the  implantation  of  electrically  neutral  species  in  semi-insulating  InP  is 
discussed.  Helium  was  chosen  as  one  such  species  because  it  should  remain 
electrically  and  chemically  inactive  in  the  crystal,  thereby  having  no  effect 
on  the  compensating  dopants  except  form  the  damage  it  creates.  It  is  possible 
for  "gas  bubbles"  or  voids  to  form  by  the  diffusion  of  insoluble  gas  atoms  to 
nucleation  sites  [28,152].  To  rule  out  redistribution  due  to  the  possible 
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Figure  7.3:  SIMS  profiles  of  Fe  in  a  semi-insulating  InP  substrate  implanted 
15  -2 

with  He  (1.06  x  10  cm  ,  51  keV)  and  annealed  for  30  min.  at  the  indicated 
temperatures  with  a  Si^N^  encapsuiant  [21]. 
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Figure  7.4:  SIMS  profiles  of  Cr  in  a  semi-insulating  InP  substrate  implanted 
13  -2 

with  He  (1.06  x  10  cm  ,  51  keV)  and  annealed  for  30  min.  at  the  indicated 


temperatures  with  a  Si^N,.  encapsuiant  [21] 


presence  of  such  voids,  a  non-gaseous  element,  boron,  was  also  implanted. 

It  is  assumed  that  boron  fills  column  III  sites,  and  is  therefore  electrically 
inactive  in  the  crystal  (although  the  possibility  of  chemical  activity  cannot 
be  completely  ruled  out) . 

15  -2 

A  He  fluence  of  1.06  x  10  cm  was  implanted  into  the  samples 
at  51  keV  and  0.5  pA  scanned  current.  Under  these  conditions  the  LSS  projected 
range  is  3640  X  with  the  profile  of  energy  deposited  into  nuclear  collisions 
peaking  at  ~  3200  X  (see  Figure  2.2).  Figures  7.3  and  7.4  show  the  Fe  and 
Cr  profiles,  respectively,  resulting  from  thirty  minute  anneals  of  this  He- 
implanted  InP.  As  in  the  previously  discussed  unimplanted  samples,  a  surface 
accumulation  zone  exists  within  the  first  1000  X  for  both  types  of  compen¬ 
sating  dopants.  Surface  accumulation  of  Fe  is  much  more  pronounced  than  that 
of  Cr,  which  is  also  consistent  with  the  previous  results.  However,  unlike 
the  profiles  in  the  unimplanted  samples,  we  find  in  the  He  implanted  samples 
an  accumulation  region  followed  by  a  depletion  region  rather  than  broad 
depletion  regions  extending  to  ~  1  pm.  It  is  clear  that  both  Fe  and  Cr 
accumulate  in  a  region  about  the  projected  range  of  He.  The  550°C,  30  min 
anneal  profiles  show  little  restructuring  except  for  small  depletion  regions 
on  the  bulk  side  of  the  projected  range  peak.  The  750°C,  30  min  annealed 
Fe-doped  sample  has  an  accumulation  zone  with  a  concentration  of  '*•  2.6  x 
10^  cm  ^  and  a  depletion  zone  with  a  concentration  of  ~  3.4  x  10*^  cm 
The  Cr-doped  sample  annealed  under  identical  conditions  shows  even  larger 
variations  of  concentration  in  the  various  regions. 

Boron  was  implanted  at  126  keV  and  0.5  pA  to  a  fluence  of  6.8  x 
10x  cm  in  Fe-doped  InP  substrates.  The  LSS  projected  range  is  3640  A  and 
the  damage  profile  peaks  at  ~  2890  X  (see  Figure  2.2).  The  as-implanted 
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Figure  7.5:  SIMS  profiles  of  Fe  in  a  semi-insulating  InP  substrate  implanted 
14  -2 

with  B  (6.8  x  10  cm  ,  126  keV)  and  annealed  for  30  min.  at  the  indicated 
temperatures  with  a  Si_N.  encapsulant  [21], 
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boron  profile  obtained  by  SIMS  is  shown  in  Figure  7.5.  As  in  the  550°C 

profiles  of  He-implanted  Fe-doped  InP,  an  Fe  depletion  zone  is  formed  on  the 

bulk  side  of  the  projected  range  peak  with  a  minimum  Fe  concentration  of 
16  -3 

~  1.6  x  10  cm  .  Accumulation  of  Fe  into  the  implant  region  and  toward 
the  surface  is  evident  in  the  650°C  and  750°C  anneals.  The  qualitative 
similarities  between  the  redistribution  of  Fe  in  the  He  and  B  implanted 
samples  appears  to  indicate  that  we  are  observing  the  gettering  of  the  compen¬ 
sating  dopant  by  implant  related  defects  rather  than  diffusion  due  to  gas 
voids  or  chemical  effects'. 

9 

7.2  Redistribution  of  Implanted  Be  in  SI  InP 

There  exists  a  need  for  controllable  acceptor  doping  of  InP  in 
several  of  the  applications  mentioned  above.  For  example,  InP  p-type  layers 
are  used  as  confining  layers  in  typical  InGaAsP/InP  double-heterostructure 
diode  lasers.  Low  dark  current  avalanche  photodiodes  can  be  constructed 
from  structures  incorporating  an  InP  p-n  junction  above  a  detecting  InGaAsP 
layer  [153].  It  is  also  assumed  that  acceptor  doping  will  play  some  part  in 
the  fabrication  of  logic  circuits.  To  date  Zn  and  Cd  have  been  the  dominant 
impurities  used  to  obtain  p-type  InP  by  growth  or  diffusion  techniques.  It 
has  been  observed,  however,  that  both  impurities  are  rapid  dlffusants  in  InP 
and  hence  difficult  to  accurately  control  in  the  fabrication  of  p-n  junc¬ 
tions  [154-157].  Because  of  these  difficulties  the  implantation  of  acceptors 

9 

into  InP  has  been  considered  as  an  alternative  doping  technique.  Be  is  a 

particularly  appealing  ion  for  such  p-type  implants.  Due  to  its  low  mass, 

9 

Be  can  be  implanted  over  a  wide  range  of  depths  with  conventional  implanters, 
and  creates  less  lattice  damage  than  heavier  p-type  ions.  Beryllium  has 
already  proven  to  be  a  useful  implanted  acceptor  with  good  electrical  activa- 

9 

tion  in  GaAs  and  GaAsn  ,Prt  .  [135-137,158].  More  recently.  Be  implantation 

U • D  Ut  4 
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has  been  successfully  employed  to  fabricate  diodes  in  InP  and  in  InGaAsP 

o 

[159-162].  For  these  reasons  the  redistribution  of  implanted  Be  resulting 

from  furnace  anneals  was  studied  in  semi-insulating  InP. 

7.2.1  Comments  on  Mass  Interferences 

On  several  occasions  we  have  observed  molecular  interference  with 

the  ^Fe  signal  in  the  implanted  region.  This  interference  was  observed  in 
15  -2  9 

the  unannealed  10  cm  Be  implant  sample,  and  is  a  problem  only  for  the 

high  fluence  implants.  Since  the  mass  56  profile  reflected  the  peak  shape 
9 

of  the  unannealed  Be  implant,  we  presume  this  interference  is  from  the 
molecule  "^(BeOP)+.  System  response  to  this  molecule  varied  between  sessions 
and  may  be  attributable  to  aging  in  the  system's  electron  multiplier.  The 
annealed  profiles  of  high  fluence  implants  presented  here  were  obtained  when 
instrumental  conditions  were  such  that  this  molecular  interference  was  not 
observed.  Profiles  taken  during  sessions  when  interference  was  observed 
contained  ^Fe  profile  peaks  even  after  subtracting  out  the  molecular  inter¬ 
ference.  In  addition,  the  reality  of  the  structures  observed  in  the  ^^Fe 

profiles  of  high  fluence  implanted  material  is  confirmed  by  similar  structures 
52 

in  the  Cr  profiles  (which  have  no  apparent  molecular  interference) . 

28  56 

Another  possible  interference  at  mass  56  is  the  Si  dimer,  S^. 

Mass  28  profiles  taken  from  samples  that  were  Si^N^  encapsulated  and  annealed 

at  800°C  for  30  min  contain  a  strong  surface  signal  well  above  background, 

apparently  resulting  from  Si  in-diffusion  from  the  encapsulant.  Our  analysis 

56  4* 

of  dimer  yield  from  crystalline  Si  indicates  that  the  Si^  yield  is  ~  2.5 

28  + 

orders  of  magnitude  lover  than  the  Si  yield  (or,  ~  0.4%  of  the  parent  ion). 
56  4*  28 

The  yield  of  S^  ^rom  dilute  Si  in  InP  should  be  even  lower  since  dimer 
production  is  roughly  proportional  to  the  square  of  the  parent  concentration. 
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Taking  this  into  account,  interference  is  possible  only  at  the  InP 

surface  in  the  data  presented  here.  Additional  facts,  however,  suggest  that 

even  at  the  surface  this  possible  dimer  interference  may  not  be  significant. 

54 

Depth  profiles  of  the  Fe  isotope  show  iron  surface  accumulation  regions 

56  56 

similar  to  those  seen  in  Fe  profiles.  Estimates  of  3  Fe  signal  levels 

54  56 

based  on  the  isotopic  abundance  of  Fe  predict  Fe  signal  levels  comparable 

to  those  observed  experimentally.  Also,  for  the  profiles  exhibited  here, 

the  first  channel  of  ^Fe  data  was  usually  acquired  at  a  depth  where  ^Si* 
interference  is  expected  on  the  basis  of  our  observations  to  be  minimal. 

9 

7.2.1  SIMS  Profiles  of  Annealed  Be  Implants 

Figure  7.6  shows  the  SIMS  profiles  resulting  from  implantation  of 

13  2  9  o 

1  x  10  cm  Be  into  Fe-doped  InP  with  subsequent  30  min  anneals.  At  550  C 

9  56 

the  Be  atomic  profile  exhibits  little  diffusion.  The  Fe  profile,  however, 

has  undergone  some  movement  even  at  this  relatively  moderate  annealing  temp¬ 
erature.  It  is  observed  that  the  Fe  concentration  is  higher  than  the  original 
background  level  both  at  the  surface  and  in  an  area  to  the  left  (surface) 
side  of  the  implanted  peak.  As  will  be  discussed  later,  this  particular 
region  coincides  closely  to  the  peak  in  the  theoretical  damage  profile  created 
by  the  implantation  process.  These  accumulation  regions  are  followed  by  an 
^Fe  depletion  region  extending  to  ~  1  urn  before  returning  to  the  background 

level.  At  650°C  and  750°C  rapid  diffusion  of  the  initial  ^Be  profile  toward 

56 

the  surface  and  the  bulk  is  seen.  In  addition,  significant  movement  of  Fe 

is  observed,  with  the  previously  noted  depletion  zone  in  the  550°C  anneal 

deepening  and  broadening.  At  750°C  the  ^Fe  profile  reaches  a  minimum 

16  -3 

concentration  of  ~  1  x  10  cm  and  returns  to  the  background  doping  level 


at  a  depth  of  ~  2.0  -  2.5  pm. 
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Figure  7.6:  SIMS  unannealed  and  annealed  profiles  of  Fe-doped  InP  implanted 
"  13  -2  9 

with  100  keV,  10  cm  Be.  Thirty  minute  anneals  were  performed  with  a 
Si_N.  encapsulant  [138]. 
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9  14 

In  Figure  7.7  the  results  of  implanting  Be  to  a  fluence  of  1  x  10 

-2  9 

cm  and  subsequent  30  min  annealing  are  shown.  The  Be  profile  resulting 

from  a  550°C,  30  min  anneal  differs  markedly  from  that  seen  in  Figure  7.6(b). 

9 

Movement  of  Be  toward  the  surface  and  the  bulk  is  seen  in  this  case,  with 
a  shoulder  appearing  on  the  bulk  side  of  the  initial  profile.  An  unusual 

feature  of  this  profile  is  the  flat,  essentially  constant  concentration  of 

9  9 

Be  ranging  from  ~  1-2  pm  following  by  an  abrupt  drop  to  the  Be  background 

detection  level.  It  can  also  be  noted  that  ^Fe  getters  into  the  implanted 

region  slightly  above  its  background  concentration.  The  ^Fe  depletion 

zone  previously  mentioned  is  seen  to  return  rapidly  to  background  in  a  region 

9  9 

coinciding  with  the  abrupt  drop  in  the  Be  atomic  profile.  The  Be  profile 

resulting  from  a  650°C,  30  min  anneal  is  similar  to  the  550°C  anneal  in  that 

a  peak  is  still  observed  about  the  implanted  range,  followed  by  a  shoulder 

q 

on  the  bulk  side  of  the  profile.  However,  the  Be  profile  appears  to  return 

to  the  background  level  more  gradually  in  this  case  than  for  the  550°C  anneal. 

At  750°C  the  ^Be  profile  is  essentially  flat  and  featureless,  exhibiting  a 

very  gradual  drop  to  the  background  detection  limit.  The  Fe  profile  for 

this  anneal  consists  of  a  surface  accumulation  region  followed  by  an  extremely 

broad,  flat  depletion  zone.  Although  Figure  7.7(d)  does  not  show  ^Fe 

returning  to  background,  it  has  been  observed  in  other  profiles  taken  from 

this  sample.  This  return  takes  place  much  more  gradually  than  in  any  of  the 

other  profiles  from  this  set  however. 

Figure  7.8  shows  profiles  of  Fe-doped  InP  implanted  with  a  fluence 
15  -2  9 

of  1  x  10  cm  of  Be.  As  this  figure  illustrates,  annealing  this  material 
at  550°C  for  30  min  results  in  gross  restructuring  of  both  the  ^Be  and  ^Fe 
initial  pr^'lles.  Referring  to  Figure  7.8(b),  note  that  in  addition  to  the 
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Be  peak  in  the  implanted  region,  a  large  peak  occurs  on  the  surface  side 

of  the  implanted  peak  with  a  smaller  peak  evident  on  the  bulk  side.  The 

position  of  this  smaller  peak  correlates  well  with  the  observed  shoulders 

at  this  depth  in  the  other  implants.  The  most  dramatic  feature  of  this 

9 

particular  profile  is  the  constant-concentration  tail  of  Be  extending  several 

microns  into  the  bulk.  “*^Fe  redistribution  appears  to  be  correlated  with 
9 

the  motion  of  Be,  with  strong  gettering  occurring  on  the  surface  side  and 

9 

the  bulk  side  of  the  initial  Be  profile.  Deeper  profiles  indicate  that  the 

9 

constant-concentration  tail  of  Be  extends  to  a  depth  of  ~  4  pm  before 
dropping  abruptly  to  the  background  detection  limit  with  a  corresponding 
sudden  return  to  the  background  doping  level  in  the  ^Fe  profile  at  this 

same  depth.  At  650°C,  ^^Fe  is  observed  to  getter  into  the  region  to  the 

9  9 

right  of  the  Be  peak  and  into  the  implanted  peak  area.  The  Be  profile  at 

this  temperature  still  exhibits  a  peak  in  the  implanted  region,  a  shoulder 
on  the  bulk  side  of  the  profile,  and  a  constant-concentration  tail  extending 
from  1-2  ym  in  the  bulk.  The  ^Be  profile  resulting  from  a  30  min,  750°C 
anneal  exhibits  a  peak  with  a  concentration  reduced  by  approximately  one 
order  of  magnitude  from  the  peak  implant  concentration.  The  "*^Fe  profile  of 
this  anneal  exhibits  the  previously  observed  surface  accumulation  and  a  deep 
depletion  zone  extending  ~  2  ym  into  the  bulk  before  returning  to  background. 

An  unusual  feature  of  these  profiles  is  the  retrograde  motion  of  the  abrupt 
drop  in  ^Be  (or  of  the  rise  in  ^Fe)  as  the  annealing  temperature  is  increased. 
Higher  anneal  temperatures  result  in  shallower  regions  of  impurity  redistri¬ 
bution,  which  is  contrary  to  simple  diffusion  mechanisms. 

To  date  most  reports  on  implantation  of  InP  have  employed  annealing 
schedules  of  T  >  700°C  for  ~  10-15  min  [18,19,76,78].  Therefore  we  examined  the 
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9Be  Implant,  100  keV,  10 18  cm  "2 


Figure  7.9:  SIMS  unannealed  and  annealed  profiles  of  Fe-doped  InP  implanted 
15  -2  9 

with  100  keV,  10  cm  Be.  Fifteen  minute  anneals  were  performed  with  a 


Si.N.  encapsulant  [138]. 


9  56 

effect  of  15  min  anneals  on  the  restructuring  of  Be  and  Fe  profiles  in 


9  15  -2 

material  implanted  with  a  Be  fluence  of  10  cm 


From  Figure  7.9  it  can 


be  seen  that  these  profiles  are  qualitatively  similar  to  those  obtained  with 


30  min  anneals.  The  redistribution  of  ^Be  and  ^Fe  in  lO^  cm  ^  implants 


appears  to  be  rapid,  with  significant  structure  and  diffusion  occurring  within 
only  15  min  of  annealing. 


Figure  7.10  contains  impurity  profiles  after  30  min  anneals  of  Cr- 
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doped  InP  implanted  with  1  x  10  cm  fluence  of  Be.  The  resultant  Be 


profiles  are  similar  in  most  respects  to  correspondingly  processed  Fe-doped 


InP.  Since  molecular  interference  is  not  expected  or  seen  at  mass  52,  the 


52  56 

observed  Cr  peaks  in  Figure  7.10  suggest  that  the  Fe  peaks  in  Figures 


7.8  and  7.9  are  real  and  not  entirely  the  result  of  molecular  interference. 


Evidence  of  retrograde  motion  of  the  Be  redistribution  front  is  again  observed 


in  these  profiles  [note  the  different  depth  scale  in  7.10(b)].  Qualitatively, 


52  56 

Cr  appears  to  redistribute  much  as  Fe  but  with  several  quantitative 


differences.  Although  near-surface  SIMS  data  is  suspect,  it  certainly  appears 


52  56 

that  Cr  does  not  getter  to  the  surface  as  strongly  as  does  Fe.  In  fact, 


52  52 

the  near-surface  Cr  signals  are  comparable  to  the  original  Cr  background 


56  52 

level.  Another  significant  difference  between  Fe  and  Cr  appears  in  the 


56t 


minimum  concentration  observed  in  the  depletion  zones.  Whereas  the  Fe 


concentration  is  observed  to  deplete  approximately  an  order  of  magnitude  in 


52 


the  depletion  zones,  Cr  depletes  approximately  two  orders  of  magnitude  in 


56  52 

this  depletion  zone.  These  differences  between  Fe  are  Cr  and  consistent 


with  the  profiles  discussed  in  Section  7.1. 


In  order  to  determine  the  role  of  the  substrate  in  the  observed  Be 


profiles,  an  undoped,  12  um-thick  vapor  phase  epitaxial  (VPE)  layer  of  InP 
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9Be  Implant,  100  keV,  1015  cm'2 


Figure  7.10:  SIMS  unannealed  and  annealed  profiles  of  Cr-doped  InP  implanted 
15  -2  9 

with  100  keV,  10  cm  Be.  Thirty  minute  anneals  were  performed  with  a 
Si.N,  encapsulant  [138]. 


Density  (cm~3) 


9Be  Implant,  100  keV,  10l5cm-2 
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9 

was  implanted  with  Be.  The  resultant  profiles  are  shown  in  Figure  7.11. 

In  most  respects  these  profiles  are  similar  to  those  observed  in  comparable 
anneals  of  Fe-doped  or  Cr-doped  material.  Note  however  in  the  550°C  anneal 

9 

shown  here  the  lack  of  the  constant-concentration  Be  tail  observed  in 
Figures  7.8(b)  and  7.10(b). 

7.2.3  Discussion :  Surface  Phenomena 

As  previously  mentioned,  one  must  be  cautious  in  attaching  too 
much  quantitative  importance  to  the  initial  surface  signals  obtained  by 
SIMS.  Surface  oxides  will  enhance  the  yield  of  positive  secondary  ions  with¬ 
in  the  first  hundred  angstroms  [52,145].  Additionally,  the  sputtering  rate 
and  data  acquisition  time  can  smear  this  effect  out  over  a  larger  distance. 
Quantitatively,  the  first  channel  of  ^Fe  data  may  to  some  degree  reflect 
56Si2  interference  as  described  earlier.  However,  although  the  exact  surface 
concentrations  are  somewhat  in  doubt,  it  is  apparent  that  in  the  majority  of 

the  anneal  profiles  ^Fe  exhibits  strong  getterlng  to  the  surface  above  the 

9  * 

background  doping  concentration  and  that  Be  also  shows  diffusion  towards  * 

the  surface,  depending  on  fluence  and  annealing.  It  is  interesting  that  the 

52  56 

surface  gettering  of  Cr  is  not  as  strong  as  that  of  Fe.  The  concentra- 
52 

tions  of  Cr  in  the  annealed  profiles  shown  in  Figure  7.10  are  either  at  the 

52 

original  doping  level  or  slightly  below  this  level.  This  behavior  of  Cr 
and  ^^Fe  is  consistent  with  previously  observed  effects  in  annealed-only 
samples  and  material  implanted  with  "neutral"  impurities  (He,  ^B)  [21]. 

Various  mechanisms  have  been  proposed  to  explain  the  observation 
of  impurity  redistribution  at  surfaces,  including:  (i)  strain-induced 
redistribution  resulting  from  thermal  mismatch  between  the  encapsulant  and 
the  substrate  [147];  (ii)  surface  states  and  surface  impurities  giving  rise  to 
a  near  surface  electric  field  [163,164];  (iii)  surface  defects  resulting  from 
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mechanical  damage  or  annealing  [148] .  Recent  work  Indicates  that  Cr  re¬ 
distribution  in  GaAs  can  occur  without  an  encapsulating  layer  [58]  and  is  in 

fact  influenced  by  the  partial  pressure  of  As  during  the  annealing  of  bare 

52  56 

GaAs  surfaces  [144,146].  Differences  in  Cr  and  Fe  accumulation  may  pos¬ 
sibly  be  understood  in  terms  of  surface  solubility  effects  and  the  original 

doping  concentration  of  the  compensating  dopant  [58].  For  all  three  impurities, 
9  52  56 

Be,  Cr,  and  Fe,  diffusion  into  the  Si^N^  must  also  be  considered.  Previous 
work  done  in  our  group  using  SIMS  has  shown  that  Be  can  be  detected  in  the 
Si^N^  encapsulating  layer  after  annealing  Be-implanted  GaAs  at  900°C  for 
30  min  [165] . 

7.2.4  Discussion:  Near  Surface  Phenomena 

The  next  region  we  wish  to  consider  extends  from  ~  0.1  -  0.6  pm. 
Several  aspects  of  the  observed  redistribution  in  this  region  may  be  under¬ 
stood  on  the  basis  of  implantation-related  damage.  Shown  in  Figure  2.2  are 

the  LSS  Gaussian,  the  Gaussian  theoretical  damage  profile,  and  the  SIMS 

15  -2  9 

profile  of  the  unannealed  10  cm  Be  implant.  The  shape  of  the  SIMS 

profile  and  its  fit  to  four-moment  distributions  will  be  discussed  in  Chapter 

8.  The  peak  of  the  Gaussian  damage  profile  occurs  at  ~  2900  X  while  the 

9  o 

LSS  mean  range  for  100  keV  Be  is  ~  3640  A  (see  Table  2.1). 

13  1  ^  -2 

At  the  fluences  of  10  and  10  cm  ,  the  position  of  the  damage 
peak  at  2900  X  agrees  well  with  the  position  of  an  Fe  peak  just  to  the  left 
of  the  implant  peak  [see  for  example  Figures  7.6(b),  7.7(b),  and  7.7(c)]. 

In  the  previous  section  we  examined  the  ^Fe  profiles  resulting  from  the 
implantation  of  He  with  the  same  scanned  current,  mean  depth,  and  peak  con¬ 
centration  as  the  100  keV,  10^  cm  ^  ^Be  implants  discussed  here.  The  ^Fe 
profiles  described  there  also  showed  gettering  into  the  damage  region  which 
peaked  at  ~  3200  X  for  51  keV  He.  Those  Fe  profiles,  however,  resemble  the 
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10  cm  profiles  shown  in  Figure  7.6  more  than  the  Fe  profiles  of  Figures 

15  -2  9 

7.8  and  7.9  which  result  from  10  cm  Be  implants.  Although  the  peak 

15  -2  9 

concentration  of  He  was  equal  to  that  of  a  100  keV,  10  cm  Be  implant, 
the  lower  mass  of  He  results  in  much  less  energy  deposition  into  the  lattice. 

For  this  reason  the  Fe  profiles  observed  in  the  He  implants  of  Section  7.1.2 

9 

resemble  those  resulting  from  low  fluence  Be  implants.  Such  gettering  of 
impurities  to  implantation-induced  damage  areas  is  frequently  observed  in 
Si  and  in  GaAs  [149,166,167]. 

The  redistribution  of  Be  and  the  compensating  dopants  differ 
15  -2 

markedly  in  the  10  cm  implanted  samples  from  the  profiles  observed  in 
the  lower  fluence  implants.  As  can  be  seen  in  Figures  7.8(b),  7.9(b), 

7.10(b),  and  7.11,  twin  peak  structures  appear  in  the  Be,  Fe,  and  Cr  profiles. 

In  all  cases  the  drop  between  the  impurity  shoulders  occurs  in  the  approxi¬ 
mate  position  of  the  damage  peak  at  ~  2900  X.  Such  "up-hill"  impurity  flow 
cannot  be  explained  by  damage  enhanced  diffusion,  which  should  lead  to 
profile  flattening  rather  than  this  type  of  redistribution  [168].  There  are 
at  least  two  possible  models  for  this  behavior. 

Such  two-peaked  structures  have  sometimes  been  attributed  to  the 
presence  of  an  amorphous  layer  [58,169].  Electron  channelling  patterns  obtained 
from  implanted,  unannealed  samples  are  shown  in  Figure  7.12.  Comparing  the 
pattern  of  Figure  7.12(a)  to  the  other  patterns  of  Figure  7.12  indicates 
that  based  on  the  loss  of  higher  order  Coates  lines,  the  implanted  samples 
are  heavily  damaged  as  expected.  The  presence  of  the  major  Coates  lines  in 
these  samples  however,  indicates  that  the  disordered  material  is  not  amorphous. 
In  view  of  these  results,  it  appears  that  amorphization  is  not  the  cause  of 


the  observed  Be  redistribution. 


♦ 


Figure  7.13:  (a)  Electron  channeling  pattern  (ECP)  from  the  surface  of  an 

unimplanted  sample,  (b)  ECP  at  a  depth  of  ~  3000  X  below  the  surface  from 

15  -2  9 

a  sample  implanted  with  100  keV,  10  cm  Be  and  annealed  for  30  min  at 
550°C  [138]. 
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Another  explanation  of  the  double  peak  structures  can  be  found  in 
the  models  proposed  by  Baruch  [168]  and  Akutagawa  ejt  al.  [150],  for  similar 
structures  in  Si.  According  to  these  models.  Impurity  shoulders  with  a 
central  dip  can  result  from:  (1)  distinct  but  interacting  currents  of  defects 
and  impurities  with  their  own  gradients  and  diffusivities;  (2)  impurity- 
defect  pairs  which  diffuse  rapidly  away  from  the  peak  of  the  localized  damage 
distribution.  These  explanations  for  the  observed  dip  and  shoulders  are  both 
consistent  with  the  fact  that  this  structure  appears  at  the  position  of  the 
theoretical  damage  peak. 

In  addition  to  the  channelling  patterns  of  unannealed  material,  we 

o  15 

have  examined  the  ECP  from  550  C,  30  min  anneals  of  10  cm  Be  implanted 
InP.  Channelling  patterns  at  depths  of  ~  0,  1500,  3000,  and  4500  X  have 
been  analyzed  and  are  virtually  indistinguishable  from  unimplanted  material. 
Only  the  pattern  shown  in  Figure  7.13  at  a  depth  of  ~  3000  8  (near  the  theo¬ 
retical  damage  peak)  showed  a  barely  perceptible  loss  of  higher  order  Coates 
lines.  It  is  obvious  from  Figure  7.12  and  7.13  that  significant  reordering 
of  the  lattice  occurs  at  this  low  temperature,  although  electrical  activation 
of  the  implanted  impurity  generally  requires  higher  annealing  temperatures. 
7.2.5  Discussion:  Bulk  Phenomena 

Several  unusual  profile  features  are  evident  in  the  bulk  regions 

of  the  samples  analyzed  in  this  work.  These  are:  (1)  constant-concentration 
9 

tails  of  Be  extending  over  distances  on  the  order  of  microns;  (2)  constant- 

52  56 

concentration  depletion  zones  of  the  compensating  impurity  (  Cr  or  Fe) 
extending  over  similar  distances;  (3)  correlation  between  the  diffusion 
profile  of  the  Be  and  the  compensating  impurity  (i.e.,  the  depth  at  which 
the  compensating  impurity  returns  to  its  doping  level  coincides  with  the 


depth  at  which  the  constant-concentration  tail  of  Be  drops  to  the  detection 

9 

limit);  and  (4)  retrograde  motion  of  the  Be  diffusion  front  (or  of  the 
compensating  impurity  depletion  zone  front)  with  higher  annealing  tempera¬ 
tures.  Of  these  features,  depletion  of  the  compensating  impurity  on  the  bulk 
side  of  the  implanted  peak  has  previously  been  observed  in  annealed  InP 
implanted  with"neutral"  elements  (He,  ^B)  as  described  in  Section  7.1.2. 

To  our  knowledge,  the  other  features  reported  here  have  not  been  previously 
observed  in  conjunction  with  implantation  of  InP  or  GaAs. 

Perhaps  the  most  striking  features  of  the  impurity  redistribution 

9 

are  the  long,  essentially  flat  Be  tails  and  the  accompanying  Fe  or  Cr 

depletion  zones.  Such  diffusion  cannot  result  from  normal  Fickian  (random 

walk)  diffusion.  Similarly,  conventional  field-enhanced  diffusion  [163,164,170] 

does  not  explain  the  correlation  between  the  ^Be  and  "^Cr  or  ^Fe  profiles. 

To  explain  such  correlations  with  field  enhanced  diffusion  would  require  that 
9  52  56 

Be,  Cr,  and  Fe  all  have  the  same  effective  diffusion  coefficients. 

When  considering  these  bulk  phenomena  the  following  relevant  factors 

emerge  from  this  work:  (1)  the  dose  dependent  location  of  the  end  of  the 

^Be  diffusion  tail  (or  of  the  ^Cr  or  depletion  zone)  at  550°C;  (2)  the 

retrograde  motion  of  these  features  with  elevated  annealing  temperature; 

(3)  the  more  gradual  sloping  of  the  ^Be  diffusion  front  at  750°C  versus  the 

constant  concentration  tails  seen  at  550°C;  and  (4)  the  lack  of  constant- 

9  15  -2 

concentration  tails  of  Be  in  10  cm  implanted  and  annealed  VPE  material. 

9 

The  fact  that  the  constant-concentration  Be  tail  does  not  appear  in  VPE 

InP  indicates  that  either  the  compensating  dopants  (Fe  or  Cr)  or  other  residual 


impurities  and/or  defects  in  semi-insulating,  bulk  InP  are  necessary  in  the 

9  52  56 

creation  of  this  tail.  Formation  of  pair  complexes  of  Be  with  Cr  or  Fe 


125 


might  explain  the  observed  correlation  between  the  Be  and  the  compensating 

9 

impurity  profiles.  It  is  unclear  if  energy  considerations  prevent  Be  from 

+3 

complexing  with  the  deep  level  acceptors  Fe  (neutral  state  Fejn)  [171,172]  or 

+2  56  52 

Cr  (neutral  state  Cr  )  [171]*  Based  on  the  depletion  of  Fe  or  Cr  in  the 

9  9 

regions  of  the  Be  tail,  it  seems  more  likely  that  Be  participates  with  the 
compensating  impurities  in  a  cooperative  diffusion  mechanism.  Such  possible 

9 

effects  alone,  however,  cannot  explain  the  retrograde  behavior  of  the  Be 

52  56 

tail  and  the  accompanying  end  of  the  Cr  or  Fe  depletion  zone.  Based  on 

the  temperature  dependence  of  this  phenomenon,  a  defect  or  defect  complex 

with  enhanced  diffusion  could  be  responsible  for  the  deep  diffusion  at  low 

temperatures  (i.e.,  550°C).  The  concentration  of  such  defects  or  defect 

complexes  is  reduced  at  higher  temperatures  (i.e.,  650-750°C),  resulting  in 

shallower  redistribution  of  the  impurities.  Annealing  effects  of  such  a 

defect  structure  with  enhanced  diffusion  can  explain  the  more  gradual  diffu- 
9 

sion  front  of  Be  observed  at  750°C  versus  the  abrupt  fronts  seen  at  550°C. 

Additionally,  by  assuming  these  defect  structures  to  be  implant  damage  related, 

9 

it  is  possible  to  account  for  the  observed  dose  dependent  depth  of  the  Be 

diffusion  front  seen  in  the  550°C  anneals.  However,  it  is  necessary  to  connect 

52  56 

such  defect  structures  with  the  presence  of  Cr  or  Fe  or  other  impurities 

9 

in  bulk  InP,  without  which  the  flat  tails  of  Be  are  not  formed.  Alternately 

9 

the  constant-concentration  Be  may  result  from  a  low  Be  solubility  in  undam¬ 
aged,  semi-insulating  InP  (published  results  for  Be-doped  LPE  layers[173] 

18  —3 

indicate  no  saturation  of  carrier  concentration  up  to  '  6  x  10  cm  )  with 
the  rate-limiting  step  of  the  process  being  the  release  of  Be  from  the 
damaged  layer.  If  it  is  further  assumed  that  the  Be  solubility  increases 


with  temperature  it  is  possible  to  understand  the  shallower  Be  profiles  at 
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TABLE  7.1 


Pauling  covalent  tetrahedral  radius[176],  r  »  Phillips  covalent  tetrahedral 

r  a 

radius [17 7] ,  rpH,  Pauling  electronegativity  [176]  e  and  Phillips  electro¬ 


negativity  [177] ,  e 


PH 


Ion 


PA 

d) 


PH 

(X) 


'PA 
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higher  temperatures.  Solubility  effects,  however,  do  not  explain  the  corre- 

9  52  56 

lation  of  the  Be  profile  to  that  of  the  Cr  or  Fe  profile  in  semi- 

insulating  InP,  so  that  this  solubility  model  also  seems  to  require  cooper- 

9 

ative  diffusion  between  Be  and  the  compensating  impurities. 

It  has  been  observed  that  an  ideal  substitutional  impurity  should 
possess:  (1)  a  covalent  tetrahedral  radius  close  to  the  host  atom  it  replaces; 

(2)  small  electronegativity  difference  with  the  host;  and,  (3)  valence  close 
to  the  host[174].  A  recent  review  of  redistribution  effects  in  III-V  compounds 
has  stressed  the  importance  of  electronegativity  and  size  effects  in  determin¬ 
ing  diffusion  behavior [ 1 75 ] .  We  have  tabulated  the  values  of  electronegativity 
and  the  covalent  tetrahedral  radii  for  relevant  atoms  as  calculated  by 

Pauling [176]  and  Phillips [177]  in  Table  7.1.  From  this  table  it  is  clear  that 

10  9 

B  is  farther  from  In  in  both  size  and  electronegativity  than  is  Be.  As 

described  in  Section  7.1.2,  ^Be  implants  of  126  keV  and  6.79  x  10*^  cm  ^ 

dose  were  performed  to  simulate  the  damage  resulting  from  a  100  keV,  10^ 

-2  9  20  -3 

cm  Be  implant  (energy  deposition  of  2  x  10  keV  cm  ).  The  resulting 

^B  atomic  profiles  show  very  little  diffusion  after  750°C,  30  min  anneals, 

g 

in  contrast  to  identical  anneals  of  any  Be  dose  considered  here  [178], 
Electronegativity  and  size  considerations  alone  are  therefore  not  sufficient 

9 

to  explain  the  rapid  diffusion  of  Be. 

7.3  Redistribution  of  Implanted  Mg  in  SI  InP 

The  need  for  a  well  behaved  (i.e.,  high  activation,  low  diffusivity) 
implanted  acceptor  in  InP  has  already  been  described  in  the  preceding  section. 

The  examination  of  high  dose  Mg  implants  in  InP  resulted  from  the  findings 

of  the  poor  annealing  characteristics  for  high  dose  Be  implants  in  InP.  In 

terms  of  atomic  mass,  Mg  is  more  appropriate  for  implantation  than  other 

acceptors  of  heavier  mass  such  as  Zn  or  Cd.  The  lower  mass  of  Mg  should 
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result  in  less  damage  and  afford  greater  variation  in  the  location  of  the 

projected  range  than  with  heavier  ions.  From  the  Phillips  values  for 

covalent  tetrahedral  radii  and  electronegativity  given  in  Table  7.1,  Mg 

should  be  a  better  substitutional  impurity  for  In  than  is  Be.  To  study 
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implanted  Mg,  the  isochronal  annealing  characteristics  of  250  keV,  10  cm 
26 

implants  of  Mg  have  been  examined  in  this  work  using  SIMS.  This  isotope 

of  Mg  was  chosen  because  it  readily  provides  the  low  ion  current  required  in 

this  study  to  prevent  beam  heating  of  the  sample.  Implantation  took  place 

_2 

at  room  temperature  using  a  Mg  current  density  of  ~  0.1  yA  cm  (see  Table  2.1). 

Shown  in  Figure  2.3  are  the  SIMS  atomic  profiles  of  the  as-implanted 

26 

Mg  and  the  theoretical  profiles  for  the  impurity  and  the  accompanying 

damage.  To  facilitate  comparisons  with  our  Be  implant  work,  the  Mg  peak 
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concentration  is  roughly  comparable  to  that  of  Be  in  the  10  cm  implant 

discussed  in  the  preceding  section.  The  mode  of  the  LSS  Gaussian  for  Mg 

occurs  at  2994  X  and  the  mode  of  the  Gaussian  damage  profiles  is  located  at 

1630  X  (see  Table  2.1).  The  fit  of  four  moment  distributions  to  the 

experimental  SIMS  profile  will  be  discussed  in  Chapter  8. 

Shown  in  Figure  7.14  are  the  SIMS  profiles  for  30  min  anneals  from 

550-750°C.  Redistribution  of  the  implanted  Mg  is  seen  to  occur  at  550°C. 

19  -3 

This  anneal  results  in  a  flat  profile  of  *  1.3  x  10  cm  of  Mg  extending 

from  the  surface  to  a  depth  of  approximately  4000  X.  At  higher  annealing 

temperatures  this  region  is  increasingly  depleted  of  Mg,  resulting  in  a 

16  -3 

minimum  concentration  of  -  3.7  x  10  cm  at  750°C.  Comparing  Figures  2.3 

and  7.14(d),  we  observe  that  the  Mg  depletion  zone  occurs  in  the  region  of 

maximum  implant  damage.  Along  with  this  depletion  zone,  an  accumulation 
region  of  Mg  occurs  on  the  bulk  side  of  the  initial  Mg  profile.  As  Mg  is 
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depleted  from  the  near  surface  region  with  increased  annealing,  the  width 
of  the  Mg  accumulation  region  increases.  This  suggests  that  Mg  is  gettering 
into  the  region  just  below  the  maximum  implant  damage.  The  compensating  Fe 
dopant  also  getters  into  the  same  region  as  the  annealing  proceeds.  This 
Fe  accumulation  zone  increases  in  width  with  increasing  annealing  tempera¬ 
ture,  essentially  matching  the  Mg  profile  shape.  This  redistribution  of 
Fe  to  closely  reflect  the  profile  of  the  implanted  acceptor  was  also  observed 
in  our  Be  implant  studies. 

The  observed  depletion  and  accumulation  zones  of  Mg  appear  to  be 

related  to  the  presence  of  an  amorphous  surface  layer  and  an  accompanying 

damaged  crystalline  region  below  the  amorphous  layers.  For  a  Mg  dose  of 
15  -2 

1  x  10  cm  implanted  at  250  keV,  the  energy  deposition  density  is  calcu- 

21  -3 

lated  to  be  "  1.7  x  10  keV  cm  (approximated  from  the  product  of  the 
dose  and  the  nuclear  energy  loss  per  unit  length) .  Assuming  the  validity 
of  the  energy  deposition  density  model  [23,179],  this  is  from  3-6  times  the 
critical  energy  deposition  density  necessary  for  the  amorphization  of  InP, 
as  calculated  for  the  cases  of  Ar  and  Si  at  room  temperature  [180,181].  Electron 
channelling  patterns  (ECP)  for  as-implanted  InP  were  obtained  from  the  surface 
and  from  a  depth  of  ~  1600  X.  Both  surfaces  were  devoid  of  Coates  lines, 
indicating  an  amorphous  implanted  region.  The  material  was  re-examined  a 
year  after  implantation  and  was  still  found  to  be  amorphous,  supporting  a 
similar  observation  by  Kennedy[i81] .  This  amorphous  surface  layer  is  located 
in  the  same  region  where  Mg  depletion  zones  are  formed  during  annealing, 
suggesting  a  causal  relationship.  Figure  2.3  illustrates  that  in  the  region 
of  gettering  (5500  -  6500  X)  the  energy  deposition  is  approximately  two 
orders  of  magnitude  less  than  at  the  peak  of  the  damage  profile.  Presumably 
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Figure  7.15:  (a)  Electron  channelling  pattern  (ECP)  from  an  unimplanted 

InP  surface,  (b-d)  ECP  at  depths  of:  (b)  surface;  (c)  ~  1600  X  below  the 
surface;  and  (d)  ~  3900  X  below  the  surface.  These  samples  were  implanted 
with  250  keV,  1015  cm  2  Mg  and  annealed  for  30  min  at  550°C  [188]. 
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this  region,  although  damaged,  is  not  driven  amorphous.  Impurity  accumula¬ 
tion  into  the  damaged  crystalline  region  below  an  amorphous  implanted  layer 
has  been  previously  observed  in  Si  and  attributed  to  gettering  by  residual 
defects [182].  A  similar  mechanism  apparently  causes  Mg  and  Fe  gettering  into 
the  damaged  crystalline  region  in  InP. 

Although  an  extensive  study  of  channelling  patterns  from  these 
samples  has  not  been  done,  some  channelling  patterns  from  annealed  material 
have  been  obtained.  Figure  7.15  exhibits  channelling  patterns  from  material 
annealed  for  30  min  at  550°C.  The  lack  of  contrast  and  of  any  Coates  lines 
suggest  than  an  amorphous  layer  still  exists  in  the  near  surface  region 
following  this  anneal.  The  results  from  material  annealed  for  30  min  at 
650 °C  are  similar  to  those  shown  in  Figure  7.15.  ECP  obtained  from  30  min 
and  60  min  anneals  at  750°C  are  shown  in  Figures  7.16(b)  and  7.16(c).  In 
these  photographs  the  backscattered  electron  signal  was  not  differentiated 
as  in  Figure  7.15  because  the  weak  bands  visible  in  7.16(b)  and  7.16(c)  were 
too  low  in  contrast  to  yield  a  visible  pattern  in  the  differentiated  signal 
mode.  The  surfaces  of  both  the  30  and  60  min,  750°C  annealed  material  revealed 
no  observable  channelling  patterns.  The  faint  channelling  patterns  of  7.16(b) 
and  7.16(c)  were  obtained  at  depths  near  the  mode  of  the  Gaussian  damage 
profile.  Although  some  recrystallization  has  apparently  occurred  at  750°C, 
the  patterns  are  of  very  poor  quality  when  compared  to  those  from  unimplanted 
material.  These  results  seem  to  be  inconsistent  with  other  published 
results [180,181]  which  indicate  that  major  reordering  of  amorphous  InP  occurs 
in  the  range  from  450-500°C.  It  has  been  reported  that,  in  addition  to  the 
degree  of  crystallinity,  factors  such  as  lattice  strain  of  a  few  percent 
can  result  in  loss  of  resolution  or  contrast  in  electron  channelling  patterns 


Figure  7.16:  (a)  Electron  channelling  pattern  (ECP)  from  an  unimplanted  InP 

surface,  (b)  ECP  at  a  depth  of  ~  1600  %  below  the  surface  of  Mg  implanted 
sample  annealed  for  30  min  at  750°C.  (c)  ECP  at  a  depth  of  ~  1900  X  below 

the  surface  of  Mg  implanted  sample  annealed  for  60  min  at  750°C.  The 
backscattered  electron  signal  has  not  been  differentiated  as  was  the  case 
in  Figure  7.15  [188]. 


[183].  For  the  Mg  implants  discussed  here,  lattice  strain  may  be  significant 

and  could  therefore  be  responsible  for  the  discrepancy  between  our  results  and 

those  of  others  regarding  the  recrystallization  temperature  for  amorphous  InP. 

The  annealing  characteristics  of  Mg-amorphized  InP  shown  in  Figure 

7.14  may  explain  the  electrical  characteristics  of  Mg-implanted  InP  reported 

by  Donnelly  and  Hurwitz  [184].  For  room  temperature  implants  of  150  keV  of  Mg 
14  -2 

fluence  above  10  cm  they  find  a  drastic  reduction  in  sheet  hole  concen- 

14  -2 

tration.  Our  calculations  indicate  that  a  10  cm  dose  of  150  keV  Mg  yields 
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a  deposited  energy  density  of  "  2  x  10  keV  cm  ,  which  is  roughly  the  thresh¬ 
old  for  amorphization  of  InP,  assuming  the  validity  of  this  amorphization 

model.  Therefore  the  reduction  in  sheet  hole  concentration  for  fluences 
14  -2 

above  10  cm  could  be  the  result  of  amorphization  and  depletion  from  the 

amorphous  layer  as  observed  here.  This  assumption  is  consistent  with  the 

15  -2 

fact  that  sheet  concentration  increases  with  dose  up  to  10  cm  for  implants 
performed  at  150°C['184],  because  this  elevated  target  temperature  is  report¬ 
edly  sufficient  to  prevent  amorphization  of  InP [180,181], 

Figure  7.17  is  a  deeper  profile  taken  from  a  sample  annealed  at  750°C. 
A  flat,  constant-concentration  tail  of  Mg  and  an  accompanying  depletion  of 

Fe  extend  over  a  considerable  distance  in  the  bulk,  similar  to  that  observed 

in  the  Be  implant  work.  Because  the  Be  implants  did  not  amorphize  the  InP, 

it  can  be  assumed  that  this  feature  is  not  related  to  amorphization  effects. 

Possible  explanations  for  this  flat  tail  feature  have  been  discussed  in  the 

preceding  section.  As  was  the  case  for  Be,  the  correlation  between  the 

flat  Mg  tail  and  the  Fe  depletion  zone  is  striking. 

28 

7.4  SIMS  Profiles  of  Si  Implants  in  SI  InP 

Because  of  the  strongly  anomalous  redistribution  of  both  Be  and  Mg, 

Si  implants  were  performed  to  determine  if  this  donor  behaves  similarly. 


138 


SI  implantation  in  InP  is  of  considerable  interest  because  of  its  use  in 
the  formation  of  n-channels  in  InP  field  effect  transistors  (FET)  [126].  It 
has  generally  been  observed  that  Si  implants  result  in  good  donor  activa¬ 
tion  [19,  76,  98].  Based  on  values  for  the  covalent  tetrahedral  radius  and 
electronegativity  given  in  Table  7.1,  Si  appears  to  be  a  reasonably  suitable 
substitutional  impurity  for  the  In  lattice  atom. 

Using  SIMS  we  have  analyzed  the  annealing  behavior  of  240  keV, 

14  —2  28 

5  x  10  cm  Si  implanted  into  Fe-doped  InP.  Figure  2.3  contains  the 
28 

as-implanted  Si  atomic  profile  together  with  the  LSS  Gaussian  approximation 
and  the  theoretical  damage  profile.  The  mode  for  the  LSS  Gaussian  occurs 
at  2464  X  while  the  mode  of  the  energy  deposition  profile  occurs  at  -  1562  X 
(see  Table  2.1). 

28 

Shown  in  Figure  7.18  are  the  Si  profiles  resulting  from  imp lan- 

o  28 

tation  and  annealing  for  30  min  at  750  C.  Unfortunately,  the  Si  profiles 
suffer  from  a  rather  large  background  signal  (corresponding  to  ~  10^  cm~^) , 

due  to  molecular  interference.  Mass  spectra  from  the  surface  and  bulk  of 
uniplanted  samples  yield  a  strong  peak  at  mass  27  and  a  weaker  peak  at 

mass  28.  This  suggests  the  interference  observed  at  mass  28  is  dominated 

27 

by  the  hydride  of  the  mass  27  background  impurity  (presumably  Al)  rather 
28 

than  by  residual  Si  content.  Recent  studies  of  liquid  encapsulated 
Czochralski  grown  InP  using  spark  source  mass  spectrometry  reveal,  however, 
that  the  Si  content  of  such  single  crystals  is  greater  than  the  residual  Al 
concentration  (-  0.1  ppma  vs.  -  0.03  ppma)  [185,  186].  It  appears  most  likely 
that  the  mass  27  background  results  from  a  "memory  effect"  in  the  SIMS 
system's  immersion  lens  caused  by  the  sputtering  of  an  Al  grid  used  in 
aligning  the  system  optics  prior  to  sample  profiling. 


28Si  Implant,  240  keV,  5xl014  cm*2 


Unannealed  28Si 
Unannealed  56Fe 
750°C,30min>  28Si 
750°C,30min,  56Fe 


0.4  0.6  0.8 

Depth  (/im) 
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Figure  7.18;  SIMS  unannealed  and  annealed  profiles  of  Fe-doped  InP  implanted 

14  -2 

with  240  keV,  5  x  10  cm  Si.  30  min  anneals  were  performed  with  a 


Si^N^  encapsulant  [188], 
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In  spite  of  the  background  interference,  the  results  in  Figure  7.18 

28 

indicate  that  little  atomic  redistribution  of  the  implanted  Si  occurs  for 

o  28 

750  C,  30  min  anneals.  The  lack  of  major  redistribution  for  Si  observed 

here  appears  consistent  with  Hall  profiling  measurements  made  by  Donnelly 

and  Hurwitz  on  Fe-doped  InP  implanted  at  room  temperature  with  10  cm  of 

Si  at  400  keV  [19].  Significant  movement  of  Fe  is  detected,  however.  An  Fe 

accumulation  zone  followed  by  a  depletion  zone  is  formed  on  the  bulk  side 
28 

of  the  implanted  Si  peak.  With  the  exception  of  the  large,  near  surface 
depletion  zone  seen  in  Figure  7.14,  this  Fe  behavior  parallels  that  observed 
in  the  Mg  implant  series. 

28 

The  lack  of  major  Si  redistribution  in  these  anneals  is  signif¬ 
icant  when  compared  to  the  annealing  characteristics  of  high  dose  Mg  or  Be 

Implants.  Kennedy  has  found  via  channelling  measurements  that  a  fluence  of 
14  -2  28 

10  cm  of  200  keV  Si  is  sufficient  to  produce  an  amorphous  InP  layer 

at  room  temperature  [181].  Based  on  this  measurement  it  is  assumed  that  the 
14  -2  28 

5  x  10  cm  ,  240  keV  implant  of  Si  has  amorphized  the  InP  examined  here. 

21 

In  fact,  the  energy  deposition  density  calculated  for  this  implant  (1.1  x  10 

_3  26  21  -3 

keV  cm  )  is  comparable  to  that  of  the  Mg  implant  (1.7  x  10  keV  cm  ). 

Despite  this  similarity,  there  are  significant  differences  between  the  anneal¬ 
ing  characteristics  of  Si  and  Mg. 

The  implanted  Si  examined  here  results  in  far  less  redistribution 

9 

than  is  seen  in  lower  fluence  Be  implants  with  much  less  implant  related 
damage.  The  redistribution  behavior  of  the  implanted  acceptors  Mg  and  Be 
seems  to  parallel  the  rapid  diffusion  exhibited  by  Zn  and  Cd  during 
growth  [155-157].  It  has  previously  been  observed  that  donor  implants  result 
in  higher  electrical  activations  than  acceptor  implants  in  InP  [184].  It  seems 
that  the  redistribution  behavior  of  acceptors  in  InP  may  be  a  fundamental 
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limitation  leading  to  low  activations. 

It  is  clear  from  these  results  that  more  information  is  needed  on 
the  amorphization  and  the  kinetics  of  recrystallization  of  amorphous  InP. 
This  is  particularly  true  since  the  measured  average  displacement  energy 
of  the  lattice  atoms  in  InP  is  less  than  that  of  GaAs  or  Si  (7.7  eV  compared 
with  9.4  and  15.8  eV  for  GaAS  and  Si,  respectively)  [187].  As  a  result, 
amorphization  effects  may  be  encountered  more  frequently  during  normal 
processing  in  InP  than  has  been  the  case  with  GaAs  or  Si. 

7.5  Summary 

The  main  results  of  the  SIMS  studies  may  be  summarized  as  follows: 

(1)  High  temperature  anneals  with  a  Si^N^  cap  result  in  the  formation  of 
Fe  accumulation  regions  within  the  first  1000  &  of  the  surface  of 
Fe-doped,  SI  InP.  Comparable  anneals  of  Cr-doped  SI  InP  indicate  that 
Cr  does  not  exhibit  such  strong  surface  gettering. 

(2)  In  Si^N^  capped  SI  InP,  depletion  zones  of  the  compensating  impurity 
are  formed  following  high  temperature  anneals.  Cr  exhibits  larger 
concentration  gradients  in  these  zones  than  does  Fe. 

(3)  In  samples  implanted  with  electrically  and  chemically  inactive  impuri- 
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ties  (  He,  B),  Fe  and  Cr  exhibit  accumulation  into  the  implant 
damage  region  on  the  surface  side  of  the  projected  range  peak.  This 
occurs  in  addition  to  the  features  noted  in  (1)  and  (2)  which  are 
associated  with  the  use  of  Si^N^  cap. 

(4)  Under  typical  annealing  conditions  for  InP  (T  _>  700°C,  t  *  15-30  min), 

9 

implanted  Be  is  a  rapid  diffusant  even  for  fluences  as  low  as 

13  -2 

1  x  10  cm  . 


(5)  Short  (15  min)  anneals  of  high  fluence  Be  Implants  result  in  sub¬ 
stantial  diffusion  and  redistribution,  comparable  to  that  for  30 
min  anneals. 

(6)  Electronegativity  and  size  considerations  alone  are  not  sufficient 

9 

to  explain  the  rapid  diffusion  of  Be  in  InP. 

(7)  Two  peaked  structures  of  Be,  Cr,  and  Fe  appear  in  the  550°C  anneals 
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of  10  cm  implanted  material.  Such  structures  may  result  from 
rapid  defect  diffusion  away  from  the  damaged  region.  Electron  channel¬ 
ling  patterns  and  comparisons  with  *^B  implants  appear  to  rule  out 

an  amorphous  layer  as  the  cause  of  such  twin  peaked  structures. 

9 

(8)  Correlations  are  observed  between  Be  and  compensating  impurity 
profiles  in  annealed  high  fluence  implanted  samples.  The  compensating 
impurities  return  to  the  background  doping  level  in  the  same  region 

9 

that  the  Be  drops  to  background  detection  limits. 

9 

(9)  Flat,  constant-concentration  tails  of  Be  extending  over  distances 

on  the  order  of  microns  are  observed  in  semi-insulating  InP.  This 

anomalous  diffusion  is  not  observed  in  undoped  VPE  InP.  It  is  possible 

that  either  the  compensating  dopant  (Fe  or  Cr)  or  other  residual 

impurities  and/or  defects  in  semi-insulating  InP  play  a  role  in  the 

formation  of  such  tails.  Alternatively,  these  constant-concentration 
9 

Be  tails  may  result  from  a  low  Be  solubility  in  undamaged,  semi- 

insulating  InP  with  the  rate-limiting  step  of  the  process  being  the 

release  of  Be  from  the  damaged  layer. 

15  -2 

(10)  Implantation  of  10  cm  ,  250  keV  Mg  at  room  temperature  results  in 
amorphization  of  InP.  Electron  channelling  patterns  from  samples 
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annealed  for  30  min  at550°C  and  650°C  exhibit  no  patterns  at  the 
surface  or  at  depths  near  the  region  of  maximum  implant  induced  damage. 
Weak  channelling  patterns  are  observed  at  depths  near  the  damage  peak 
in  material  annealed  at  750°C  for  30  and  60  min. 

(11)  Drastic  redistribution  of  Mg  and  Fe  occurs  during  30  min  anneals  of 
such  implanted  material.  Depletion  zones  are  formed  in  the  region  of 
maximum  implant  damage,  followed  by  accumulation  zones  on  the  bulk 
side  of  the  initial  Mg  profile.  The  widths  of  the  accumulation  zones 
increase  with  temperature,  suggesting  gettering  of  impurities  from  the 
surface  into  this  bulk  region.  As  in  the  Be  implant  studies,  the  motion 
of  Fe  observed  here  closely  reflects  the  profile  of  the  implanted  Mg 
acceptor. 

(12)  Flat,  constant-concentration  tails  of  Mg  extending  over  distances  of 
1  um  or  more  are  observed  in  the  bulk  after  annealing.  The  depth  at 
which  Fe  returns  to  its  background  level  and  its  concentration  gradient 
at  this  point  are  closely  correlated  to  the  Mg  profile.  These  features 
have  been  observed  also  in  the  Be  study.  Since  the  Be  implants  did  not 
amorphize  InP,  it  may  be  assumed  that  these  features  are  not  related 

to  amorphization  effects. 

28 

(13)  Little  redistribution  of  Si  occurs  for  30  min,  750°C  anneals  of  240 
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keV,  5  x  10  cm  Si  implants.  Based  upon  calculations  and  avail¬ 
able  experimental  data,  this  implant  schedule  should  result  in  amorphous 

28 

InP.  This  behavior  of  Si  differs  significantly  from  that  of  Mg  in 

amorphous  InP.  Redistribution  of  Fe  does  occur  in  these  Si  implanted 

samples,  however,  resulting  in  an  accumulation  region  followed  by  a 

28 

depletion  zone  on  the  bulk  side  of  the  implanted  Si  peak. 
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8.  THEORETICAL  AND  EMPIRICAL  DISTRIBUTIONS  FOR  InP  ION  IMPLANTATION  PROFILES 


8. 1  Introduction 

In  the  course  of  the  Implantation  studies  in  InP  the  shortcomings 
of  the  LSS  Gaussian  approximation  became  apparent  for  light  mass  projectile 
implantations  dominated  by  electronic  stopping.  Although  the  theroetical 
LSS  projected  ranges  agreed  adequately  with  experiment,  profile  shapes  were 
encountered  which  differed  significantly  from  the  LSS  Gaussian.  There  has 
been  substantial  interest  [189-191]  in  finding  analytical  distributions  which 
when  combined  with  the  appropriate  parameters  will  predict  the  shape  of  ion 
implanted  profiles  with  greater  sophistication  than  the  common  Gaussian 
approximation  using  the  penetration  theory  of  Lindhard,  Scharff,  and  Schiott 
(LSS)  [31].  Distributions  such  as  the  joined  half-Gaussian  [192],  the  Edge- 
worth  [37],  the  Gram-Charlier  [35],  and  the  Pearson  system  [35,190,191]  have 
been  proposed  to  fit  asymmetrical  implant  profiles.  The  Pearson  IV  distribu¬ 
tion  in  particular  has  been  shown  capable  of  excellent  fits  to  experimental 
implant  profiles  in  Si  [190,191]. 

In  this  chapter,  four  implantation  profiles  involving  ^Be,  ^B,  ^Mg 

28 

and  Si  are  considered  in  the  intermediate  energy  range  3  <  e  <  7  (where  e 
is  the  LSS  dimensionless  energy  parameter  [31,37]),  The  purpose  of  this  chapter 
is  to  examine  the  fit  to  experimental  profiles  from  distribution  functions 
which  employ  four  moments  rather  than  only  two  moments  as  with  the  Gaussian. 

8. 2  Experimental  and  Gaussian  Profiles 

Qualitatively  the  experimental  (SIMS)  profile  shapes  shown  in 
Figures  8.1  -  8.4  can  be  understood  on  the  basis  of  electronic  stopping 
and  the  mass  ratio  (u)  of  the  average  target  mass  to  the  projectile  mass. 
Electronic  stopping  begins  to  dominate  over  nuclear  stopping  [24]  in  the  region 


—  Gaussian 
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Figure  8.1;  SIMS  profile  of  ^Be  implant  (100  keV,  1 
InP  and  theoretical  profiles  [195]. 


Figure  8.2:  SIMS  profiles  of  implant  (126  keV,  6.79  x  lO1^  cm  in  (100) 
InP  and  theoretical  profiles  [ 195 J . 
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e  *  3.  This  fact,  the  values  of  e,  and  the  fractional  energy  loss  to  elec¬ 
tronic  stopping  given  in  Table  8.1  indicate  that  electronic  stopping  dominates 

over  nuclear  stopping  in  all  the  implantations  discussed  here,  except  for 
28 

possibly  the  Si  implant.  Winterbon  has  noted  that  implants  dominated  by 

nuclear  stopping  result  in  profiles  containing  a  tail  extending  to  large 

depths,  whereas  in  implants  dominated  by  electronic  stopping  this  tail  is 

cut  off  and  the  profile  is  skewed  toward  the  surface  [35].  The  sign  and 

magnitude  of  the  third  moment  (skewness)  y^,  shown  in  Table  8.2  therefore 

indicates  the  role  of  electronic  stopping  in  these  implants  and  correctly 

predicts  the  skewing  of  the  profiles  toward  the  surface.  Another  factor 

9  10 

relevant  in  the  case  of  low  mass  projectiles  (high  u)  such  as  Be  and  B  is 

backscattering  from  target  atoms,  which  should  skew  the  distributions  even 

more  toward  the  surface.  These  observations  are  borne  out  by  our  experimental 

9  10 

profiles,  particularly  in  the  case  of  the  Be  and  B  implantations. 

Shown  in  Figures  8.1  -  8.4  are  the  theoretical  LSS  Gaussian  distri¬ 
butions  for  the  implants  considered.  The  mean  range,  <x>,  along  the  projected 
range  path  and  the  standard  deviation  along  the  projected  range  path,  a, 
were  obtained  from  a  computer  program  utilizing  the  Gibbons,  Johnson,  and 
Mylroie  formulation  of  the  LSS  theory  [37].  Table  8.2  lists  these  LSS  theo- 
rectical  values  as  <x>„  and  a_.  From  Figure  8.1  -  8.4  it  can  be  seen  that  the 
experimental  modal  depths  are  within  ~  +  10%  of  the  LSS  mean  range  values. 
However,  it  appears  from  Figure  8.1  and  8.2  that  the  LSS  value  of  a  can  differ 
significantly  from  experimental  values.  It  has  been  noted  by  others[24, 35,190] 
that  based  on  experimental  results  the  uncertainty  in  the  LSS  electronic 
stopping  power  is  ~  20-30%  (due  to  approximations  of  linear  energy  loss  with 
velocity  and  in  an  associated  proportionality  constant  [31,35]).  In  particular. 


TABLE  8.2 


Gibbons' 
tion  aG> 


Ion 


mean  range  <x>G>  Winterbon's  mean  range  <x>w»  Gibbons'  standard  devia- 
Winterbon's  standard  deviation  ow>  skewness  kurtosis 


ENERGY  <x>G  <x>w  c w  oG  Yj.  &2 

(keV)  d)  d)  (X)  (£) 
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Schiott  [193]  has  shown  that,  due  to  such  underestimation  of  the  magnitude 
of  electronic  stopping  for  some  projectile-target  combinations,  LSS  can 
overestimate  experimental  values  of  o  by  20-30%.  As  a  result,  the  large 
deviations  of  the  LSS  Gaussian  from  the  experimental  profiles  might  be 

9 

expected  because  of  the  dominant  role  electronic  stopping  plays  in  the  Be  and 

Implants.  Yet,  a  comparison  of  the  theoretical  values  of  a  in  Table  8.2 

to  the  empirical  values  of  a  in  Table  8.3  reveals  the  differences  to  be 

9  10  26  28 

less  in  the  cases  of  Be  and  B  than  in  the  Mg  and  Si  cases.  Therefore, 

the  differences  between  the  LSS  Gaussian  and  the  experimental  profiles  in 

Figures  8.1  and  8.2  do  not  result  from  an  over estimation  of  a.  Rather,  these 

9  10 

differences  in  the  cases  of  the  Be  and  B  implants  seem  more  simply  to  be 

due  to  the  large  values  of  the  third  and  fourth  moments  (Table  8.3).  Such 

large  values  of  Yj^  and  result  in  highly  asymmetric  profiles  which  the 

symmetric  Gaussian  is  incapable  of  adequately  representing. 

8.3  Edgeworth  and  Pearson  I  Distributions 

Motivated  by  these  deficiencies  of  the  LSS  Gaussian  approximation 
9  10 

for  the  Be  and  B  profiles,  we  have  examined  the  fits  obtained  with  Edgeworth 
and  Pearson  distributions  using  four  moments  based  on  theoretical  results. 
Mathematical  definitions  of  these  first  four  moments  of  a  distribution  are 
listed  here  for  completeness  [190]: 


<x>  (mean) 


r 

I  xf(x)  dx 
—  00 


(skewness) 


J  (x-<x>)3  f(x)  dx 


3 

a 


(x-<x>)^ 


f(x)  dx 


$2  (kurtosis) 


4 


a 

A  computer  program  was  written  which  utilizes  the  quantities  listed  in 
Table  8.2  to  generate  the  Edgeworth  and  Pearson  I  distribution  solutions. 

A  comment  concerning  the  definition  of  <x>  and  <x>  in  Table  8.2 

V*  W 

must  be  made  at  this  point.  The  first  moment  as  obtained  from  either  the 
Gibbons  [37]  or  Winterbon  [35]  tables  is  the  average  or  mean  range  of  the 
distribution  along  the  projected  range  path.  This  moment  has  often  been 
referred  to  as  the  "projected  range".  The  experimentally  relevant  quantity, 
however,  is  the  mode  or  peak  of  the  distribution  f(x),  which  occurs  at  the 
point  where  df(x)/dx  *  0.  The  mean  and  modal  values  are  equal  only  when  the 
distribution  is  symmetrical,  as  in  the  case  of  the  Gaussian  [190].  The 
approximate  position  of  the  modal  range  can  be  obtained  from  the  first  four 
moments  of  the  distribution  by  use  of  the  Pearson  approximation  given  in 
reference  194  as  equation  (3.87): 


mode 


<x>  +  a 


where 


a 


°y1^62  + 

10B2  -  12y*  -  18 


oy1(62  +  3) 
A 


The  quantities  a  and  A  are  defined  by  Hofker  et  al.  [190].  As  can  be  seen 
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from  this  equation*  for  distributions  with  negative  values  of  skewness  » 
the  modal  range  is  deeper  than  the  mean  range  and  vice-versa  for  distribu¬ 
tions  with  positive 

The  values  of  <x>  and  o  listed  in  Table  8.2  are  the  mean  range 

G 

and  standard  deviation  along  the  projected  range  path  obtained  from  our 
computer  program  utilizing  the  Gibbons,  Johnson,  and  Mylroie  formulation  of 
the  LSS  theory[37].  The  moments  <x>w»  o w,  y^»  and  B2  listed  in  Table  8.2 
have  been  calculated  from  Winterbon's  tables  of  data  [35]  with  the  assumption 
of  mass  scaling  and  the  substitution  of  the  compound  target  with  a  monatomic 
target  of  average  charge  and  mass. 

8 . 31  Edgeworth  Distribution 

The  Edgeworth  distribution  is  a  series  expansion  in  which  the 
series  terms  are  related  to  successively  higher  order  moments  [36,194].  This 
distribution  may  be  expressed  as  the  product  of  a  Gaussian  function,  g(x), 
and  a  polynomial ,  p (x) : 

f(x)  *  p(x)  g(x) 

with  the  polynomial  terminated  to  include  the  highest  order  moment  desired. 

In  this  work  the  following  expression  for  the  polynomial,  terminated  to 
include  all  moments  up  to  the  fourth,  has  been  used[37]. 

2 

Y  Y 

p(x)  ■  1  +  ~  (x^  -  3x)  +  -yj  (x^  "  15x^  +  45x^  -  15) 

0  -3 

2  J  4  2 

+  ~2^~  (x  ~  6x  +  3) 

The  calculated  Edgeworth  distributions  are  shown  in  Figures  8.1  - 
8.4.  These  particular  curves  have  been  constructed  using  the  moments  <x>q, 
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<jg,  and  $2*  In  all  four  cases  the  curve  is  terminated  at  the  point  where 
the  distribution  function  begins  to  exhibit  oscillations.  Such  oscillations 
have  been  previously  encountered  and  attributed  to  the  termination  of  the 
polynomial  p(x)  [37].  In  all  cases  the  Edgeworth  series  peaks  at  a  deeper 
point  than  the  LSS  Gaussian.  This  type  of  behavior  has  been  previously 
reported  in  conjunction  with  the  Edgeworth  series  [36]  and  should  be  expected 
since  we  are  modifying  a  Gaussian  to  take  into  account  the  predicted  distri¬ 
bution  skewness.  Since  the  theoretical  values  of  in  Table  8.2  are  all 
negative,  the  Edgeworth  modal  range  is  deeper  than  the  mode  of  the  LSS 
Gaussian  which,  due  to  the  Gaussian's  symmetry,  is  also  the  mean  range  <x>. 

9 

In  addition,  because  the  magnitude  of  the  skewness  is  greatest  in  the  Be 
and  implants,  the  largest  difference  between  mean  and  modal  ranges  is 
seen  in  Figures  8.1  and  8.2. 

8.3.2  Pearson  I  Distribution 

The  Pearson  distributions  are  a  family  of  frequency  functions 
based  on  the  solutions  of  the  differential  equation: 

df (x)  m  (x  -  a)  f(x) 

dx  ,  ,  ,  ,  ,  2 

b  +  b.x  +  b„x 
o  1  2 


where  f(x)  is  the  frequency  function  [190,194]. 
b^,  and  b ^  are  [190,194]: 


The  expressions  for  a,  b  , 

o 


aYl^82  + 

A 

c2(4S2  -  3yJ) 


b 


1 


a 
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(2B2  -  3y ‘  -  6) 


where 


A  -  1082  -  12y*  -  18 

In  contrast  to  series  expansion  distributions  such  as  the  Edgeworth,  the 

Pearson  distributions  are  determined  entirely  by  the  first  four  moments. 

the  particular  Pearson  distribution  which  is  the  solution  of  the  differential 

equation  is  determined  by  the  nature  of  the  roots  [194]  of  the  denominator 
2 

(b  +  b,x  +  b_x  ).  The  nature  of  the  moments  <x>  ,  a„,  y, ,  and  were  such 
o  1  2  w  w  1  2 

as  to  yield  real  roots  of  opposite  sign  for  the  characteristic  equation. 

Because  the  Pearson  IV  distribution  is  valid  only  when  the  roots  of 
2 

(bQ  +  bjX  +  b2x  )  are  imaginary,  we  were  excluded  from  using  this  particular 
Pearson  distribution  with  the  theoretical  moments  considered  here.  The  nature 
of  these  roots  demanded  a  Pearson  I  distribution.  The  expression  for  f(x) 
that  we  have  employed  in  our  computer  program  is  given  in  reference  194  as 
equation  (6.12): 


,+m„+l 


(aj+a^ 


B.(mj+1»«2+1) 


(  1  +2-  ) 
al 


where  x  is  measured  from  the  mode.  The  constants  in  this  equation  are  defined 


p  (nij+1)  p  (m2+l) 


B(V1-  Vl>  ■  r  (a^+mpj)- 


1  b2(ai+a2) 


2 


b2(al  +  a2} 


where  a^,  satisfy  the  relations 


al*  a2  >  ° 


bQ  +  b^x  +  b2x  -  b2(x  +  a1)(x  -  a2> 


The  range  of  the  distribution  is  - a ^  <  x  <  a2  where  again  x  is  measured  from 

the  mode.  This  is  in  contrast  to  the  Pearson  IV  distribution  which  has 

unlimited  range  in  both  directions. 

The  calculated  Pearson  I  distributions  are  shown  in  Figures  8.1  - 

8.4.  In  all  four  cases  the  Pearson  I  distribution  yields  a  slightly  improved 

fit  over  the  LSS  Gaussian  approximation  in  the  near  surface  region.  The 

rapid  drop-off  of  the  curves  beyond  the  mode  is  characteristic  of  the 

distribution  as  it  approaches  the  range  validity  limit  [194].  In  the  case  of 
26  28 

the  Mg  and  Si  implants,  the  Pearson  I  mode  coincides  quite  closely  to 

9 

the  experimental  mode.  However,  in  the  profiles  of  primary  concern  (  Be  and 
10B  implants)  the  Pearson  I  is  similar  to  the  Edgeworth  in  that  the  distri¬ 
bution  mode  is  several  hundred  angstroms  deeper  than  the  experimental  mode. 

9  10 

The  experimental  mode  for  the  Be  and  B  implantations  is  in  fact  closer  to 
the  theoretical  range  value  <x>  than  to  the  expected  modal  range,  which  is 
fortuitous  for  the  Gaussian  approximation. 

Therefore,  with  the  exception  of  a  slight  improvement  in  fit  in 
the  near  surface  region,  it  can  be  seen  that  Edgeworth  and  Pearson  I  distri¬ 
butions  constructed  with  four  moments  derived  from  theoretical  values  are 


not  significantly  better  representations  of  implant  profiles  in  InP  than  is 
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the  standard  LSS  Gaussian.  Based  on  these  results  it  appears  that  for 
certain  implants  in  InP  accurate  prediction  of  the  implant  profile  using 
moments  derived  from  theory  (with  the  substitution  of  a  monatomic  target) 
will  be  difficult. 

8.4  Pearson  IV  Distribution  and  Emprically  Determined  Moments 

In  collaboration  with  R.  G.  Wilson  of  Hughes  Research  Laboratory, 
the  first  four  moments  of  the  implants  considered  here  have  been  emprically 
determined  by  iteratively  fitting  the  Pearson  IV  distribution  solution  to  the 
experimental  profile  [191].  The  Pearson  IV  distribution  is  employed  for  these 
empirical  curves,  rather  than  other  four  moment  distributions  such  as  the 
Edgeworth,  because  of  its  established  success  in  describing  implant  profiles 
in  Si  and  in  GaAs  [190,191].  The  moments  obtained  with  this  procedure  are  a, 
Y^,  and  82*  The  first  moment  <x>  is  obtained  by  using  the  Pearson  approxima¬ 
tion  [194]  given  above,  together  with  the  experimental  mode  and  the  other 
empirically  derived  moments.  These  first  four  moments  are  given  in  Table  8.3. 
It  can  be  noted  by  comparing  the  theoretical  four  moments  in  Table  8.2  to  the 
empirical  four  moments  in  Table  8.3  that  the  largest  discrepancies  occur  in 
the  values  of  skewness,  y^,  and  of  kurtosls,  $2*  A  similar  result  concerning 
calculated  and  experimental  moments  has  also  been  reported  by  Hofker  ejt  al. 
[190],  It  is  not  the  intent  here  to  determine  the  cause  of  these  discrepancies 
between  the  theoretical  and  empirical  moments;  however,  Hofker  et  al.[190] 
have  found  that  when  larger  values  of  the  electronic  stopping  coefficient 
are  used  in  their  theoretical  calculations  the  resulting  theoretical  moments 
correspond  reasonably  well  to  the  experimental  moments.  Also  note  that,  as 
expected  with  Pearson  IV  distributions,  the  moments  in  Table  8.3  yield  imag¬ 
inary  roots  for  the  characteristic  equation  mentioned  above.  The  resulting 
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TABLE  8.3 


Empirically  determined  Pearson  IV  moments:  The  mean  range  <x>,  the  standard 
deviation  a,  skewness  ,  and  kurtosis  $2 
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Figure  8.5:  SIMS  experimental  points  for  Be  implant  (100  keV, 
and  empirically  derived  Pearson  IV  distribution  [195] . 
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Pearson  IV  distributions  and  their  fit  to  experimental  points  are  shown  in 
Figures  8.5  -  8.8.  As  with  Si  [190]  and  GaAs  [191],  these  empirically  derived 
Pearson  IV  distributions  excellently  fit  the  implantation  profiles  in  InP 
over  at  least  two  orders  of  magnitude.  Figures  8.5  -  8.8  illustrate  that  if 
the  appropriate  set  of  four  moments  is  employed  the  Pearson  IV  distribution 
is  quite  capable  of  describing  implantation  profiles  in  InP. 

8.5  Summary 

We  have  observed  in  InP  that  light  mass  projectile  implantations 
highly  dominated  by  electronic  stopping  cannot  be  accurately  represented  by 
standard  LSS  Gaussian  profiles.  With  the  intention  of  obtaining  better  pre¬ 
dictions  of  implantation  profiles,  we  have  examined  profiles  obtained  with 
statistical  distributions  which  require  higher  order  moments.  The  moments 
for  these  distributions  have  been  calculated  from  Winterbon's  tables  of 
data  [35],  which  are  based  on  theoretical  considerations.  In  deriving  these 
moments  we  have  substituted  for  the  compound  target  a  monatomic  target  of 
average  charge  and  mass,  and  we  have  assumed  the  validity  of  mass  scaling [35]. 

Using  theoretical  values  for  the  first  four  moments,  profiles 
have  been  constructed  using  the  Edgeworth  and  Pearson  I  distribution  functions. 
The  Pearson  IV  distribution  would  be  preferable  to  the  Pearson  I  because  of 
its  reported  success  in  describing  implantation  profiles  in  Si  and  in  GaAs; 
however,  the  values  of  the  theoretical  moments  in  our  implants  are  such  as  to 
preclude  the  Pearson  IV  and  require  the  Pearson  I  solution  of  the  Pearson 
family  differential  equation.  It  is  found  that  Edgeworth  and  Pearson  I 
distributions  constructed  with  four  moments  derived  from  theoretical  vaHese 
are  not  significant  improvements  over  the  standard  LSS  Gaussian  representa¬ 
tion  of  implant  profiles  in  InP. 
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As  noted,  the  theoretical  values  for  the  first  four  moments  have 

been  obtained  by  substituting  the  compound  target  with  a  monatomic  target 

of  average  charge  and  mass.  Such  an  approximation  may  be  appropriate  for 

GaAs,  in  which  the  ratio  of  target  masses  is  ”  1:1  but  in  InP  this  ratio  is 

9 

~  3.7:1.  In  this  regard  the  backscattering  of  Be  in  InP  should  differ  from 
that  in  the  almost  equivalent  monatomic  target,  Ge.  Winterbon  has  observed 
that  for  4:1  target  mass  ratios  the  moments  derived  with  the  monatomic  target 
approximation  differ  from  those  derived  without  this  assumption  [35].  It 
remains  to  be  determined  whether  more  accurate  profiles  could  be  constructed 
from  the  Edgeworth  or  Pearson  I  distribution  using  theoretical  moments  derived 
without  the  substitution  of  an  equivalent  monatomic  target. 

We  have  empirically  determined  the  four  moments  of  the  experimental 
implantation  profiles.  The  values  of  these  moments  are  such  as  to  permit 
the  use  of  the  Pearson  IV  distribution.  These  empirically  determined  moments 
together  with  the  Pearson  IV  distribution  provide  excellent  fits  to  the 
experimental  profiles  over  at  least  two  orders  of  magnitude.  If  the  appro¬ 
priate  set  of  four  moments  is  obtained,  the  Pearson  IV  distribution  is 
capable  of  accurately  describing  implantation  profiles  in  InP. 

In  light  of  these  results  it  appears  that  £  priori  predictions  of 
implantation  profiles  in  InP  based  on  theoretical  moments  will  not  be  a 
trivial  matter.  At  the  present,  accurate  descriptions  of  implantation  pro¬ 
files  in  InP  will  best  be  obtained  using  the  Pearson  IV  distribution  with 
tables  of  empirically  determined  four  moments. 
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9.  LUMINESCENCE  ASSOCIATED  WITH  DAMAGED 
InP  SURFACES 

9.1  Introduction 

During  the  photoluminescence  studies  discussed  in  Chapter  6,  a 

radiative  transition  was  observed  in  the  PL  spectra  of  damaged  InP  surfaces 

which  was  not  seen  in  the  spectra  of  InP  with  polished  surfaces.  Several 

publications  have  drawn  attention  to  the  unique  surface  properties  of  InP 

in  comparison  with  other  III-V  compounds.  Capacitance-voltage  studies 

on  metal-insulator-semiconductor  (MIS)  n-type  In?  capacitors  report  (almost 

Independently  of  the  dielectric  employed)  interface  state  densities  from 

2  -  5  x  10^  cm  ^  eV  ^  113-15,  196].  This  density  is  lower  than  that  for 

12  -2  -1 

n-type  GaAs  (~  2  x  10  cm  eV  )  and  makes  n-type  InP  an  attractive 

material  for  integrated  logic  circuits  employing  normally-off 

devices  [11,12,197].  Additionally,  several  articles  have  focused  interest 

on  the  surface-related  optical  properties  of  InP.  Of  particular  interest 

is  the  report  of  a  low  surface  recombination  velocity  for  n-type  InP 
3-1  7-1 

of  ~  10  cm  sec  (versus  10  cm  sec  in  GaAs)  and  a  photoluminescence 
(PL)  intensity  ~  100  times  greater  than  that  obtained  from  GaAs  of  similar 
free  electron  concentration  [140].  In  this  chapter  another  unusual  optical 
feature  of  the  InP  surface  is  discussed. 

9.2  Experimental  Results 

Three  types  of  surfaces  have  been  investigated  and  these  are 
shown  in  Figure  9.1.  Figure  9.1(c)  shows  the  surface  of  a  sample  prior  to 
any  polishing.  This  surface  is  the  result  of  the  sawing  operation  used  to 
obtain  (100)  oriented  wafers  from  the  single  crystal  boule.  Such  surfaces 
will  be  referred  to  here  as  "as-sawn".  Figure  9.1(a)  shows  the  results  of 
chemically  polishing  an  "as-sawn"  surface  with  a  1%  bromine-methanol  solution 


Figure  9.1;  Optical  micrographs  of  the  three  surfaces  discussed: 
(a)  "polished",  (b)  "semi-polished",  (c)  "as-sawn"  [211]. 
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and  a  rotating  polishing  pad.  The  results  of  this  operation  are  monitored 
with  an  optical  microscope  and  the  polishing  is  continued  until  scratches, 
digs,  or  pits  are  no  longer  visible  on  the  surface.  The  mirror-like  surface 
of  Figure  9.1(a)  is  referred  to  here  as  "polished".  The  surface  displayed 
in  Figure  9.1(b)  results  from  chemical  polishing  similar  to  that  employed 
for  "polished"  surfaces,  but  in  this  case  the  polishing  is  terminated 
before  a  perfect  mirror  surface  is  obtained.  The  scratches  visible  in 
Figure  9.1(b)  are  the  result  of  incomplete  polishing  of  "as-sawn"  surfaces. 
Surfaces  similar  to  that  in  Figure  9.1(b)  are  referred  to  here  as  "semi- 
polished"  . 

Figure  9.2  displays  the  photoluminescence  spectra  obtained  from  the 
three  surfaces  described  above.  Referring  to  Figure  9.2(a)  the  spectral 
features  characteristic  of  standard  Fe-doped  InP  are  seen.  The  origins  of 
these  features  have  already  been  discussed  in  Section  6.2. 

An  expected  result  of  this  work  is  the  observation  of  a  decrease  in 
luminescence  intensity  when  the  surface  is  less  than  perfectly  polished. 

The  effects  of  polishing  damage  on  the  PL  intensity  of  GaAs  have  been 
reported  in  the  past  [198-200].  These  studies  indicate  that  surface  damage 
can  decrease  the  radiative  efficiency  by  at  least  an  order  of  magnitude  from 
that  obtained  from  the  best  surface.  A  similar  result  is  observed  here  for 
InP  surfaces.  Comparing  the  integrated  intensitites  of  the  "semi-polished" 
and  "as-sawn"  spectra  to  that  of  the  "polished"  spectra  we  obtain  20%  and 
40%  respectively.  Presumably  such  decreases  in  luminosity  can  be  attributed 
to  surface  damage  on  a  macroscopic  scale  (e.g.,  scratches)  and  on  a  micro¬ 
scopic  scale  (e.g. ,  dislocations,  vacancy-complexes)  which  introduce  non- 
radiative  recombination  channels  [201].  The  difference  in  luminosity  be¬ 
tween  "semi-polished"  and  "as-sawn"  surfaces  may  appear  contradictory; 
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Figure  9.2:  Typical  photo luminescence  spectra  from  the  three  surfaces  shown 
in  Figure  9.1.  The  sample  temperature  for  these  measurements  was  5°K  [211]. 
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however,  the  majority  of  this  effect  can  be  explained  in  terms  of  the  latter 
surface  having  a  much  larger  absorption  of  the  incident  radiation.  Crude 
measurements  of  the  reflected  power  from  the  surfaces  of  Figure  9.1  have  been 
made.  "Polished"  and  "semi-polished"  surfaces  yield  comparable  reflection 
while  the  reflected  power  from  "as-sawn"  surfaces  is  ~  20-40  times  weaker. 
These  measurements,  together  with  the  equation  for  luminescent  intensity 
given  in  reference  140,  suggest  that  absorption  effects  are  the  predominant 
cause  of  the  observed  luminosity  differences  seen  in  Figures  9.2(b)  and 
9.2(c).  Several  publications  have  shown  that  the  PL  intensity  from  InP  is 
strongly  affected  by  surface  band-bending  effects  leading  to  a  space  charge 
layer  and  an  optically  dead  layer  [202-204].  Differences  in  surface  band¬ 
bending  between  "as-sawn"  and  "semi-polished"  surfaces  may  also  play  a 
role  in  the  observed  luminosity  differences. 

An  tuusual  result  of  this  work  is  the  appearance  of  a  new 
spectral  feature  at  1.30(6)  eV.  As  can  be  seen  from  Figure  9.2(a)  "polished" 
surfaces  exhibit  no  PL  peaks  in  this  energy  range,  whereas  spectra  obtained 
from  "semi-polished"  surfaces  exhibit  a  weak  luminescence  in  this  region. 
Greater  damage  to  the  surface  in  the  "as-sawn"  case  results  in  the  strongest 
1.30(6)  eV  feature.  Figure  9.2  suggests  that  the  evolution  of  this 
transition  is  related  to  the  degree  of  mechanical  surface  damage.  This 
transition  has  also  been  observed  under  other  circumstances  which  further 
support  this  view  of  its  origin.  "Polished"  surfaces  degraded  by  exposure 
to  concentrated  HC1  exhibit  weak  spectral  features  in  the  1.30(6)  eV  region 
which  were  not  present  in  the  spectra  prior  to  the  HC1  treatment.  ,  This 
1.30(6)  eV  feature  has  also  been  observed  when  the  excitation  beam  is 
incident  on  an  is ->lated  pit  of  an  otherwise  perfect  "polished"  surface. 

These  results  support  the  surface  damage  origin  of  this  transition  proposed 


by  Street  e£  al.  who  were  to  our  knowledge  the  first  to  observe  this 
feature  [204] .  They  observed  this  PL  band  from  samples  which  cleaved  poorly 
under  vacuum  and  resulted  in  surfaces  with  a  very  large  density  of  surface 
steps.  It  should  be  pointed  out  that  our  material  and  processing  techniques 
differ  in  several  respects:  (1)  Street  et  aL.  performed  experiments  on 
(110)  InP  while  ours  were  on  (100)  In£;  (2)  their  crystals  were  lightly 
doped  n  and  p-type  while  ours  are  semi- insulating;  (3)  their  samples  were 
cleaved  and  examined  in  ultrahigh  vacuum  while  ours  have  been  exposed  to 
air.  Despite  these  differences  we  observe  the  identical  transition  reported 
by  Street  et  al. 

In  a  publication  following  their  initial  report  of  the  1.3  eV 
transition.  Street  and  Williams  [205]  reported  observing  a  decay  in  the 
radiative  luminescence  of  the  transition  as  a  result  of  room  temperature 
annealing  in  a  period  of  days.  They  stated  that  this  may  explain  why  this 
feature  is  not  routinely  observed  in  the  PL  spectra  of  InP.  Although  we 
made  no  effort  to  examine  such  annealing  effects,  we  have  observed  the 
transition  in  "as-sawn"  surfaces  months  after  the  sawing  operation.  This 
may  be  the  result  of  the  much  greater  damage  of  these  surfaces  in  comparison 
to  those  examined  by  Street  and  Williams. 

9.3  Summary 

Surface  studies  performed  in  the  last  few  years  indicate  that  no 
intrinsic  surface  states  exist  in  the  fundamental  gap  of  vacuum  cleaved 
GaAs  or  InP  (110)  surfaces  [206,207].  It  is  now  believed  that  previous 
reports  of  Fermi  level  pinning  were  due  to  extrinsic  states.  In  light  of 
these  findings  and  considering  the  extent  of  damage  to  our  surfaces  (together 
with  their  exposure  to  air) ,  it  is  believed  that  transitions  due  to  intrinsic 
surface  states  can  be  ruled  out.  Using  a  value  of  1.423  eV  for  the  baudgap 


energy  at  5°K  1139]  would  place  the  observed  transition  at  117  meV  from 
the  bulk  conduction  band  edge;  however,  the  exact  location  of  this  level 
with  respect  to  the  surface  conduction  band  is  unknown  due  to  probable 
band  bending  at  the  surface.  In  this  regard  the  1.30(6)  eV  feature  may  be 
a  transition  to  an  effectively  shallow  level  (bent  band  edge  to  surface 
defect  state)  or  a  117  meV  transition  to  an  optically  deep  level.  The 
detailed  nature  of  the  defects  giving  rise  to  the  observed  transition  is 
not  known.  In  this  regard  the  influence  of  the  anion  should  not  be  over¬ 
looked,  since  it  appears  to  affect  Schottky  barrier  heights  and  the  surface 
recombination  velocity  in  other  multicomponent  semiconductors  [208-210]. 
Experimental  evidence  suggests  that  chemical  adsorption  of  oxygen  on  InP 
directly  affects  the  P  surface  atoms  but  not  the  In  surface  atoms  [206]. 

The  chemical  volatility  of  the  phosphorus  component  of  InP  would  appear 
to  make  phosphorus-related  defects  likely  candidates  for  understanding  the 
detailed  nature  of  the  1.30(6)  eV  transition. 


10.  BORON  CONTAMINATION  AND  PRECIPITATION 
IN  LEC  GROWN  InP 

In  this  chapter,  polycrystalline  InP  and  single  crystal,  liquid 
encapsulated  Czochralski  (LEC)  grown  InP  have  been  examined  for  evidence  of 
unintentional  boron  contamination  due  to  growth  procedures.  Precipitates 
of  boron  or  a  boron  compound  have  been  found  in  InP  grown  by  the  LEC  method 
with  ooric  oxide  encaPsulation  and  pyrolitic  boron  nitride  (pBN) 

crucibles. 

10.1  LEC  Growth  Technique 

The  single-crystal  InP  substrates  analyzed  here  are  the  result 
of  two  separate  growth  steps,  each  of  which  may  introduce  unintentional 
impurities.  The  first  growth  procedure  is  the  synthesis  or  compounding  of 
indium  and  phosphorus  to  obtain  the  polycrystalline  InP  charge  used  in  the 
subsequent  single  crystal  growth  [20,212,213].  During  synthesis  a  boat 
containing  an  In-rich  melt  is  located  within  a  furnace  zone  at  ~  1060- 
1100°C  And  is  in  contact  with  phosphorus  vapor  from  solid  phosphorus 
situated  in  a  cooler  furnace  zone  (~  400-500°C) .  Synthesis  is  accomplished 
by  slowly  pulling  the  In-rich  melt  through  a  temperature  profile  leading 
to  directional  solidification  of  the  melt. 

The  second  step  involves  the  growth  of  single  crystal  InP  from 
the  synthesized  polycrystalline  ingot.  Crystalline  InP  has  been  prepared 
by  a  number  of  techniques  including  zc.ie  melting,  gradient  freeze,  horizontal 
Bridgman,  and  liquid  encapsulated  Czochralski  (LEC)  pulling  [16,20,117,212-215] 
As  a  result  of  the  large,  round  boules  (~  100-300  g,  ~  1.5"  diameter) 
obtainable  by  the  LEC  method,  this  presently  appears  to  be  the  most  practical 
technique  for  growing  device-grade  InP.  In  the  LEC  technique  the  crucible 
containing  the  melt  is  housed  in  a  high  pressure  chamber  which  has  been 


evacuated  and  backfilled  with  an  inert  gas  (e.g.,  N2  or  Ar)  [216,217]. 

A  liquid  encapsulant  (usually  boric  oxide,  B^^)  covers  the  melt  in  the 
crucible.  The  encapsulant  ideally  serves  as  a  liquid  seal  preventing  vapor 
loss  of  volatile  elements  if  the  pressure  of  the  inert  gas  above  the 
encapsulant  is  greater  than  the  dissociation  pressure  of  the  melt.  As 
the  crystal  emerges  from  the  melt  and  the  liquid  encapsulant  layer,  a 
film  of  encapsulant  adheres  to  the  outside  of  the  boule,  thereby  preventing 
loss  of  volatile  components  from  the  growing  crystal.  Upon  completion  of 
pulling,  this  B2°3  ^ayer  *s  removed  with  hot  water. 

The  possible  sources  of  contamination  in  this  growth  process  are: 
(i)  the  boat  containing  the  In  melt  during  synthesis;  (ii)  the  encapsulating 
layer  used  during  LGC  pulling;  (ill)  the  crucible  containing  the  melt 
during  LEC  pulling.  B203  is  used  almost  universally  as  the  liquid  encapsulant 
in  LEC  growth  of  various  III-V  compounds  because  of  its  chemical  stability, 
density,  low  melting  point,  viscosity,  and  transparency  [216-218].  The 
composition  of  the  synthesis  boat  and  of  the  melt  crucible  varies  from  lab 
to  lab,  but  the  most  common  materials  are  silica  or  pyrolytic  boron  nitride 
(pBN) .  It  has  been  reported  that  silica  crucibles  introduce  silicon 
impurities  in  LEC  crystals  due  to  dissolution  of  the  crucible  by 
B2<>3  [185,219-221].  Additionally,  silica  crucibles  containing  tend 

to  fracture  when  cooled  to  room  temperature  [20,217].  To  avoid  these  problems 
the  use  of  pBN  crucibles  has  become  common. 

10.2  Experimental  Results 

Tables  10.1  and  10.2  summarize  our  study  of  boron  contamination 


in  the  growth  of  InP.  Impurity  concentrations  were  measured  by  secondary- 
ion  mass  spectromety  (SIMS)  using  a  CAMECA  IMS-3f  ion  microprobe. 


Table  10.1 


SAMPLE 


SYNTHESIS  LEC 

BOAT  CRUCIBLE 


“B  CONCENTRATIONS 
(atoms/cnr) 


PRECIPITATION 


Metals  Res.  Corp 

Boule  -  ? 

Slice  -  ? 

Si02 

Si02 

(3.2  +  1.6) 

x  101 

Crys. 

Boule 

Slice 

Comm.  Inc. 

-  136  (CCIPS) 

-  Top 

sio2 

Si02 

(11.0  +  4.7) 

x  10 

Crys. 

Boule 

Slice 

Comm.  Inc. 

-  136  (CCIPS) 

-  Bottom 

sio2 

Si02 

(5.3  +  3.9) 

x  101 

NRL 

Boule 

Slice 

-  2-25 (H) 

-  1 

Si02 

pBN 

(3.8  +  2.1) 

x  101 

NRL 

Boule 

Slice 

-  2- 25(H) 

-  55 

Si02 

pBN 

~  3  x  1018 

NRL 

Boule 

Slice 

-  2-10 (H) 

-  Bottom  1/3 

pBN 

pBN 

~  1  x  1018 

The  results  of  our  analysis  of  various  NRL  (Naval  Research 


Laboratory)  polycrystalline  samples  are  shown  in  Table  10.1.  Our  results 
are  in  general  agreement  with  earlier  reports  of  boron  concentrations 
measured  in  polycrystalline  ingots  from  other  labs  [17,185].  We  have 
examined  samples  from  both  ends  of  polycrystalline  ingots  synthesized  in 
silica  and  in  pBN  boats.  No  variation  of  boron  concentration  along  the 
length  of  a  polycrystalline  ingot  is  seen  in  ingots  grown  with  either  silica 
or  pBN  boats.  Additionally,  from  Table  10.1,  no  significant  differences  in  the 
B  levels  can  be  detected  between  silica  boat  ingots  and  ingots  synthesized 
in  pBN.  Our  data  indicate  that  within  detection  limits,  pBN  boats  have 
little  influence  on  boron  contamination  of  the  melt. 

Several  single  crystals  of  semi- insulating  InP  were  examined,  and 
the  results  are  summarized  in  Table  10.2.  Except  in  samples  where  precipitation 
was  detected,  the  measured  boron  concentrations  are  comparable  to  those 
reported  by  other  groups  for  LEC  grown  InP  [185,216,217,222,223].  Three 
samples  shown  in  Table  10.2  were  synthesized  and  pulled  exclusively  from 
silica  ware,  and  employed  a  encaPsulaat  during  pulling.  The  average 

level  of  boron  present  in  these  samples  is  slightly  higher  than  that  found 
in  the  polycrystalline  samples  of  Table  10.1,  however,  within  the  experi¬ 
mental  uncertainty  of  the  SIMS  technique,  the  concentrations  still  overlap. 

From  these  results  we  conclude  that,  at  most,  the  employed  in  LEC 

pulling  contributes  only  a  small  amount  of  boron  to  the  crystal.  Hopkins 
et  al.  [224]  have  recently  reported  the  same  result  for  LEC  grown  GaAs. 

This  3tudy  was  initiated  because  of  the  discovery  of  boron-rich 
regions  or  precipitates  while  depth  profiling  boron  implants  in  single 
crystals  from  NRL  boule  2-10(H).  In  the  course  of  depth  profiling  to  2.5  um. 


Sputtering  Time  (min) 


LP-1739 


Figure  10.1:  concentration  SIMS  profiles  for  slices  from  top  (slice  1) 

and  bottom  (slice  55)  of  NRL  boule  2-25(H)  [178]. 


approximately  2-3  precipitates  of  varying  size  would  be  encountered.  Between 
precipitates  the  boron  background  in  these  samples  is  comparable  to  samples 
listed  in  Table  10.2  which  did  not  contain  precipitates.  Boule  2-10(H) 
was  synthesized  and  pulled  entirely  in  pBN  growth  ware  with  the  analyzed 
slices  taken  from  the  botton  (last-to-f reeze)  third  of  the  boule. 

Two  slices  of  semi-insulating  InP  from  different  portions  of  the 
same  NRL  boule  (2-25 (H))  were  analyzed.  To  isolate  possible  sources  of 
boron  contamination,  the  polycrystalline  InP  charge  for  this  boule  was 
synthesized  in  silica.  Subsequently, the  single  crystal  ingot  was  pulled 
from  a  pBN  crucible.  Comparison  of  the  SIMS  B  profiles  for  these  two 
slices  is  shown  in  Figure  10.1.  Slice  1,  which  is  from  the  top  (first-to- 
freeze)  end  of  the  boule,  has  a  boron  concentration  comparable  to  other 
nc  -precipitate  containing  samples  as  seen  from  Table  10.2.  No  large 
precipitates  are  observed  in  slice  1.  Slice  55  is  taken  from  the  bottom  of 
boule  2-25(H).  As  seen  from  Figure  10.1  the  boron  concentration  in  this  sample 
differs  drastically  both  in  nature  and  concentration  from  that  of  slice  1. 
Whereas  slice  1  has  an  essentially  uniform  level  of  boron  contamination, 
slice  55  contains  large,  distinct  precipitates  with  boron  concentrations 
2-3  orders  of  magnitude  higher  than  the  boron  level  in  slice  1.  This  pre¬ 
cipitation  phenomenon  can  be  more  clearly  observed  in  a  series  of  photo¬ 
graphs  (Figure  10.2)  taken  at  successive  depths  using  the  ^  ion  image 
from  the  CAMECA  IMS-3f.  In  these  photographs,  clusters  containing  high 
concentrations  of  boron  appear  and  are  then  replaced  by  other  clusters  as 

the  sputtering  proceeds.  A  first  order  estimate  of  the  precipitate 

18  —3 

density  indicates  ~  1  x  10  precipitates  cm  in  slice  55.  Averaging 

the  boron  signal  over  the  precipitate  fluctuations  gives  a  boron  concentra- 
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tion  of  ~  3  x  10  cm  ,  in  sharp  contrast  with  the  ~  4  x  10  cm  observed 


T 


for  slice  1  from  the  same  boule.  we  have  observed  in  various  samples  that 
the  number  of  precipitates  encountered  during  a  profile  depends  on  the 
sample's  location  in  the  boule  with  the  density  of  precipitates  increasing 
towards  the  bottom  (last-to-f reeze)  end  of  the  boule. 

The  exact  chemical  composition  of  the  observed  boron  clusters  is 

27  28  52  56 

not  known.  Profiles  have  been  taken  monitoring  Al,  Si,  Cr,  and  Fe 
simultaneoulsy  with  ^B.  No  correlation  of  boron  precipitates  with  any  of 
these  impurities  has  been  observed.  Components  of  the  precipitates  such  as 
■^N,  160,  or  still  remain  possible  due  to  poor  detection  capability  and 
background  masking  effects. 

It  has  been  reported  by  Mullin  £t  al^.  [17]  that  a  peak  at  1.32 
eV  is  found  in  the  photoluminescence  spectra  of  pulled  crystals  only  and  never 
in  the  spectra  of  polycrystalline  starting  material.  They  noted  that  this 
peak  is  clearly  associated  with  liquid  encapsulated  growth.  It  was  not 

determined  if  the  peak  is  directly  due  to  elemental  boron  or  due  to  an 
interaction  of  B2®3  other  contaminants  [17,222].  We  have  analyzed 

the  photoluminescence  spectra  of  the  sample  shown  in  Tables  10.1  and  10.2. 

With  the  samples  held  at  5°K,  spectra  were  obtained  for  the  energy  range 
from  1.42  -  1.15  eV  using  a  system  resolution  of  1  meV.  No  strong  peak  in 
the  region  about  1.32  eV  is  detected  in  either  slice  55  or  slice  1,  and 
their  spectra  in  this  region  are  comparable  to  those  obtained  from  the 
polycrystalline  samples  of  Table  10.1.  Comparison  of  the  spectra  from  slice 
55  to  that  from  slice  1  yielded  no  significant  differences.  It  appears  that 
boron  does  not  produce  any  luminescence  lines  in  the  range  from  1.42  -  1.15 


eV,  supporting  a  similar  observation  by  Williams  et  al.  [102]. 
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10. 3  Discussion 

Our  work  to  date  has  not  unambigously  established  the  source  of 

boron  contamination  and  precipitation  in  LEC  grown  InP.  We  have  observed, 

however,  that  pBN  in  the  absence  of  ^2^3  ^oes  not  cause  detectable  boron 

contamination  above  what  is  normally  measured  in  samples  grown  in  silica. 

The  polycrystalline  ingots  analyzed  here  were  each  grown  over  approximately 

two  week  periods  while  the  NRL  single  crystals  were  grown  in  1-2  days.  The 

temperatures  involved  for  the  synthesis  and  LEC  growth  of  the  analyzed 

samples  are  comparable.  LEC  growth  proceeds  at  the  melting  point  (1062°C) 

while  the  synthesis  boat  for  the  NRL  polycrystalline  samples  sits  in  a 

o 

linearly  ramped  temperature  zone  of  1060-1100  C.  If  pBN  were  the  source  of 
boron  contamination  then  taking  into  account  the  times  and  temperatures 
involved,  one  would  conclude  that  boron  contamination  similar  to  that  seen 
in  pBN  crucible  grown  single  crystals  should  certainly  be  present  in  pBN 
boat  synthesized  ingots.  Our  results  indicate,  however,  that  the  boron 
levels  at  either  end  of  pBN  synthesized  polycrystalline  ingots  are  com¬ 
parable  to  those  synthesized  in  silica  boats.  On  the  other  hand,  we  observe 
that  in  the  absence  of  a  pBN  crucible  does  not  lead  to  precipitation. 

Although  the  average  level  of  boron  present  in  silica  crucible  grown 
samples  is  slightly  higher  than  that  found  in  the  polycrystalline  samples, 
the  concentrations  are  clustered  in  a  range  of  less  than  an  ordet  of  mag¬ 
nitude.  No  precipitates  have  been  observed  in  ^2^3*  silica  crucible  grown 
single  crystals.  These  results  lead  us  to  postulate  that  the  source  of 
boron  contamination  most  consistent  with  our  observations  is  the  result  of 
interaction  between  B2O2  and  pBN. 


To  study  the  diffusivity  of  boron  in  InP  we  have  examined  the  boron 

atomic  profiles  resulting  from  thirty  minute  of  implanted  InP  encapsulated 

with  2300  X  of  Si„N. .  Very  little  diffusion  of  boron  could  be  detected  in 
3  4 

material  annealed  at  750°C.  Computer  simulation  of  profiles  resulting  from 
various  diffusion  coefficients  has  allowed  us  to  determine  the  upper  limit 
for  the  boron  coefficient  in  InP  to  be  D(750°C)  <  1  x  10  cm  sec 

To  our  knowledge  no  literature  exists  on  the  electrical  properties 
of  boron  in  InP,  although  it  is  commonly  assumed  to  be  electrically  in¬ 
active  since  it  is  isoelectronic  with  In  [216].  Aside  from  the  question  of 
elemental  boron's  electrical  activity  in  InP,  the  boron  precipitates  can  be 
expected  to  cause  problems  in  devices  fabricated  on  such  substrates.  Lat¬ 
tice  distortion  in  precipitate  regions  could  affect  epitaxial  growth  and  could 
also  change  the  "effective"  dielectric  constant  of  the  material  which,  upon 
application  of  an  external  field,  would  lead  to  local  inhomogeneities  of 
electric  field.  Such  precipitates  may  also  create  localized  space-charge 
regions  which  will  influence  carrier  mobility  [225].  These  effects  will 
have  serious  consequences  on  the  performance  of  devices  fabricated  in  InP 
substrates  rich  in  boron  precipitates. 

10.4  Summary 

We  have  discovered  precipitates  of  boron  or  a  boron  compound  in 
single  crystal  InP  grown  by  the  LEC  method  with  an<*  P®^  crucibles. 

The  den*,  -y  of  precipitates  depends  on  the  sample's  location  in  the  boule, 
with  the  number  of  precipitates  increasing  towards  che  last-to-f reeze  end 
of  the  boule.  The  boron  background  concentrations  in  poly crystal line  InP 
and  in  encapsulated,  silica  crucible,  LEC  grown  InP  single  crystals 


lie  in  the  range  from  10*^  -  10*^  cm  We  have  observed  that  pBN  in  the 
absence  of  ^2^3  does  not  lead  to  detectable  boron  contamination  above  what 
is  normally  seen  in  samples  grown  in  silica.  There  are  indications  that 

may  be  responsible  for  possibly  higher  boron  concentrations  than  those 
normally  found  in  polycrystalline  samples,  but  the  increase  is  less  than  an 
order  of  magnitude.  in  the  absence  of  a  pBN  crucible  has  not  been 

observed  to  lead  to  boron  precipitation.  We  postulate  that  the  likely  source 
of  boron  contamination  in  crystals  with  boron  precipitates  is  the  result  of 
an  interaction  of  and  pBN. 

Photoluminescence  spectra  were  obtained  for  the  SIMS  analyzed 
polycrystalline  and  single  crystal  samples.  On  the  basis  of  photo luminescence 
spectra  it  was  not  possible  to  distinguish  between  single  crystals  with 
precipitates  and  those  without.  Boron  does  not  appear  to  produce  any  lumines¬ 
cence  in  the  energy  range  of  1.42  -  1.15  eV. 

On  the  basis  of  boron  atomic  profiles  resulting  from  thirty 
minute,  750°C  anneals  of  *^B  implanted  InP,  a  diffusion  coefficient  of 
D(750°C)  <  10  ^  cm^  sec  ^  has  been  determined  for  boron  in  InP. 


11.  SUMMARY 


If  Ion  implantation  is  to  become  a  standard  doping  technique  in  InP 
it  will  be  necessary  to  develop  reliable  annealing  procedures.  Yet  because 
of  the  relatively  limited  research  effort  expended  to  date  on  implantation 
studies  of  InP  a  large  number  of  basic  issues  remain  to  be  examined.  Three 
such  topics  have  been  addressed  in  this  work.  These  are:  Which  of  the 
commonly  employed  encapsulants  (SiO^,  PSG,  or  Si^N^)  but  suppresses  In  or  P 
outdif fusion  from  the  InP  surface  during  high  temperature  anneals?  What 
are  the  anneal  time  and  temperature  required  for  lattice  recovery  from  im¬ 
plant  damage  and  for  activation  of  implanted  impurities?  What  are  the 
redistribution  characteristics  of  implanted  and  other  intentionally  present 
impurities  during  such  anneals? 

Encapsulation  studies  described  here  Indicate  that  a  PSG  cap  with 
'  5  at.Z  P  is  superior  to  either  a  SiO^  cap  or  a  Si^N^  cap  for  annealing  InP 
in  the  temperature  range  from  450*-800°C.  Using  photoluminescence,  P  vacancy 
related  luminescence  has  been  observed  in  the  spectra  of  SiO^  capped  and 
annealed  InP,  indicating  P  outdif fusion  through  SiO^  encapsulants  of  InP. 

SIMS  analysis  indicates  that  600eC,  60  min  anneals  with  a  SiO^  cap  result  in 
Si  levels  in  the  InP  above  those  found  in  unprocessed  InP  to  depths  of  ~  1  pm 
In  addition,  the  usefulness  of  SiO^  as  an  InP  cap  is  limited  since  it  has 
been  found  to  crack  when  annealed  in  the  temperature  range  from  650-750°C. 

It  has  been  concluded  from  Auger  studies  that  factors  other  than  thermal 
expansion  mismatch  (e.g.,  chemical  Interactions)  play  a  role  in  this  SlO^ 
cracking  phenomenon.  On  the  basis  of  PL  measurements  it  appears  that  little, 
if  any,  P  outdif fusion  occurs  through  a  Si.N,  cap.  It  has  been  observed  with 
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PL  that  high  temperature  anneals  with  a  Si^N^  cap  result  in  luminosity 

levels  greater  than  those  associated  with  unprocessed  material  and  in  the 

emergence  of  a  new  spectral  feature  of  1.378  eV.  Indium  has  been  detected 

by  both  AES  and  SIMS  on  the  surfaces  of  Si^N^  caps  following  30  min  and  60 

min,  750*C  anneals.  SIMS  measurements  of  samples  annealed  for  60  min  at  750°C 

with  SijN^  indicate  Si  levels  above  those  found  in  unprocessed  InP  to  depths 

of  ~  0.75  pm.  No  such  Si  contamination  of  InP  has  been  observed  following 

anneals  with  a  PSG  cap.  PL  measurements  indicate  that  in  comparison  to 

SIO2  or  Si^N^  caps,  PSG  encapsulation  least  affects  the  luminosity  levels  or 

spectral  features  of  the  capped  InP.  Some  SIMS  results  suggest  that  In  is 

outdif fusing  through  the  PSG  caps,  but  further  work  is  required  on  this  point 

before  reaching  a  definite  conclusion.  On  the  basis  of  the  results  discussed 

here  it  is  concluded  that  PSG  encapsulation  with  ~  5  at.Z  P  is  better  than 

either  Si02  or  in  preserving  the  characteristics  of  encapsulated  InP 

during  furnace  anneals  in  the  range  from  450-750°C. 

PL  measurements  of  the  integrated  bandege  luminescence  from 

13  15  -2  4 

annealed  samples  following  51  keV,  10  -  10  cm  Implants  of  He  indicate 

that  30  min  anneals  at  750°C  are  sufficient  for  the  lattice  to  recover  from 

Implant  damage  and  to  restore  pre-implant  levels  of  optical  quality. 

Identical  conclusions  are  obtained  from  PL  examination  of  the  integrated 

13 

bandedge  luminescence  of  annealed  samples  implanted  with  100  keV,  10 
15  -2  9 

10  cm  Be.  From  these  studies  it  is  estimated  that  the  Be  acceptor 
ionization  energy  in  InP  is  41.3  ±  3  meV.  Maximum  Be  optical  activation 
occurs  for  30-60  min  anneals  at  750°C. 
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Using  SIMS,  depletion  zones  of  the  compensating  impurity,  Cr  or 

Fe,  have  been  observed  In  Si^N^  capped  and  800°C  annealed  SI  InP.  Cr 

exhibits  larger  concentration  gradients  in  these  depletion  zones  than  does 

Fe.  Under  typical  annealing  conditions  for  InP  (T  ^  700°C,  t  *  15-30  min) 

9 

it  is  found  from  SIMS  analysis  that  implanted  Be  is  a  rapid  diffusant  in 

15  -2  9 

SI  InP.  High  dose  (10  cm  ),  100  keV  implants  of  Be  into  Si  InP  result 

in  structured  and  correlated  Be  and  Cr  or  Fe  atomic  profiles.  Flat  tails 
9 

of  Be  extending  over  distances  of  a  few  microns  are  observed  in  high  dose 

9 

implanted  SI  InP  but  not  in  comparably  processed  VPE  InP.  Also,  Be  is 

found  to  diffuse  deeper  (~  3.5  ym)  at  annealing  temperatures  of  550°C  than 

15  -2 

in  750°C  anneals  (~  2.3  ym)  of  10  cm  implanted  SI  InP.  SIMS  studies  of 
15  -2 

10  cm  ,  250  keV  Mg  implanted  SI  InP  show  drastic  redistribution  of  Mg  and 

Fe  during  30  min  anneals  for  temperatures  from  550-750°C.  As  in  the  Be 

implant  studies,  the  redistribution  of  Fe  is  correlated  with  the  diffusion 

9  26 

profile  of  the  implanted  Mg.  In  contrast  to  Be  or  Mg,  little  redistri- 
28 

bution  of  Si  is  seen  following  750°C,  30  min  anneals  of  SI  InP  implanted 
with  240  keV,  5  x  1014  cm"2  28Si. 

Although  the  present  results  answer  the  basic  questions  raised 
here,  important  new  questions  have  arisen  regarding  InP  implantation.  Most 
importantly,  further  work  is  required  to  understand  the  physics  governing 
the  redistribution  behavior  of  implanted  impurities  in  SI  InP.  The  behavior 
of  Be  and  Mg  seen  here  seems  to  parallel  the  rapid  diffusion  shown  by  Zn  and 
Cd  during  diffusion  doping  of  InP  [155,157].  Furthermore,  others  have 
pointed  out  that  donor  implants  result  in  higher  electrical  activations  than 
acceptor  lmpants  in  InP  [1841 .  Is  the  redistribution  behavior  of  acceptors 
in  InP  a  fundamental  limitation  leading  to  low  activations?  This  question 
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is  important  not  only  from  a  scientific  viewpoint  but  also  technologically 
since  it  will  severely  limit  the  applicability  of  acceptor  implantation  in 
InP.  Furthermore,  the  frequency  of  unusual  or  unexpected  results  in  this 
work  (i.e..  Chapters  8-10)  illustrates  the  current  limited  understanding 
of  InP  properties.  In  conclusion,  the  results  presented  here  clearly 
indicate  the  need  for  further  studies  before  ion  implantation  of  InP  is  a 
well  understood  technique. 
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